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FOREWORD 


The  US  Army  Electronic  Proving  Ground  was  responsible  for  the  execution  of 
this  project.  Mr.  C.  M.  Giorgi  was  the  Activity  Project  Officer  with  assist¬ 
ance  from  Major  W.  M.  Decker,  IV.  Technical  guidance  as  well  as  the  idea  for 
the  project  was  provided  by  the  Test  and  Evaluation  Command  Project  Officer, 
Dr.  N.  E.  Pentz. 

Dr.  E.  C.  Crittenden  and  Dr.  E.  A.  Milne,  of  the  Naval  Postgraduate  School 
were  responsible  for  conducting  the  specific  technical  investigations  with 
assistance  from  Mr.  S.  R.  Crager  and  Mr.  J.  H.  Conner,  Jr.,  of  that 
organization.  The  material  presented  in  this  report  are  the  results  of  those 
investigative  efforts. 
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SECTION  1 .  SUMMARY 


1 . 1  BACKGROUND 

a.  The  spot  patterns  produced  by  laser  designators  suffer  from  a  number 
of  defects.  The  patterns  are  broadened  and  exhibit  wander  and  intensity  fluc¬ 
tuation  because  of  turbulence  in  the  atmosphere.  Designator  spot  patterns 
also  exhibit  similar  defects  because  of  internal  laser  instability  and  motion 
of  the  transmitter  platform.  In  order  to  correct  the  latter  defects,  a  system 
is  needed  that  is  capable  of  separating  the  contributions  of  the  several 
effects. 


b.  Previous  work  at  the  Naval  Postgraduate  School  (ref  1,  2,  3,  app  E) 
and  elsewhere  (ref  4,  app  E)  has  demonstrated  that  the  effects  of  atmospheric 
turbulence  on  designator  spot  patterns  can  be  expressed  in  terms  of  the 
Optical  Transfer  Function  (OTF)  of  the  atmosphere.  This  quantity  is  the 
Fourier  transform  of  the  spot  profile  due  to  the  atmosphere  alone.  If  other 
causes  of  spot  broadening  are  present,  such  as  laser  instability  and  platform 
motion,  the  spot  profile  will  be  the  convolution  of  the  spot  profiles  due  to 
each  of  the  effects  s.  The  spot  profile  broadening  due  to  laser  and  platform 
instability  can  then  be  separated  from  the  atmospheric  effects  by  dividing  the 
Fourier  transform  of  the  spot  profile  by  the  OTF  of  the  atmosphere  (point-by- 
point  as  a  function  of  spatial  frequency).  The  necessary  OTF  of  the  atmos¬ 
phere  can  be  obtained  by  means  of  a  slit-scanning  telescope  system,  which 
views  a  point  laser  in  the  target  vicinity  from  a  location  near  the  designator 
transmitter.  An  additional  imaging  system,  viewing  the  designator  spot  on  a 
target  screen,  from  a  location  near  the  target,  yields  the  composite  spot 
profile.  Data  from  this  imager,  together  with  that  from  the  telescope  imager 
at  the  transmitter  site,  can  separate  the  various  effects.  As  will  be 
described  later  in  more  detail,  spot  motion  and  spot  broadening  are  separated 
by  computer  tracking  (image  centering)  techniques,  applied  to  the  data  from 
the  spot  profile  TV  camera. 

c.  Because  designator  lasers  use  short  pulses,  the  previously  developed 
mechanical  slit-scan  techniques  cannot  be  used.  Storage  imaging  systems, 
using  TV,  CCD,  or  CID  techniques,  are  needed  to  sense  the  image  for  the  very 
short  pulse  periods,  and  provide  the  equivalent  of  a  slit-scan  during  the 
interpulse  period.  TV  was  chosen  over  CCD  or  CID  techniques  for  this  purpose 
because  of  higher  interpixel  uniformity  for  TV. 

1.2  OBJECTIVES 

The  principal  objective  of  the  work  described  in  this  report  was  to 
demonstrate  the  quantitative  separation  of  the  relative  contribution  of  the 
effects  of  atmospheric  turbulence,  laser  instability,  and  platform  motion,  by 
utilization  of  TV  imaging  systems.  The  basic  principles  of  this  separation 
have  been  established  by  previous  work  at  the  Naval  Postgraduate  School,  in 
terms  of  Fourier  transform  data  processing.  However,  the  previous  experimen¬ 
tal  work  there  has  been  carried  out  with  mechanical  slit-scanning  optics.  A 
principal  objective  was  thus  the  verification  that  television  imaging  systems 
can  yield  the  same  results  as  slit-scanning  systems.  The  techniques  were  to 
be  tested  in  the  laboratory,  but  be  capable  of  implementation  in  the  field. 

The  relevant  field  conditions  to  which  this  should  apply  are  to  include  the 
use  of  a  uniform  reflectivity  target  for  the  designator  spot.  However,  the 
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techniques  should  also  be  adaptable  to  the  case  of  a  realistic  nonuniform 
reflectivity  field  target,  such  as  a  tank.  Also,  although  not  part  of  the 
required  objective,  the  techniques  employed  were  to  be  chosen,  where  possible, 
to  be  extendable  later  to  the  use  of  moving  targets.  An  additional  objective 
of  the  work  was  to  verify  the  effectiveness  of  the  equipment  chosen  for  the 
experimental  program. 

1.3  SUMMARY  OF  PROCEDURES 

a.  Techniques  were  developed  and  tested  in  a  135  meter  long  medium- 
turbulence  optical  tunnel  for  the  purpose  of  obtaining  the  Modulation  Trans¬ 
form  Function  (MTF),  and  Cn2  (a  measure  of  the  severity  of  turbulence)  for 
the  path  by  means  of  a  television  imaging  system.  The  system  imaged  a  point 
laser  source  at  the  far  end  of  the  tunnel,  simulating  a  laser  mounted  in  the 
designator  target  during  field  tests.  The  television  signals  were  recorded  on 
analog  tape,  then  digitized  and  processed  later.  The  recording  techniques 
were  tested  and  found  to  be  sufficiently  linear.  A  computer  data  processing 
system  developed  for  the  TV  data  reduction  was  tested  on  the  observed  TV  image 
signals. 

b.  A  previously  proven  slit-scan  system  for  measuring  MTF  and  Cn 2 
and  the  TV  imaging  system  described  above  were  then  used  simultaneously  to 
view  identical  images  of  a  point  laser  source  at  a  distance.  The  images  were 
obtained  by  means  of  a  half-transmitting  mirror  splitter  in  the  telescope 
output.  The  signals  for  both  systems  were  recorded  and  data  processed  by 
Fourier  transform  techniques  to  yield  values  of  the  turbulence  structure 
constant,  Cn2,  as  a  quantitative  means  of  comparing  the  results  obtained 

by  the  two  systems.  The  results  are  summarized  below. 

c.  The  procedures  for  digital  subtraction  of  residual  background  (per¬ 
sisting  through  the  laser  line  filter),  in  TV  images  of  the  laser  spot  pattern 
on  a  uniform  screen,  were  tested,  using  a  Quantex  DS-30  digital  image  proces¬ 
sor.  The  data  processing  to  yield  the  MTF  from  the  resulting  background-free 
laser  designator  spot  was  then  tested. 

d.  Extension  of  the  techniques,  to  the  case  of  laser  designator  targets 
with  nonuniform  reflectivity,  was  analyzed  in  principle,  and  found  to  be  prac¬ 
tical.  Actual  performance  testing  of  this  process  awaits  the  availability  of 
a  larger  computer  planned  for  future  acquisition  at  Fort  Huachuca. 

1.4  SUMMARY  OF  RESULTS 

Measurements  of  the  OTF  of  the  atmosphere  made  by  means  of  TV  imaging 
systems  are  in  agreement  with  measurements  made  through  identically  the  same 
atmosphere  with  mechanical-line  scanning  equipment.  Also,  test  of  the 
recording  systems  needed  for  proposed  field  test  measurements  show  that  the 
recording  equipment  is  sufficiently  linear  to  permit  data  reduction  to  be 
carried  out  later.  The  general  techniques  are  capable  of  separating  the 
effects  of  atmospheric  turbulence,  internal  laser  designator  instability,  and 
effects  of  platform  motion. 
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1.5  ANALYSIS 


a.  The  Phase  I  analytical  effort  and  laboratory  experiments  were 
designed  to  establish  the  feasibility  of  applying  the  video  image  analysis 
technique  and  to  lay  the  groundwork  for  the  Phase  II  field  experiments.  The 
feasibility  assessment  required  the  evaluation  of  the  linearity  of  the  compo¬ 
nents  of  the  video  recording  system  and  also  required  the  development  of  a 
computerized  data  analysis  program  to  generate  the  Cn^  data  from  the  video 
data.  System  linearity  was  empirically  evaluated  by  comparing  the  output  of 
the  video  camera  with  the  final  output  from  the  disc  unit.  The  reproducibil¬ 
ity  achieved  demonstrated  the  basic  feasibility  of  replacing  the  mechanical 
slit  scan  system  with  the  video  system. 

b.  The  original  scope  of  the  Phase  I  study  did  not  include  the  setting 
up  of  a  complete  video  data  acquisition  and  analysis  instrumentation  configu¬ 
ration.  However,  the  availability  during  the  time  frame  of  the  study  of  a 
digital  video  analyzer,  video  disc  unit  and  a  desk  top  computer  made  it 
possible  to  set  up  a  complete,  self  contained  Cn^  measurement  system.  The 
measurement  system  utilized  the  same  components  used  by  most  test  ranges  to 
record  and  analyze  laser  designator  spot  characteristics.  The  availability  of 
the  instrumentation  configuration  accelerated  the  generation  of  the  supporting 
software  program  to  derive  the  C data  from  the  spot  imagery  and  the 
results  of  that  work  are  covered  in  detail  in  Appendix  C.  The  results  show 
that  a  comparison  of  the  measured  versus  predicted  patterns  will  establish  the 
level  of  jitter  attributable  to  the  laser  platform. 

c.  Improvements  to  the  data  processing  capability  can  be  made  by  con¬ 
trolling  the  video  disc  and  the  tape  unit  through  an  interface  bus  instead  of 
manually  as  done  for  this  analysis  program.  However,  most  of  the  currently 
operational  designator  spot  analysis  systems  use  such  a  bus  and  the  software 
in  Appendix  C  can  be  adapted  to  take  advantage  of  that  capability.  Addition¬ 
ally,  a  significant  increase  in  the  speed  of  the  data  analysis  operation  can 
be  achieved  if  a  computer  system  with  more  speed  and  greater  memory  than  the 
HP-9825  is  used.  It  is  currently  planned  to  use  an  HP-1000  for  the  field  test 
phase  of  the  investigation. 

d.  The  use  of  the  beam  splitter  technique  to  compare  the  performance  of 
the  video  system  and  the  NPS  mechanical  slit  scanning  system  represented  an 
ideal  solution  to  the  path  control  problem.  The  commonality  of  both  the 
spatial  and  temporal  inputs  to  the  two  systems  greatly  simplified  the  analysis 
of  the  data.  The  data  obtained  during  the  comparison  tests  is  covered  in 
detail  in  Appendix  D  and  it  established  that  the  video  system  did  not  produce 
any  degradation  in  the  signal  data. 

e.  Due  to  a  prior  commitment  for  the  video  disc  and  digitizer  equip¬ 
ments,  these  items  were  not  available  to  use  in  the  analysis  of  the  comparison 
test  data.  As  a  result,  an  alternate  data  analysis  method  was  used  and  a 
considerably  modified  computer  program  prepared.  The  alternate  method,  as 
covered  in  detail  in  Appendix  D,  demonstrated  that  a  storage  oscilloscope  can 
be  used  for  analyzing  the  data  and  that  a  good  approximation  of  the  point 
spread  function  can  be  made  by  recording  a  single  TV  line  through  the  laser 
spot  and  using  it  to  calculate  the  one-dimensional  line  spread  function.  The 
computer  program  generated  for  this  data  analysis  is  presented  in  algorithm 
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form  in  order  to  facilitate  its  transfer  to  other  computer  systems.  The  use 
of  the  digital  storage  oscilloscope  demonstrated  that,  if  data  reduction 
turnaround  time  is  not  critical,  the  video  data  analysis  can  be  performed 
using  a  $6,000  instrument  in  place  of  a  $60,000  video  digitizer  and  disc 
system.  For  path  weighted  measurements  not  involving  the  high  data 
rate  laser  designator  spot  profile  analyses,  the  digital  storage  oscilloscope 
system  represents  a  modest  investment  for  the  capability  provided. 

f.  The  conduct  of  the  laboratory  type  system  comparison  tests  repre¬ 
sented  an  interim  step  between  the  originally  planned  Phase  I  study  effort  and 
the  Phase  II  field  tests.  The  successful  completion  of  the  interim  step 
laboratory  measurements  has  increased  the  probability  for  the  successful 
conduct  of  the  field  tests  and  has  also  eliminated  the  requirement  for  a  major 
portion  of  the  preliminary  technique  verification  tests  originally  planned  for 
the  Phase  II  field  test  program.  The  proposed  field  tests  are  covered  in 
Appendix  A. 

g.  By  generating  the  analysis  programs  for  the  experimental  tests  for 
two  levels  of  instrumentation  support,  the  end  products  of  Phase  I  study  and 
experimental  effort  can  be  applied  to  a  range  of  both  laboratory  and  field  0TF 
and  Cn^  measurement  requirements.  This  includes  applying  the  techniques 

to  standard  testing  of  E-0  systems  and  devices  in  the  field  as  well  as  for 
characterizing  atmospheric  path  conditions  for  degraded  environment  testing. 

1.6  CONCLUSIONS 

The  proposed  field  experiments  to  separate  the  effects  on  laser  desig¬ 
nators  of  atmospheric  turbulence,  internal  laser  designator  instability,  and 
platform  motion  can  be  successfully  carried  out  with  the  techniques  and  equip¬ 
ment  recommended. 

1.7  RECOMMENDATIONS 

It  is  recommended  that  the  proposed  field  test  experiments  be  carried  out 
to  separate  the  effects  of  turbulence,  internal  laser  instability,  and  plat¬ 
form  motion  in  laser  designators. 
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SECTION  2.  DETAILS  OF  INVESTIGATION 


2.1  INTRODUCTION 

a.  The  experimental  and  analytical  program  carried  out  has  been  that 
necessary  to  establish  the  validity  of  experiments  proposed  for  a  later  field 
test  program.  That  test  program  is  summarized  in  appendix  A.  Briefly,  the 
proposed  field  test  program  addresses  the  problem  of  separating  the  types  of 
instability  in  laser  designators,  i.e.  ,  turbulence,  the  internal  stability  of 
the  laser  designator  and  its  optics,  and  the  platform  instability. 

b.  The  effects  of  atmospheric  turbulence  are  broadening  of  the  image 
profile,  wander  of  the  center  of  the  beam  and  scintillation  within  the  beam 
profile.  The  scintillation  is  the  most  difficult  to  handle  in  detail.  How¬ 
ever,  it  can  be  circumvented  by  dealing  with  the  profile  as  observed  in 

the  average  of  a  number  of  individual  pulse  profiles. 

c.  The  effects  of  turbulence  on  beam  propagation  have  been  investigated 
at  a  number  of  laboratories,  with  perhaps  the  longest  and  roost  complete  sus¬ 
tained  program  being  at  the  Naval  Postgraduate  School  (ref  1,  app  E).  This 
work  has  led  to  a  well-verified  analytical  model  for  the  beam  profile  to  be 
expected  after  a  beam  has  traversed  an  optical  path  through  the  atmosphere. 
This  is  expressed  in  terms  of  the  OTF  or  Optical  Transfer  Function  of  the 
atmosphere.  Two  forms  of  this  are  equally  useful.  One  expresses  the  OTF  for 
the  average  of  a  number  of  beam  profiles  without  centering  (or  tracking)  the 
beam.  The  other  applies  to  a  beam  that  is  centered  (or  tracked)  before  aver¬ 
aging  a  number  of  profile  measurements.  In  the  course  of  centering  before 
averaging,  a  tally  of  the  offsets  required  gives  the  variance  of  the  beam 
"wander."  Both  these  quantities  are  functions  of  the  turbulence  constant  for 
optical  index,  Cn^»  and  the  range,  as  well  as  the  wavelength. 

d.  The  beam  spread  resulting  from  a  number  of  causes,  such  as  the  turbu¬ 
lence,  designator  instability,  and  platform  instability,  is  the  convolution  of 
the  profiles  from  each  cause.  This  allows  the  Fourier  transform  convolution 
theorem  to  be  utilized.  This  says  that  the  transform  of  the  convolution  of 
several  functions  is  the  product  of  the  Fourier  transforms  of  those  functions. 
This  can  be  reversed  and  the  functions  unfolded  by  dividing  the  transforms 
point-by-point.  This  technique  has  been  used  to  separate  the  types  of 
broadening. 

e.  For  the  case  of  wander  where  the  statistical  behavior  is  Gaussian, 
the  variance  of  the  displacement  is  the  sum  of  the  variances  of  the  two  com¬ 
ponent  wanders.  This  allows  separation  of  the  wander  terms. 

2.2  THEORETICAL  BACKGROUND 

a.  Separation  of  the  atmospheric  contribution  to  the  beam  wander  and 
broadening  involves  measuring  the  turbulence  structure  constant  for  optical 
index,  Cn2.  This  can  be  determined  in  a  number  of  ways.  For  example,  it 
can  be  determined  by  measuring  the  temperature  at  points  along  the  optical 
path.  Although,  in  principle,  it  could  be  measured  with  thermal  sensors,  this 
has  not  been  very  practical  in  most  instances.  One  difficulty  is  that  the 
optical  properties  of  interest  here,  the  broadening  and  wander,  depend  on  the 
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integrated  value  of  Cn^  along  the  optical  path.  Also  the  relative  impor¬ 
tance  of  Cn^  varies  with  the  position  along  the  path,  and  this  weighting 
factor  depends  on  the  particular  type  of  optical  property  involved.  The  rela¬ 
tive  path-position  weighting  of  Cn^  is  shown  in  figure  1,  for  several 
different  optical  situations.  The  situation  for  an  optical  beam  formed  by  a 
projector  lens  system  is  shown  in  curve  A,  for  the  case  where  the  lens  system 
is  on  the  right  and  the  target  on  the  left.  By  a  reciprocity  theorem,  a 
telescope  system  behaves  in  the  same  way  as  the  projector  if  the  telescope  is 
on  the  right  and  a  point  source  is  on  the  left. 


b.  Other  methods  of  measuring  Cn^ ,  for  example,  by  scintillation, 
produce  a  path-position  weighting  distribution  as  shown  in  curve  B.  This 
emphasizes  the  center  of  the  path  and  does  not  give  the  proper  distribution  to 
predict  the  behavior  of  a  laser  designator.  In  this  case  then,  it  is  impor¬ 
tant  to  use  a  telescope  to  measure  Cn2  to  obtain  the  weighting  needed  to 
predict  the  behavior  of  the  designator. 


2 . 3  EXPERIMENTAL  PROGRAM 


The  experimental  program  to  verify  the  practicality  of  the  proposed  field 
test  program  involved  several  facets.  The  general  objective  was  verification 
of  the  general  applicability  of  the  use  of  TV  imagers  in  place  of  mechanical 
line  slit-scanners  for  the  purpose  of  separating  the  types  of  spot  instabil¬ 
ity.  This  also  involved  the  verification  of  the  suitability  and  linearity  of 
individual  components  of  the  proposed  overall  test  system.  We  will  look  first 
at  the  tests  of  system  components,  followed  by  an  experimental  comparison  of 
TV  and  mechanical  line-scanner  techniques. 


a.  System  Components 


(1)  Silicon  Vidicon.  A  silicon  focal  plane  screen  is  needed  for  the 
vidicon  viewing  the  designator  spot  pattern  in  order  to  reach  the  1.06  micro¬ 
meter,  near-IR,  wavelength.  The  usual  TV  vidicons,  or  other  TV  cameras  such 
as  orthicons,  are  optimized  for  the  visual  range  of  wavelengths  and  do  not 
reach  this  wavelength.  A  TV  technique  is  needed  because  the  pulse  length  of 
the  designator  is  too  short  (a  few  nanoseconds)  for  a  split-scanning  system 
to  function.  The  image  of  the  designator  spot  is  stored  on  the  semiconducting 
sensitive  surface  of  the  vidicon.  The  electron  beam  raster  scan  interrogates 
the  Image  to  produce  the  usual  composite  TV  video  format.  It  was  found  on 
trial  that  each  of  the  two  interlaced  "fields"  of  the  raster  completely  reset 
the  sensitive  surface.  The  storage  time  is  thus  the  1/60  second  between  field 
scans.  A  test  of  the  performance  of  the  vidicon  was  carried  out  in  the  labor¬ 
atory  by  imaging  a  broad  spot  produced  by  GaAs  laser  at  a  wavelength  of  0.905 
micrometers . 


(2)  Vidicon  for  Atmospheric  OTF  Measurement.  A  second  vidicon  is 
used  to  view  a  point  laser  source  located  in  the  target  screen,  in  order  to 
provide  the  data  concerning  the  atmospheric  OTF.  (See  fig.  1,  app  A.)  This 
system  is  less  demanding  as  far  as  pulse  length  is  concerned,  and  the  wave¬ 
length  can  be  any  value  available  in  a  convenient  small  laser,  as  long  as  it 
differs  from  the  1.06  micrometers  of  the  designator.  In  practice  a  HeNe 
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Figure  I.  Relative  Weighting  of  Cn2  as  a  Function  of  Position  Along  the  Path,  fo 
and  for  Scintillation.  The  Telescope  End  of  the  Path  is  at  the  Right. 


laser  was  used,  emitting  .6328  micrometers.  A  power  of  1-3  milliwatts  is 
sufficient.  The  viewing  vidicon  has  a  narrow-band  interference  filter  in  its 
optical  train  to  reduce  the  background  signal  and  exclude  the  designator  wave¬ 
length.  In  order  to  have  sufficient  magnification  to  resolve  the  spread  of 
the  image  due  to  the  atmospheric  turbulence,  this  vidicon  viewed  the  source 
through  a  6-inch  diameter  Cassegrain  telescope,  with  a  90-inch  focal  length. 
This  proved  to  be  too  low  a  magnification  so  a  "Barlow"  negative  lens  was 
inserted  just  ahead  of  the  vidicon  to  increase  the  magnification  by  a  factor 
of  3.  This  caused  the  typical  image  to  occupy  a  <~<dth  of  the  order  of  10 
pixels,  sufficient  to  analyze  the  image. 

(3)  Video  Tape  Recording 

(a)  The  vidicon  signal  is  continuously  recorded  on  a  Panasonic, 
U-Vision,  type  NV9240  TV  tape  recorder  system.  One  matter  of  immediate  con¬ 
cern  in  whether  the  tape  recorder  is  sufficiently  linear  to  be  used  for  data 
recording.  Not  only  must  a  linear  analog  recorder  be  used  in  the  field,  but 
it  is  necessary  to  later  identify  and  play  back  with  precise  identification, 
several  Individual  fields.  This  requires  a  disc  recorder  with  field  and  frame 
labeling.  For  this  purpose  a  Eigen,  model  16-10,  video  disc  recorder  was 
used. 


(b)  To  test  the  linearity  of  the  complete  chain  of  amplifiers 
and  two  recording  systems,  images  of  a  point  source  broadened  by  the  atmos¬ 
phere  were  compared.  Two  types  of  signals  are  illustrated  in  figure  2. 

Figure  2a  shows  a  single  TV  line  that  passes  through  the  center  of  the  re¬ 
corded  spot.  The  upper  trace  is  direct  from  the  TV.  The  lower  trace  is  after 
the  complete  amplifier  chain  and  two  recording  systems.  The  patterns  are 
sufficiently  alike  to  give  comparable  shape  parameter  measurements.  In  figure 
2b,  a  number  of  successive  TV  lines  passing  through  the  image  are  displayed. 
Again,  the  top  figure  is  direct  and  the  lower  figure  is  after  the  complete 
amplifier  chain.  The  reproducibility  is  excellent. 

(4)  Background  Subtraction.  Although  both  vidicon  systems  have 
narrow-band  filters  to  reduce  background  light,  it  is  not  completely  elimi¬ 
nated.  In  order  to  remove  this,  the  signals  are  digitized  in  a  Quantex  DS-30 
TV  frame  digitizer.  This  system  can  subtract  the  signals  of  one  field  from 
that  of  another.  This  is  carried  out,  using  the  field  recorded  signals,  to 
obtain  the  image  signals  of  the  designator  and  OTF  vidicons  alone. 

(5)  Data  Processing  for  Atmospheric  Cn^ 

(a)  In  order  to  carry  out  a  calculation  of  the  atmospheric 
Cn^  from  an  image  of  a  point  source  at  a  distance,  the  OTF  of  the  atmos¬ 
phere  is  measured.  A  complete  curve  of  the  OTF  is  also  needed  for  the  separa¬ 
tion  of  the  various  causes  of  spot  broadening  and  wander.  Complete  details  of 
the  video  data  processing  techniques  are  given  in  appendixes  C  and  D. 

(b)  To  measure  the  OTF,  a  single  TV  image  field  is  first  trans¬ 
ferred  to  the  Quantex  DS-30  from  the  video  disc.  A  background  field  is  then 
transferred  from  the  video  disc.  This  is  subtracted  from  the  image  field  in 
the  DS-30.  The  background  is  ordinarily  very  small  as  the  image  has  been 
taken  through  a  narrow-band  interference  filter  matched  to  the  laser  source 
wavelength. 
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Fig.  la.  Signals  for  the  image  of  a  point  sou: 
for  the  TV  line  passing  through  the  image  sen: 
Too  trace :  Direct  signal  from  Yidicon 
Bottom  trace:  Signal  after  recording  in  tape  : 
di s c  recorder ,  DS-50  digitizer  and  recorder,  , 
to  analong  conversion  for  viewing. 


HI, 


Fig.  1  b .  Signals  for  the  image  of  a  point  source 
tor  a  sequence  of  TV  lines  covering  the  image.  T 
envelope  of  these  peaks  is  processed  to  give  ina 
profile. 

Top  trace:  Direct  signal  from  Yidicon 
Bottom  trace :  Signal  after  recording  in  tape  rec 
disc  recorder ,  DS-50  digitizer  and  recorder,  dig 
to  analog  conversion  for  viewing. 


(c)  After  subtraction  in  the  DS-30  of  a  complete  background 
field  from  a  complete  image  field  that  includes  an  image  of  the  point  laser 
source,  a  single  raster  line  at  a  time  is  transferred  to  the  HP-9825  calcu¬ 
lator.  The  line  transferred  is  limited  in  length  to  about  twice  the  number 
of  pixels  required  to  express  the  image  shape  -  usually  20  to  50  pixels.  The 
line  signal  is  then  summed  in  the  HP-9825.  This  sum  corresponds  to  the  integ¬ 
rated  light  passing  through  the  linear  slit  of  a  slit-scanner.  Successive 
lines  are  then  summed  similarly  and  the  sequence  of  sums  becomes  the  "line- 
spread  function."  This  corresponds  to  the  signal  function  obtained  by  a 
slitscanner.  This  line-spread  function  is  stored  in  the  HP-9825  memory.  The 
next  field  is  then  analyzed  similarly.  To  obtain  the  "long-term"  line-spread 
function,  successive  line-spread  functions  are  averaged  in  the  HP-9825. 

(d)  To  obtain  the  "short-term”,  or  "image-centered  line-spread 
function,"  the  center  of  area  of  each  line-spread  function  is  calculated  and 
stored  in  memory.  The  line  center  is  then  shifted  to  a  predetermined  position 
and  the  average  of  a  number  of  such  centered  functions  is  the  image-centered 
line-spread  function.  The  RMS  average  of  the  center  displacements  needed  to 
center  is  the  RMS  "image-wander"  value.  Its  square  is  the  variance  of  the 
wander. 


(e)  Next,  the  Fourier  transforms  of  the  two  line-spread  func¬ 
tions  are  carried  out  in  the  HP-9825.  Each  of  these  is  then  divided  point- 
by-point  by  the  transform  of  the  line-spread  function  of  the  imaging  telescope 
in  the  absence  of  turbulence.  This  removes  the  effects  of  the  instrument. 

The  result  is  two  curves  of  the  OTF  of  the  atmosphere,  the  line-term,  and  the 
image-centered,  OTF  values.  The  curves  can  be  plotted  out  when  desired. 

Next,  the  OTF  curves  are  fitted  to  the  Fried  model  to  yield  Cn 2.  Ordin¬ 
arily,  the  two  values  of  Cn2  are  identical  if  the  image  motion  is  due  to 
turbulence  alone,  and  the  turbulence  is  not  unusual.  The  RMS  wander  displace¬ 
ment  also  yields  a  value  of  Cn2 . 

(f)  Determination  of  the  values  of  Cn2  is  not  actually 
needed  for  much  of  the  work  to  be  carried  out  here,  but  that  value  is  a  con¬ 
venient  single  number  that  characterizes  the  level  of  atmospheric  turbulence 
along  the  optical  path.  The  actual  shape  of  the  line-spread  functions  will  be 
what  is  used  in  the  data  processing  to  separate  the  various  types  of  line- 
spreads.  For  the  wander,  the  actual  values  of  the  variance  are  also  the 
values  that  will  be  used  in  the  separation  of  various  wander  causes. 

(6)  Data  Processing  for  the  Designator  Spot  Pattern 

(a)  The  details  of  this  data  reduction  are  provided  in  appen¬ 
dixes  C  and  D.  Before  discussion  of  the  data  analysis,  a  few  features  of  the 
equipment  operation  need  to  be  mentioned.  It  would  be  desirable  for  the  tim¬ 
ing  of  the  TV  field  initiation  to  be  synchronized  with  the  designator  flash 
occurrence,  in  order  to  have  the  flash  occur  during  the  TV  blanking  and  fly¬ 
back.  If  the  designator  flash  occurs  randomly  relative  to  the  vidicon,  it 
could  occur  during  the  period  in  which  the  electron  beam  is  scanning  the  area 
of  the  image.  This  would  make  the  data  processing  much  less  dependable.  In 
practice,  an  occasional  TV  field  would  be  lost.  Synchronization  would  make 
every  flash  useful. 
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(b)  The  image  of  the  designator  spot  taken  by  the  close-up 
vidicon  is  next  analyzed.  Basically,  the  same  data  processing  procedure  is 
used  as  to  obtain  the  OTF  of  the  atmosphere.  The  TV  image  is  digitized  and 
the  background  subtracted  in  the  DS-30.  The  average  of  a  number  of  images  is 
obtained;  both  the  image-centered  average,  and  the  uncentered  average.  The 
process  of  obtaining  the  image-centered  average  also  yields  the  variance  of 
the  image  wander. 

(c)  The  image-centered  line-spread  function  is  the  convolution 
of  the  actual  spread  of  the  designator  and  the  spread  due  to  the  atmosphere. 

To  separate  the  two  spreads,  the  Fourier  transform  is  taken  off  the  image- 
centered  line-spread  function.  This  is  then  divided,  point-by-point,  by  the 
Fourier  transform  of  the  atmospheric  line-spread  function  (the  OTF  of  the 
atmosphere).  The  resulting  function  is  then  reinverted  to  obtain  the  spread 
due  to  the  designator  alone.  This  is  plotted  out  on  the  plotter.  This  curve 
is  one  of  the  primary  goals  of  the  measurement.  This  can  be  used,  then,  as 
the  information  needed  to  correct  errors  in  the  optics  of  the  laser  designator 
system. 


(d)  The  variance  of  the  image  wander,  obtained  above,  is  the 
result  of  both  the  actual  wander  of  the  designator  beam  and  the  atmospheric 
wander.  The  two  can  be  easily  unfolded,  assuming  that  the  statistics  are 
reasonably  Gaussian.  The  atmospheric  wander  of  the  designator,  due  to  inter¬ 
nal  instability,  would  also  be  Gaussian.  The  variance  of  the  combined  wander 
is  then  simply  the  sum  of  the  wanders  of  the  atmosphere  and  the  designator. 
Hence,  the  atmospheric  wander  variance  is  simply  subtracted  from  the  observed 
total  wander  variance  to  obtain  the  wander  of  the  designator.  This  also  is 
one  of  the  primary  goals  of  the  measurement  as  this  information  can  help  in 
the  readjustment  of  the  laser  and  optics  to  reduce  the  wander. 

(e)  Analysis  of  the  pattern  projected  on  a  uniform  target  has 
been  demonstrated  in  work  at  the  Naval  Postgraduate  School,  using  a  pulsed 
GaAs  laser  source.  The  images  have  been  digitized,  subtracted  from  back¬ 
ground,  and  the  spot  profile  processed  to  yield  a  line-spread  function,  as 
well  as  a  spot  wander  value.  These  have  then  been  separated  from  the 
atmospheric  spread  and  wander. 

b.  TV  versus  Line-Scanner  Comparison 

(1)  Direct  comparison  of  OTF  values  from  the  NPS  mechanical  line- 
scanner  with  those  from  the  TV  system  was  accomplished  by  use  of  a  beam  split¬ 
ter  in  the  output  optics  of  a  270-inch  focal  length  Cassegrain  telescope.  The 
two  images  analyzed  were  thus  identical.  The  telescope  viewed  a  point  (HeNe) 
laser  source  at  the  far  end  of  a  tunnel,  132  meters  in  length.  Turbulence  in 
the  tunnel  was  produced  by  overhead  heat  duct  ports  and  amounted  to  moderate 
turbulence. 

(2)  The  video  signals  from  the  TV  system  were  recorded  with  the 
system  proposed  here  for  field  use.  The  recorded  signals  were  analyzed  later 
by  the  techniques  previously  described  in  this  report,  with  the  output  in  the 
form  of  a  plotted  curve  of  OTF  as  a  function  of  angular  spatial  frequency. 
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(3)  The  NPS  mechanical  line-scanner  signals  were  recorded  on  a 
Precision  Data  frequency  modulated  tape  recorder,  with  a  recording  bandwidth 
of  100  kHz.  This  recorder  is  frequently  used  for  this  purpose  with  the  NPS 
line-scanner,  although  recently  the  online  data  processor  has  made  this  re¬ 
corder  unnecessary.  In  this  case,  it  made  careful  processing  possible  for 
identically  the  same  time  interval  as  that  used  for  the  TV  data. 

(4)  A  pair  of  curves  of  OTF  are  shown  in  figures  3  and  4  for  the  NPS 
line  scanner  and  the  TV  system,  respectively.  The  data  is  for  the  long-term 
average.  The  output  of  the  NPS  line-scanner  yields  the  actual  data  points  as 
well  as  the  best-fit  curve  from  the  Fried  model  (ref  4,  app  E).  The  output  of 
the  TV  system  displayed  only  the  best-fit  Fried  curve,  at  the  time  of  that 
data  reduction.  In  the  future,  it  will  also  display  the  actual  data  points. 
The  two  OTF  curves  can  differ  only  in  their  horizontal  scale  because  they  are 
curves  for  the  same  model.  A  customary  measure  for  comparison  is  the  angular 
spatial  frequency  at  which  the  OTF  has  decreased  to  a  fraction,  1/e,  of  its 
initial  value  of  unity.  This  point  is  often  used  as  a  measure  of  the 
coherence  "length."  For  the  NPS  scanner,  the  1/e  point  occurs  at  26.7  lines 
per  railliradian.  For  the  TV  system,  the  1/e  point  occurs  at  29.6  lines  per 
milliradian. 

(5)  The  above  values  are  in  agreement  within  approximately  10  per¬ 
cent.  This  corresponds  to  an  agreement  within  10  percent  for  Cn,  or  within 
20  percent  for  Cn^.  This  degree  of  agreement  indicates  that  the  TV  system 
is  adequate  to  replace  the  mechanical  line-scanner  for  measurements  of  the 
OTF.  That  agreement  is  quite  good,  and  since  much  of  the  use  of  the  two  TV 
systems  for  OTF  involves  use  of  the  results  of  one  system  in  terms  of  the 
other,  the  matched  pair  of  two  TV  systems  should  reduce  the  effects  of  system¬ 
atic  error. 
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Figure  3.  OTF  of  the  Atmosphere  for  the  NPS  Line-Scanner 
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OTF  of  the  Atmosphere  for  the  TV  System 
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SECTION  3.  APPENDIXES 
APPENDIX  A.  PROPOSED  FIELD  TEST  PROGRAM 

1 .  INTRODUCTION 

The  proposed  test  program  has  three  general  phases.  In  the  first  phase, 
a  uniform  reflectivity  target  screen  would  be  used.  These  techniques  would  be 
adapted  in  the  second  phase  to  permit  use  of  a  real  target  such  as  a  tank. 

The  techniques  utilized  in  the  first  two  phases  would  be  chosen  so  that  they 
are  adaptable  to  a  third  possible  phase  in  which  real  targets  would  be  in  mo¬ 
tion. 

2.  PHASE  ONE,  UNIFORM  TARGETS 

In  the  first  phase  of  the  program,  field  measurements  would  be  made  with 
a  target  screen  with  uniform  reflectivity  at  a  representative  range.  The 
equipment  is  planned  to  be  arranged  as  shown  in  figure  1.  The  measurements 
and  data  reduction  would  be  carried  out  in  three  successive  steps. 

a.  Step  One 

(1)  In  the  first  step,  field  measurements  would  be  made  with  a 
strapped-down  laser  designator  with  its  spot  trained  onto  a  uniform  target 
screen.  These  measurements  would  serve  to  determine  the  stability  of  the  de¬ 
signator  laser  and  optical  system  by  separating  the  atmospheric  effects  from 
the  overall  effects. 

(2)  The  target  screen  would  be  viewed  with  a  silicon  vidicon  located 
a  short  distance  from  the  screen  along  the  line  between  the  screen  and  the  de¬ 
signator,  but  slightly  off  that  line  to  avoid  obscuring  the  screen  from  the 
designator  location  as  shown  in  figure  1. 

(3)  A  silicon  focal  plane  screen  is  needed  for  this  vidicon  in  order 
to  reach  the  1.06  micrometer,  near-IR,  wavelength.  The  usual  TV  vidicons,  or 
other  TV  cameras  such  as  orthicons,  are  optimized  for  the  visual  range  of 
wavelengths  and  do  not  reach  this  wavelength. 

(4)  To  reduce  the  background  signal,  the  close-up  vidicon  would  be 
equipped  with  a  narrow-bandpass  optical  filter  matched  to  the  designator 
wavelength.  The  vidicon  would  view,  alternately,  the  designator  spot  plus 
background,  and  the  background  alone.  The  video  signals  would  be  tape  re¬ 
corded  during  tests.  Later,  the  signals  for  the  background  image  would  be 
subtracted  from  the  signals  for  the  background  plus  designator  spot  to  give 
the  signals  for  the  spot  alone.  This  subtraction  would  be  carried  out  in  a 
Quantex  DS-30  video  signal  processor. 

(5)  In  addition  to  the  close-up  vidicon,  another  vidicon,  with  long 
focal  length  (telescope  optics),  would  be  located  near  the  designator  trans¬ 
mitter  and  would  view  a  point  source  laser  located  in  the  target  screen  as 
shown  in  figure  1.  That  laser  would  operate  at  a  different  wavelength  from 
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the  designator,  e.g.,  HeNe,  at  6028  micrometers.  The  telescope  would  have  a 
narrow  band  filter  to  match  this  laser  line.  This  will  eliminate  the  desig¬ 
nator  signals,  but  not  completely  the  daylight  continuum  background.  Separate 
recording  of  signal  and  background  will  permit  subtraction  of  background  here 
in  the  same  manner  as  for  the  close-up  vidicon.  This  vidicon  can  be  a  more 
common  visual  range  vidicon  than  the  silicon  vidicon  that  is  required  for  the 
1.06  micrometer  designator  signals. 

(6)  Processing  of  the  taped  data  from  the  telescope  vidicon  would 
give  a  quantitative  determination  of  the  effects  of  the  atmospheric  turbu¬ 
lence.  This  data  would  then  be  used,  together  with  the  taped  data  from  the 
close-up  vidicon  that  views  the  designator  spot  on  the  target  screen,  to  un¬ 
fold  the  atmosphere  effects  from  the  designator  spot  width  and  wander,  using 
Fourier  transform  techniques.  The  diffraction  spread  and  any  intentional  geo¬ 
metrical  spread  of  the  designator  spot  would  also  be  removed.  The  remaining 
spread  (beam  profile)  and  wander  would  then  be  entirely  that  due  to  improper 
laser  adjustment,  or  instability  of  the  laser  and  associated  optics.  This 
data  should  then  make  it  possible  to  correct  or  minimize  the  defects  in  the 
designator  system.  Even  if  not  corrected,  the  system  would  then  be  charac¬ 
terized,  and  is  ready  for  use  in  the  next  step. 

b.  Step  Two 

In  the  second  step,  the  laser  designator  would  be  placed  on  the  un¬ 
stable  platform  (mini  RPV).  The  mini  RPV  would  be  mounted  in  as  realistic  a 
manner  as  possible,  with  engine  and  all  systems  running  to  simulate  real  oper¬ 
ational  conditions,  particularly  with  respect  to  vibration,  but  with  the  mini 
RPV  still  on  a  ground-based  platform.  This  is  necessary  so  that  the  telescope 
vidicon  can  be  mounted  next  to  it,  to  view  the  point  source  in  the  target 
screen  and  thereby  obtain  a  measurement  of  the  effects  of  atmospheric  turbu¬ 
lence.  This  would  evaluate  the  actual  beam  pattern  of  the  designator,  as  it 
would  be  in  the  air,  by  removing  the  atmospheric  effects.  The  beam  pattern 
would  now  be  characterized  sufficiently  that  it  can  be  used  in  the  third  step. 

c.  Step  Three 

(1)  In  the  third  step,  the  mini  RPV  would  be  airborne,  directing  its 
laser  designator  to  the  uniform  target  screen.  The  close-up  silicon  vidicon 
would  view  the  screen  as  before.  In  this  step,  the  distant  telescope  vidicon 
viewing  the  point  source  is  not  necessary,  nor  could  it  contribute  to  the 
information,  if  used.  It  would,  of  necessity,  be  ground-based  and  would  not 
look  through  the  same  atmospheric  path  as  the  beam  from  the  mini  RPV  to  the 
target.  The  beam  from  the  designator  on  the  mini  RPV  is  now  well  enough  known 
that  it  can  be  used  to  evaluate  the  atmospheric  turbulence  at  the  instant  of 
designator  firing.  The  data  tape  from  the  close-up  vidicon  would  yield  a 
sequence  of  beam  profiles.  The  center  of  each  of  these  is  calculated  and  the 
profile  centered  before  averaging.  The  averaged  profile;  so  obtained,  is  then 
Fourier  processed  to  remove  the  known  designator  profile;  that  has  been  deter¬ 
mined  in  step  two  above.  From  this  reduced  profile,  the  atmospheric  turbu¬ 
lence  is  determined.  Knowing  this  turbulence,  the  expected  wander  of  the  spot 
center  is  calculated.  This  is  then  unfolded  from  the  directly-measured  spot 
wander  to  give  the  motion  of  the  spot  due  to  the  platform  motion  alone. 
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(2)  All  of  the  above  general  processes  have  been  proven  by  applica¬ 
tion  in  a  slightly  different  situation,  namely  determining  the  optical  proper¬ 
ties  of  lasers  transmitted  from  shipboard.  A  new  feature  in  the  present 
experiments  is  the  use  of  vidicons  for  imaging.  Scanning-slit  telescopes  have 
been  used  previously  to  carry  out  the  same  type  of  measurements.  In  the  re¬ 
sults  reported  in  this  report,  data  reduction  for  imaging  with  the  vidicons 
has  now  been  tested  in  laboratory  work  at  the  Naval  Postgraduate  School. 

3.  PHASE  TWO,  NONUNIFORM  TARGETS 

a.  In  the  second  phase  of  the  program,  field  measurements,  similar  to 
phase  one  would  be  carried  out,  but  with  nonuniform  reflectivity  targets,  such 
as  a  tank  with  natural  landscape  as  background. 

b.  The  measurements  would  differ  from  those  of  phase  one  in  that  an  ad¬ 
ditional  laser,  with  the  same  wavelength  as  the  designator,  would  be  used  to 
provide  pulse  illumination  of  the  target  from  a  close-up  location.  This 
pulsed  illumination  would  alternate  with  the  designator  pulses.  The  close-up 
vidicon  would  view  the  scene  under  pulse  illumination  plus  background  light, 
and  then  under  background  light  alone.  Both  signals  would  be  tape  recorded. 
Later,  the  i_«o  images  would  be  subtracted  to  yield  the  signals  for  illumina¬ 
tion  by  laser  light  alone.  The  close-up  vidicon  would  then  view  the  scene 
when  the  designator  spot  is  pulsed.  This,  again,  would  be  tape  recorded.  The 
difference  of  this  scene  and  the  background  would  be  used  later  for  the  sig¬ 
nals  for  the  designator  alone.  Later,  in  data  reduction,  the  designator  spot 
signals  would  be  divided  pixel-by-pixel  by  the  signals  for  uniform  laser  illu¬ 
mination.  This  process  normalizes  the  spot  profile  to  be  the  same  as  if  it 
had  reflected  from  a  uniform  reflectivity  target. 

c.  Steps  one,  two  and  three  above  would  then  be  carried  out,  after  the 
normalization,  to  analyze  the  behavior  of  the  laser  and  platform. 

4.  PHASE  THREE,  FUTURE  MOVING  TARGETS 

The  systems  for  use  in  phases  one  and  two  would  be  set  up  so  as  to  be  ex¬ 
tendable  to  laser  designator  analysis  for  moving  targets.  This  will  involve 
all  the  techniques  of  phases  one  and  two,  including  steps  one,  two  and  three, 
plus  provision  of  background  laser  illumination  by  flashes  that  follow,  or 
lead,  the  designator  pulses  as  closely  as  possible.  In  practice,  this  means 
following  or  leading  by  one  TV  field  time  interval.  This  is  1/60  second  for 
standard  American  TV.  Use  of  a  more  rapid  scan  TV  system  would  improve  the 
situation.  However,  the  normalization  operation  would  be  improved  by  a 
lateral  shift  of  the  uniform  illumination  calibration  image,  to  more  exactly 
coincide  with  the  image  at  the  time  of  the  designator  flash.  The  proper  dis¬ 
placement  could  be  determined  by  carrying  out  a  correlation  of  the  two  images 
as  a  function  of  offset,  with  the  final  offset  determined  bv  the  maximum  cor¬ 
relation  between  the  images.  This  would  be  a  time-consuming  operation  if  both 
lateral  and  vertical  offset  were  involved.  If  the  direction  of  offset  is 
known,  for  example,  for  a  vehicle  that  is  constrained  to  move  on  a  road  sur¬ 
face  of  known  inclination,  then  the  correlation  operation  would  be  greatly  re¬ 
duced.  Carrying  out  this  type  of  measurement  seems  well  within  the  capability 
of  the  equipment  to  be  provided  for  phases  one  and  two,  except  that  a  larger 
computer  will  probably  be  needed.  This  should  be  borne  in  mind  in  selecting 
the  computer  for  use  in  phases  one  and  two. 
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APPENDIX  B.  METHODOLOGY  INVESTIGATION  PROPOSAL 


1.  TITLE .  Optical  Turbulence  Measurement 

2.  CATEGORY .  Electro-Optics 

3.  INSTALLATION.  US  Army  Electronic  Proving  Ground,  Fort  Huachuca,  Arizona 
85613 

4.  PRINCIPAL  INVESTIGATOR.  Colin  M.  Giorgi,  Test  Operations  Office, 
STEEP-MT-T,  AUTOVON  879-6016/6417 

5.  STATEMENT  OF  THE  PROBLEM.  Optical  turbulence  in  the  atmosphere  can  signi¬ 
ficantly  degrade  the  field  performance  of  both  active  and  passive  optical  sys¬ 
tems.  Measurement  techniques  currently  being  used  to  characterize  this  uncon¬ 
trolled  and  highly  variable  test  parameter  are  inadequate.  A  measurement 
technique  is  required  which  will  quantitatively  and  accurately  account  for  the 
contribution  of  optical  turbulence  to  variances  in  equipment  performance  dur¬ 
ing  field  testing. 

6.  BACKGROUND .  The  effects  of  atmospheric  turbulence  on  the  performance  of 
imaging  and  beam  forming  optical  systems  have  been  well  known  for  many  years. 
Measured  variances  in  beam  intensity  and  beam  position  due  to  atmospheric 
turbulence  have  shown  that  the  field  performance  of  electro-optical  systems, 
and  active  systems  in  particular,  can  be  severely  biased  by  the  turbulence 
characteristics  of  the  optical  path.  Because  of  the  critical  time  and  speci¬ 
fic  path  sensitivity  of  the  turbulence  parameter,  existing  measurement  methods 
are  deficient  in  establishing  the  true  path  turbulence  structure  affecting  a 
test  items  performance.  During  a  recent  five  year  investigative  effort,  the 
Naval  Postgraduate  School  (NPS)  applied  optical  transfer  function  (OTF) 
measurement  technology  to  the  turbulence  measurement  problem  using  a  laser 
source  and  a  slit  scanning  telescope  to  collect  field  data.  The  results  of 
the  study  and  empirical  evaluation  indicate  that  the  technique  developed 
eliminates  the  deficiencies  of  previous  turbulence  measurement  methods.  The 
key  to  the  measurement  is  the  determination  of  a  properly  path-position 
weighted  Cn^ ,  the  turbulence  structure  constant  for  optical  index.  Dur¬ 
ing  the  review  of  the  investigative  effort  by  the  TECOM  Laser  Technical  Com¬ 
mittee,  it  was  determined  that  the  technique  had  the  potential  of  being  ap¬ 
plied  to  the  test  data  correction  problem  by  utilizing  the  video  imaging  and 
automated  data  analysis  instrumentation  currently  used  to  collect  spot  data 
for  laser  designators.  The  use  of  a  video  data  recording  system  and  an  auto¬ 
mated  video  data  reduction  capability  has  definite  advantages  over  the  mechan¬ 
ical  slit  scanning  device  used  in  the  NPS  investigation  but  application  of  the 
technique  will  primarily  depend  on  the  feasibility  of  utilizing  the  test  item 
laser  as  the  probe  source. 

7.  GOAL .  To  develop  an  improved  test  methodology  for  measuring  atmospheric 
turbulence  for  application  to  the  testing  of  electro-opt ical  systems  in  the 
field. 
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8.  DESCRIPTION  OF  INVESTIGATION. 


a.  This  investigation  will  develop  a  test  method  for  collecting  and  an¬ 
alyzing  the  modulation  transfer  function  of  atmospheric  paths  for  the  purpose 
of  deriving  the  turbulence  structure  constant  for  optical  index  (Cn^).  The 
thrust  of  the  effort  will  be  directed  at  implementing  the  NPS  image  scanning 
technique  by  applying  the  laser  spot  image  recording  and  analysis  instrumen¬ 
tation  currently  being  used  by  TECOM  in  support  of  laser  designator  testing. 
The  investigation  will  cover  both  a  theoretical  and  empirical  evaluation  of 
the  application  and  will  include  the  methodology  essential  to  applying  the 
Cn2  parameter  data  to  analyze  test  item  performance  data. 

b.  The  US  Army  Electronic  Proving  Ground  will  conduct  the  investigation 
in  two  phases  as  follows: 

(1)  Survey  and  Analysis  Phase. 

(a)  The  initial  effort  will  cover  a  survey  followup  to  the  NPS 
investigation  to  determine  what  additional  study  and  experimental  work  has 
been  accomplished  since  the  conclusion  of  the  basic  investigation. 

(b)  The  major  effort  during  the  analysis  phase  of  the  investi¬ 
gation  will  address  the  feasibility  of  applying  the  NPS,  MTF  image  analy¬ 
sis  technique  to  TECOM's  electro-optical  system's  test  programs  using  the 
laser  spot  image  recording  and  data  reduction  system.  The  priority  applica¬ 
tion  will  be  the  laser  designator  test  programs.  As  the  NPS  measurement  set¬ 
up  uses  a  one-way  path  and  a  reference  laser  source  which  is  configured  for 
the  measurement  requirement,  a  determination  will  be  made  with  regard  to  the 
feasibility  of  using  the  unmodified  test  item  laser  as  the  probe  source  and  a 
folded  optical  path. 

(c)  Following  a  determination  of  the  feasibility  of  applying 
the  NPS  measurement  technique  to  laser  designator  testing,  a  test  plan  for  the 
field  experimental  phase  of  the  investigation  will  be  prepared.  The  test 
planning  effort  will  include  the  identification  of  the  instrumentation,  real 
estate  and  personnel  resource  requirements  to  implement  the  field  tests. 

(2)  Empirical  Testing  Phase. 

(a)  The  purpose  of  this  phase  will  be  to  empirically  verify  the 
feasibility  of  the  technique  application.  The  initial  effort  will  require  the 
acquisition  of  instrumentation  support  as  specified  in  the  test  plan.  The 
testing  may  be  conducted  at  EPG  or  at  another  test  installation  depending  upon 
the  availability  of  representative  test  items  and  measurement  instrumentation. 
Regarding  the  latter,  EPG  is  planning  to  acquire  a  laser  spot  data  acquisition 
and  analysis  capability  by  the  end  of  FY81  in  order  to  support  the  RPV  test 
program. 
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(b)  The  data  acquisition  and  analysis  phase  will  include  eval¬ 
uating  data  from  both  the  experimental  setup  as  well  as  data  from  on-going 
laser  test  programs  on  a  non-interfering  basis. 

(c)  With  the  conclusion  of  the  field  experiments  and  subsequent 
data  analysis,  a  final  report  will  be  prepared  covering  the  complete  effort. 

If  feasibility  of  the  technique  is  established,  the  report  will  include  a  de¬ 
tailed  description  of  the  test  procedure  and  the  data  analysis  program  and 
recommendations  for  instrumentation  development  acquisition  as  appropriate. 

c.  Investigation  Schedule. 

MILE  STONE /ACTIVITY  SCHEDULE 


FY 

81 

(Qtrs) 

FY 

82 

(Qtrs) 

1 

2 

3  4 

1 

2 

3  4 

Analytical  Phase 

Updating  Survey 

X 

Error  Analysis 

X 

X 

Feasibility  Study 

X  X 

X 

Test  Design 

X 

X 

X 

Empirical  Evaluation  Phase 

Resource  Acquisition 

X 

X 

X 

Lab  and  Field  Testing 

X 

X 

X 

Data  Analysis 

X 

X 

X 

Report  Preparation 

X 

d.  This  investigation  will  result  in  the  development  of  an  improved 
optical  turbulence  test  procedure  and  the  determination  of  the  supporting  data 
acquisition  and  processing  instrumentation  and  software  requirements. 

e.  Environmental  Impact  Statement.  Execution  of  this  task  will  not  have 
an  adverse  impact  on  the  quality  of  the  environment. 

f.  Health  Hazard  Statement.  Execution  of  this  task  may  involve  health 
hazards  to  personnel.  This  task  may  require  project  personnel  to  participate 
in  field  tests  involving  the  operation  of  laser  equipment.  Safety  standards 
for  laser  field  testing  are  well  documented  in  TB  MED  279  titled  "Control  of 
Hazards  to  Health  from  Laser  Radiation".  All  field  experiments  will  be  con¬ 
ducted  in  full  compliance  with  the  regulation. 

9.  PROGRESS.  This  task  received  partial  funding  in  FY81  (35k).  The  updating 
survey  was  initiated  in  January  1981  with  a  followup  review  of  the  NPS 
investigations.  The  major  portion  of  the  analysis  and  test  design  phases  of 
the  investigation  are  being  performed  by  NPS  staff  scientists.  The  theoreti¬ 
cal  analysis  and  empirical  measurements  made  to  date  have  established  the 
basic  feasibility  of  replacing  a  telescope/mechanical  slit  scanner/detector 
unit  with  a  video  camera  and  a  video  digitizing  subsystem.  A  prototype 
Cn^  measurement  system  has  been  assembled  by  the  NPS  using  EPG  GFP. 

Subsystem  interfacing  has  been  completed  and  controlled  laboratory  tests  are 
underway. 
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10.  JUSTIFICATION. 

a.  Mission  and  Impact  Statements. 

(1)  Association  with  mission.  TECOM  and  EPG  are  assigned  the  mis¬ 
sion  of  conducting  the  development  testing  of  a  wide  range  of  electro-optical 
devices  and  systems.  These  include  both  active  (laser)  and  passive  (sensor) 
equipments . 

(2)  Present  capability,  limitations,  improvement,  and  impact  on  test 
if  not  performed  in  the  proposed  fiscal  year.  The  most  common  method  used  to 
describe  atmospheric  turbulence  is  simply  a  qualitative  assessment  of  the 
"seeing"  conditions  as  good,  fair  or  poor.  This  imprecise  and  subjective 
criteria  was  developed  primarily  as  a  general  measure  of  optical  attenuation 
and  it  does  not  account  for  the  time  and  specific  path  dependent  nature  of 
atmospheric  turbulence.  The  second  method  employs  point  temperature  measure¬ 
ments  which  also  do  not  provide  the  key  time  and  integrated  path  data 
essential  to  correlate  turbulence  effects.  Consequently,  for  most  EO  testing, 
the  turbulence  factor  remains  a  highly  variable  and  uncontrolled  test  condi¬ 
tion  which  can  produce  serious  equipment  performance  data  biasing.  The  appli¬ 
cation  of  the  NPS  technique  will  provide  a  turbulence  measurement  which  is 
time  and  optical  path  coincident.  This  coincidence  will  allow  the  turbulence 
variable  to  be  fully  monitored  and  accounted  for,  thereby  eliminating  the 
biasing  of  test  data.  If  this  task  is  not  completed  in  FY83,  the  capability 
will  not  be  available  to  support  the  RPV  laser  designator  test  scheduled  for 
FY83  or  other  key  TECOM  laser  tests  sched-uled  for  the  FY83-84  time  frame. 

b.  Dollar  Savings.  The  successful  implementation  of  a  Cn2  measure¬ 
ment  capability  is  expected  to  reduce  the  extent  of  field  testing  required  for 
electro-optical  devices  and  systems.  This  reduction  will  result  from  the 
ability  to  monitor  and  account  for  an  important  test  variable  which  now  re¬ 
quires  many  test  runs  and  testing  at  different  test  locations  to  establish  a 
meaningful  statistical  base.  The  successful  monitoring  of  the  variable  should 
also  lead  to  using  the  data  to  develop  performance  predictions  for  an  extended 
range  of  atmospheric  conditions.  This  modeling  capability  should  reduce  the 
requirement  for  some  field  testing  now  being  conducted.  However,  the  genera¬ 
tion  of  specific  dollar  saving  estimates  will  have  to  await  the  outcome  of  the 
feasibility  evaluation  and  the  determination  that  image  data  now  collected  on 
laser  tests  can  be  directly  used  for  the  turbulence  measurement. 

c.  Workload.  The  methodology  will  be  applied  initially  to  all  laser 
equipment  field  testing  to  include  designators,  range  finders,  radar  and 
illuminators.  If  successful,  the  techniques  will  also  be  applied  to  the  field 
testing  of  all  electro-optical  sensors. 
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Examples  of  items  anticipated  for  testing  at  EPG  include: 


Test  Item 


FY  83  84  85  86 


Chaparral  Night  Vision  Subsystem  X 

RPV  Designator  X 

RPV  Television  Sensor 
RPV  Mini-FLIR 
ELOCARS 


X 

X 

X 

X 


NOTE:  A  number  of  laser  and  passive  EO  sensor  tests,  scheduled  at  other  TECOM 
facilities  during  the  1983-86  time  frame,  will  also  utilize  the  developed 
methodology. 


d.  Recommended  TRMS  Priority.  Refer  to  the  workload  paragraph  (lO.c) 
and  the  ODCSOPS  priority  listing. 

e.  Association  with  Requirements  Documents.  The  ROC's  for  most  laser 
and  E-0  systems  specify  that  the  systems  will  perform  in  a  "degraded  atmos¬ 
phere"  environment  but  do  not  quantify  the  degradation  conditions  or  level. 
The  development  of  the  turbulence  measurement  methodology  can  be  an  important 
step  towards  quantifying  the  "degraded  atmosphere"  requirement. 


11.  RESOURCES. 


a.  Financial. 


( 1)  Funding  Breakdown. 

Dollars  (Thousands) 

FY  83 

In-House  Out-of-House 


Personnel  Compensation 

20,000 

Travel 

Consultants 

&  Other  Svcs 

3,000 

34,000 

Materials  & 

Supplies 

3,000 

Equipment 

Subtotals 

5,000 

31,000 

34,000 

FY  TOTALS 

65,000 

NOTE:  No  funds  were  provided  for  this  project  in  1983.  Funding  requirements 
and  funding  provided  were  as  follows: 

Requirement  Funding 
FY  81  $54,000  $35,000 
FY  82  $61,000  $29,400 
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Optical  Turbulence  Measurement  (Cont'd) 


(2)  Explanation  of  Cost  Categories. 


(a)  Personnel  Compensation.  Covers  in-house  labor  costs  for 
the  principal  investigator  and  other  in-house  project  support  personnel. 

(b)  Travel .  Travel  is  required  to  coordinate  the  task  with 
other  TECOM  installations.  Travel  in  FY83  may  be  revised  if  the  field  testing 
is  scheduled  for  WSMR  or  YPG. 

(c)  Consultants  and  Other  Services.  Consultant  support  will  be 
utilized  to  assist  in  performing  the  detailed  technical  feasibility  analysis 
and  to  evaluate  the  field  test  data. 


(d)  Materials  and  Supplies.  Recording  media  and  incidental 
supplies  will  be  required  to  support  the  field  testing. 

(e)  Equipment .  The  empirical  verification  field  tests  will  re¬ 
quire  a  designator  source.  Current  shortage  of  fielded  items  or  spare  test 
items  indicates  that  a  simulator  unit  will  have  to  be  procured.  FY83  funding 
covers  equipment  contingencies  including  servicing  and  repair. 


b.  Anticipated  Delays.  A  delay  in  EPG's  acquisition  of  the  programmed 
laser  spot  recording  and  data  reduction  instrumentation  capability  will  extend 
the  initiation  of  the  field  verification  testing.  However,  an  alternate  sup¬ 
port  capability  may  be  made  available  at  WS MR  or  YPG  to  alleviate  the  problem. 


c.  Obligation  Plan  (FY  83). 


FQ  1  2 _ 3 _ 4  TOTAL 

Obligation  Rate  3.0  38.0  14.0  10.0  65.0 

(Thousands) 


d.  In-House  Personnel. 

(1)  Requirements .  _ FY  83 _ 

_ Manhours _ 

Number  Required  Available 

Physicist,  GS-1310  1  500  500 

Elec  Engr ,  GS-0855  1  500  500 

Totals  1000  1000 


(2)  Resolution  of  Non-available  Personnel. 


Not  applicable. 
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Optical  Turbulence  Measurement  (Cont'd) 


12.  INVESTIGATION  SCHEDULE. 

FY  83 

ONDJFMAMJJAS 

In-house  -  -  -  -  -  -  -  -  R 

Contract  A  —  ——  —  —  — 

Symbols:  -  -Active  investigation  work  (all  categories) 

...Contract  monitoring  (in-house  only) 

A  Award  of  Contract 
R  Final  report  due  at  HQ,  TECOM. 

13.  ASSOCIATION  WITH  TOP  PROGRAM .  It  is  expected  that  the  turbulence 
measurement  methodology  initially  developed  for  specific  commodity  items  such 
as  laser  designators  and  rangefinders  will  be  incorporated  into  the  existing 
TOPs  covering  those  items.  As  more  generalized  procedures  are  developed  that 
apply  to  a  broad  class  of  E-0  test  items,  then  a  new  TOP  will  be  generated 
that  can  be  referenced  in  the  E-0  commodity  item  TOPs. 

FOR  THE  COMMANDER: 


GARY  L.  KOSMIDER 

Colonel,  SigC 

Director  of  Materiel  Test 
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APPENDIX  C 

VIDEO  DATA  ANALYSIS  USING  A  DIGITAL  VIDEO  PROCESSOR  SYSTEM 


C- 


ABSTRACT 


This  appendix  presents  a  computer  program  designed  to  analyze  the  data 
from  a  TV  camera  for  investigating  laser  beam  propagation  through  the  atmos¬ 
phere.  It  uses  aspects  of  Fourier  optical  theory  to  analyze  the  TV  image  to 
measure  the  effects  of  atmospheric  disturbances  and  platform  stability  on  the 
target  spot.  The  computer  program  is  written  for  a  data  analysis  configura¬ 
tion  consisting  of  a  video  recorder,  video  disk  unit,  digital  video  image 
processor  (Quantex  DS-30),  and  a  desk-top  computer  unit  (HP  9825). 
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I.  INTRODUCTION 


A.  PROBLEM 

A  problem  has  occasionally  been  experienced  in  concen¬ 
trating  laser  energy  from  a  laser  designator  onto  a  target 
and  sufficiently  illuminating  that  target  so  a  weapon 
sensitive  to  that  illumination  will  home  in  on  the  target. 
Two  circumstances  have  possibly  caused  this  to  happen;  the 
energy  is  either  too  diffused  at  the  target  to  illuminate  it 
properly  or  the  beam  partially  or  completely  misses  the 
target.  The  aim  of  this  study  is  to  develop  a  computer 
program  to  aid  in  determining  which  part  of  these  effects  is 
atmospheric  and  which  is  the  instability  of  the  laser  and 
stabilization  system  in  the  illuminating  aircraft. 

B.  METHOD  OF  APPROACH 

The  analysis  of  the  problem  requires  two  general  types 

of  systems.  One  system  evaluates  the  performance  of  the 

designator  using  some  techniques  from  Scott  [Ref.  1].  A 

silicon  vidicon  views  the  successive  pulsed  designator  spots 

on  the  target  from  a  location  near  the  target.  A  line 

spread  function  of  the  spot  on  the  target  and  the  standard 

deviation  of  the  wander  of  the  laser  beam  is  produced.  The 

2 

other  system  provides  measurement  of  Cn  for  atmospheric 
turbulence  along  an  optical  path  adjoining  the  designator 
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optical  path.  The  laser  on  the  target  is  viewed  by  a 

silicon  vidicon  located  near  the  designator's  optics.  The 
2 

value  of  Cn  obtained  is  inherently  properly  path  weighted 

l.0  express  the  atmospheric  effects  on  the  designator  spot. 

Techniques  from  Fourier  optics  theory  use  the  measured  value 
2 

of  Cn  to  predict  the  intensity  pattern  of  the  laser  beam  on 
a  distant  target.  The  results  from  the  two  systems  can  then 
be  analytically  compared  to  determine  the  causes  of  spot 
wander  and  broadening. 

The  sponsor-supplied  hardware  includes  an  HP-9825 
calculator  with  23K  bytes  of  internal  memory  for  equipment 
control  and  data  processing;  an  HP-9885  disc  memory  used  for 
additional  storage  space;  and  a  Quantex  DS-30  Digital  Video 
Analyzer  used  for  digitizing  the  analog  video  input  from  an 
Eigen  video  disc;  a  Panasonic  NV-1240  video  tape  recorder 
for  initial  recording  of  video  from  the  target  on  a  silicon 
vidicon  camera.  After  analysis,  output  is  in  the  form  of 
plots  produced  on  an  HP-9862  plotter. 

The  sequence  of  analysis  is  currently  as  follows. 
Modulation  transfer  functions  (MTF)  of  the  laser  output  and 
optics  are  measured  or  calculated  and  stored  on  disc.  A 
video  recording  of  the  target  being  illuminated  by  the 
designator  is  made  for  a  number  of  video  frames,  the  target 
alone  is  also  recorded,  for  later  subtraction  from  the  total 
picture  to  produce  an  image  of  the  laser  spot  alone.  The 
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video  recording  is  then  played  through  the  DS-30,  controlled 
by  the  HP-9825,  to  digitize  the  information  for  use  by  the 
HP-9825.  The  HP-9825  takes  the  digitized  video  of  the 
target,  produces  a  line  spread  function  (LSF)  of  the  image, 
subtracts  the  background,  and  averages  a  number  of  frames  to 
produce  a  short  term  measurement  of  the  laser  energy  on  the 
target.  The  amount  that  each  laser  spot  wanders  from  the 
others  is  also  stored.  In  the  second  phase,  the  HP-9825 
uses  the  MTF's  of  the  laser,  optics,  and  atmospheric 
measurements  to  produce  a  total  system  predicted  MTF,  LSF, 
and  wander.  These  results  of  the  target  spot  measurement 
and  the  baseline  studies  can  be  compared  for  possible 
correlation  . 

The  laser  pulse  is  timed  so  that  it  occurs  on  the  fly¬ 
back  of  the  TV  signal.  This  will  cause  the  laser  spot  to 
appear  on  one  field  of  a  frame.  For  accurate  results  the 
DS-30  needs  to  digitize  a  single  field  for  each  pulse  of  the 
laser.  The  videotape  unit  currently  used  will  not  allow 
this  capability.  For  this  reason  the  tape  recording  video 
is  transferred  to  an  Eigen  video  disc ,  which  has  the 
capability  of  displaying  each  frame  or  field  (1/2  video 
frame  interlace)  individually  and  of  selecting  any  one  frame 
or  field  for  viewing. 
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II.  BACKGROUND  AND  PROGRAM  DESCRIPTION 


A.  THEORETICAL  DESCRIPTION 

In  order  to  combine  the  individual  optical  components  of 
the  laser  to  the  target  system,  some  elements  of  Fourier 
theory  must  first  be  remembered.  The  Fourier  transform 
represents  the  one  dimensional  position  variable  g(y) 
expressed  in  the  spatial  frequency  domain  U(v).  The  inverse 
transform  repeats  the  operation  in  the  opposite  direction. 
Their  forms  are  commonly  represented  as  follows: 


U(v)s 

g(y)  = 


g ( y ) exp( -2 

U(v)exp(2 


ir 


ir  ivy)dy= 


ivy)dv  = 


^Lg(y)] 

«^'1[U(v)] 


The  program  calculates  the  value  of  the  integral  at  a  preset 
number  of  points  and  yields  a  discrete  Fourier  transform.  A 
theorem  from  Fourier  theory  that  allows  us  to  compute  the 
line  spread  function  as  presented  in  Fried  [Ref.  2]  for  one 
of  the  elements  in  the  system  given  the  others  is  the  con¬ 
volution  theorem  stated  as  follows: 

If  &  Cg(y) ]  =  G ( v ) 
and  &  C h ( y )  ]  =  H(v) 
then  &  [g*h]  =  GH 

(*  -  convoluted  with) 
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The  first  element  in  the  system  we  are  analyzing  is,  of 

course,  the  laser  itself.  The  radial  intensity  output  of 

o  2  2 

the  laser  is  generally  Gaussian  of  the  form  A=A  exp(-r  /2 a  ) 
as  shown  in  Figure  1 .  This  can  be  calculated  for  each  laser 
for  approximate  results  but  was  usually  measured  directly  in 
this  study. 

The  above  radial  distribution  of  intensity  must  then  be 
integrated  over  x  to  produce  a  one  dimensional  line  spread 
function  for  the  source.  A  Fourier  transform  is  then  used 
on  this  LSF  to  yield  the  MTF  of  the  laser.  The  MTF  is  a 
measure  of  the  spatial  frequency  response  of  a  system 
compared  to  the  input.  A  "perfect”  system  response  would  be 
1.0  out  to  its  limit  (large  spatial  frequencies)  and  then 
dropping  to  zero. 

The  optics  associated  with  the  laser  is  basically 
diffraction  limited.  The  "Airy  function"  is  used  to  cal¬ 
culate  the  diffraction  point  spread  function  for  the  laser 
optics.  This  point  spread  function  is  then  converted  to  a 
line  spread  function,  Fourier  transformed  to  get  the  MTF  and 
then  combined  by  the  convolution  theorem  into  a  laser  system 
Fourier  transform. 

The  next  step  is  to  calculate  the  effects  of  atmospheric 

turbulence  on  the  laser  beam.  The  value  of  C  ,  "the  index 

n 

of  refraction  turbulence  structure  constant,"  as  expressed 
by  the  relationship  from  Tatarski  [Ref.  31  and  Ochs  et.  al. 
[Ref.  4]: 
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FIGURE  1.  IDEAL  LASER  INTENSITY  OUTPUT  LINE  SPREAD  FUNCTION 


<(Trl"Tr2)2  >=r2/3cT2  and  Cn =79x 1 0‘6p/T2CT 

where  T  , -T  _  is  absolute  temperature  at  two  points 
r  l  r  2 

separated  by  a  distance  r  and  p  is  pressure  in  millibars;  is 
used  to  predict  the  MTF  of  the  atmosphere  over  the  path 
between  the  laser  and  the  target.  Measurements  made  along 
the  path  are  inputs  to  the  Cn  equation  which  predicts  the 
turbulence  effects  on  the  MTF  and  wander  of  the  beam.  In 
the  field,  however,  plans  are  to  measure  Cn  by  optical 
methods  decribed  in  Crittenden  [Ref.  5]. 

The  theory  then,  according  to  Fried  [Ref.  2],  shows  that 
the  total  effect  of  the  atmosphere  and  the  laser  can  be 
found  by  the  product  of  the  two  transform  functions: 

G  total  (v)=G  optics(v)  x  G  atra(v) 

The  inverse  transform  then  will  yield  the  line  spread 
function  of  the  total  system,  the  LSF  as  predicted  on  the 
target . 

Crittenden  et .  al ,  [Ref.  53  describe  the  process  of 
converting  the  one  dimensional  LSF  to  a  circular  symmetric 
PSF  using  the  Abel  transform,  from  Griem  [Ref.  63.  The 
power  inside  a  circle  of  radius  R  can  then  be  obtained  by 
integrating  the  PSF  out  to  R.  For  the  purposes  of  compari¬ 
son  in  this  study,  however,  the  resulting  system  LSF  is  used 
to  compare  the  prediction  with  measured  values. 
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B.  HARDWARE  USE  AND  LIMITATIONS 

The  HP-9825  calculator  used  for  this  program  had  an 
internal  capacity  of  23K  bytes.  This  was  quite  adequate  for 
program  storage  but  because  of  the  great  volume  of  digitized 
video  data  that  was  involved,  the  inclusion  of  an  HP-9885 
disc  drive  system  was  necessary.  The  Quantex  DS-30  can 
store  and  rapidly  transfer  a  maximum  of  512  lines  of  video 
with  512  pixels  (picture  elements)  per  line,  easily 
exceeding  the  HP-9825’s  memory  capability.  The  large 
overflow  of  data  was  reduced  by  storing  each  image  into  a 
quarter  of  the  memory,  thus  limiting  the  amount  of  data  that 
needed  to  be  transferred.  This  also  made  disc  storage  more 
reasonable  and  sped  up  program  running  time. 

The  DS-30  was  also  capable  of  taking  the  difference  on  a 
pixel  by  pixel  basis  between  a  reference  image  and  the  input 
and  storing  that  in  memory.  Differencing  was  necessary  in 
order  to  remove  the  background  from  the  recorded  image  for 
processing  the  laser  spot  alone.  Normally,  it  is  necessary 
to  divide  the  image  by  the  background  to  offset  the  effect 
of  the  non-uniform  reflectivity  of  the  target.  This  was  not 
possible  with  the  DS-30  so  these  steps  need  to  be  performed 
by  the  computer.  The  division  was  not  possible  to  do  at 
this  time  because  of  the  memory  limitations  of  the  HP-9825. 
This  must  be  done  pixel  by  pixel  since  division  by  the  LSF 
is  not  equivalent.  The  memory  limitation  was  alsc  he 
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reason  why  the  HP-9825  used  the  LSF  instead  of  pixel  by 
pixel  subtraction  for  this  process. 

All  the  peripheral  equipment  to  the  HP-9825  was 
controlled  through  an  IEEE  standard  8  bit  interface  bus. 

This  allowed  control  of  all  aspects  of  data  acquisition  and 
processing  to  be  modified  by  software.  The  video  recorders 
were  the  exception  to  this  and  provided  the  only  real  manual 
manipulations  required  after  the  program  had  begun.  Each 
frame  of  video  needed  by  the  program  had  to  be  selected  when 
requested  by  the  program  because  of  the  lack  of  interface 
control . 

The  linearity  of  the  recording  system  was  checked  and 
found  to  be  almost  distortion  free.  A  signal  from  an  image 
was  passed  directly  to  an  oscilloscope  from  the  vidicon  and 
compared  side  by  side  to  a  signal  from  the  same  image  which 
had  been  recorded  in  the  tape  recorder ,  transferred  to  the 
video  disc  and  then  passed  through  the  DS-30  circuitry  to  a 
digital  analog  converter  for  viewing.  Crittenden  et.  al . 
[Ref.  71  present  photographs  of  the  two  images  for 
comparison  and  further  discussion. 

Variable  usage  became  a  problem  as  the  program  grew, 
which  wa3  due  to  the  limitations  of  the  HP-9825.  The 
HP-9825  has  available  twenty-six  variables  (A-Z),  twenty-six 
variables  for  arrays  (A-Z),  a  subscripted  variable  (r)  with 
as  many  elements  as  the  memory  has  space  for  plus  a 
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subprogram  variable  (p)  for  use  in  passing  parameters 
through  "called"  subprograms.  One  of  the  usage  problems  -as 
using  identifiable  variables  in  passing  from  main  programs 
to  subprograms.  A  good  number  of  the  subprograms  originally 
were  designed  for  other  programs  and  the  variables  often 
conflicted.  The  r  variable  is  nice  to  use  here  because  of 
its  practically  unlimited  numbers.  A  more  serious  problem, 
though,  was  the  inability  of  HPL ,  Hewlett  Packard's  version 
of  the  Basic  language,  to  allow  an  entire  array  to  be  used 
as  a  parameter  to  pass  to  a  subprogram.  This  at  first  made 
necessary  extra  internal  storage  for  additional  arrays  to  be 
used  only  within  a  certain  subprogram,  with  the  variable 
name  being  changed  through  a  lengthy  process  upon  entering 
or  leaving  that  subprogram.  A  much  faster  and  easier  way 
was  using  separate  storage  files  on  the  disc  for  each  major 
array  manipulation.  When  entry  into  a  subprogram  was  needed 
the  array  could  be  stored  in  the  file  assigned  to  that 
subprogram,  which  could  then  be  read  out  and  labelled  as  any 
array  needed  for  each  subprogram. 

The  programming  to  use  the  DS-30  also  was  primarily  a 
language  problem.  The  DS-30  was  very  sensitive  to  the 
format  in  which  it  received  data  or  commands.  All  leading 
zeroes  and  carriage  return/line  feed  commands  had  to  be 
suppressed.  The  address  of  command  locations  within  the 
DS-30  also  had  to  be  encoded  in  hexadecimal  nibbles  (one- 
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half  byte)  and  sent  in  reverse  order  to  load  the  registers 
properly.  Commands  to  the  DS-30  were  sent  as  ASCII 
character  strings.  Data  was  output  from  the  DS-30  to  the 
HP-9825  ’'ia  a  fast  read  buffer  which  allowed  a  very  fast 
data  transfer  rate.  In  use,  a  block  of  256  bytes  at  a  time 
were  transferred  representing  one  horizontal  TV  line  of  the 
quartered  memory  space.  This  was  done  in  a  loop  which  went 
back  and  took  256  lines,  completely  transferring  the  image. 

It  was  decided  that  output  would  be  on  the  HP-9825 
plotter  so  that  direct  visual  comparisons  could  be  made 
between  measured  and  predicted  intensity  LSF's.  This  was 
much  more  dramatic  and  easier  to  see  than  tabular  output 
which  then  would  have  to  be  compared  on  a  point  by  poiwt 
basis  for  a  correlation. 


15 


1 


III.  SUMMARY 


This  program  can  be  useful  in  solving  the  problem  tnat 
spurred  this  study.  It  will  be  able  to  show  how  much  wander 
can  be  expected  out  of  a  laser  designator  beam  as  it  travels 
through  the  atmosphere  from  the  laser  to  the  target.  Based 
on  the  video  recording  of  the  actual  image  on  a  target  it 
will  show  comparatively  the  amount  of  wander  actually 
present  in  the  laser.  Plans  are  to  unfold  the  results 
analytically,  but  at  present,  a  visual  comparison  of  the 
measured  and  predicted  patterns  is  made  to  determine  if 
there  is  any  excess  wander  or  "jitter"  attributable  to  the 
laser  platform  itself.  This  determination  will  be  helpful 
in  setting  specifications  for  the  stability  needed  in  the 
system  to  remove  as  much  excessive  spot  movement  as 
possible.  It  will  also  allow  testing  to  determine  if  those 
specifications  are  being  met. 

This  program  can  also  be  used  to  verify  atmospheric 
problems  when  both  the  laser  and  target  spot  are  fixed  and 
known.  Essentially,  any  of  the  laser  to  target  elements  can 
be  determined  by  knowing  the  others  and  calculating  tne 
unknown  using  adaptations  of  the  program.  The  beginning 
part  of  this  program  can  be  used  to  digitize  and  store  any 
video  image  and  produce  an  LSF  of  the  image.  This  may  be 
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useful  in  recording  and  testing  the  output  of  several 
different  lasers  and  optical  aperatures  for  comparative 
analysis . 

This  computer  program  was  designed  in  the  initial  stages 
of  a  continuing  project  studying  the  described  problem.  If 
desired,  changes  can  be  made  to  the  program  to  adapt  to 
changing  field  conditions  and  analysis  requirements.  In 
particular,  the  program  now  requires  a  greater  understanding 
of  its  internal  operations  and  operator  interaction  than  is 
necessarily  needed.  More  complete  cueing  prompts  and  data 
entry  parameters  may  be  worthwhile  changes.  Adaptability 
to  other  hardware  may  provide  more  program  control,  especi¬ 
ally  if  control  of  the  video  disc  or  tape  can  be  handled 
through  an  interface  bus  instead  of  manually.  Also, 
improvements  in  output  design  can  be  foreseen  to  expand  its 
usefulness  and  adaptability. 

Finally,  the  speed  of  program  operation  could  be 
increased  through  use  of  a  faster  computer  with  a  greater 
memory  capacity.  The  HP-9825  is  comparatively  slow  and  has 
insufficient  memory.  The  use  of  a  faster  system  such  as  the 
HP-1000  is  recommended  for  actual  measurements. 
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ANNEX  A 


USING 


:0MPU7I3  PROGRAM 


The  analysis  program  requires  considerable  preliminary 
set  up  and  operator  interaction.  3efore  the  program  is 
started  there  should  be  several  items  already  recorded  on 
videotape.  The  first  item  is  a  TV  camera  shot  of  the  laser 
output.  This  is  normally  taken  with  a  long  focal  length 
lens  and  filters  to  reduce  the  intensity  without  losing  the 
beam  pattern.  A  part  of  this  recording  should  include  a 
known  diffraction  grating  through  which  the  beam  passes  in 
order  to  measure  the  scale  of  the  data  for  each  system  of 
lenses  used.  The  diffraction  pattern  should  extend  in  the 
vertical  as  seen  on  a  TV  screen  for  proper  output  on  the 
plotter.  Also  needed  to  be  recorded  is,  of  course,  a 
sequence  of  laser  spot  images  on  the  target.  Part  of  this 
recording  should  have  the  target  only  as  a  background 
reference  for  later  calculations. 

The  next  step  that  needs  to  be  taken  is  to  compute  and 
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record  the  square  atmospheric  turbulence  index  (Cn  ),  then 
record  the  following  system  data  for  later  program  entry; 
wavelength  of  the  laser,  diameter  of  the  optics  objective 
lens,  ratio  of  the  obscuration  to  the  diameter  of  the 
objective  lens,  distance  to  the  target,  and  the  extinction 
coefficient . 
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The  first  cue  the  program  asks  is  for  the  input/output 
parameters.  Table  1  lists  the  parameters  and  the  different 
function  each  one  is  used  for. 

If  the  laser  source  pattern  is  not  to  be  computed  the 
next  cue  requests  that  the  laser  pattern  be  displayed  from 
the  video  disc  for  transfer  to  HP-9825  memory.  A  cue  is 
also  presented  asking  for  the  horizontal  width  of  the  images 
to  be  taken  from  the  DS-30.  The  first  value  should  be  the 
left  most  pixel  desired,  the  second  being  the  furthest  right 
pixel  to  be  used.  The  first  value  must  be  less  than  the 
second,  with  the  first  no  less  than  1  and  the  second  no  more 
than  256.  This  cue  is  also  presented  below  before  either 
the  scale  image  or  the  target  image  is  taken. 

If  the  scale  of  the  source  data  is  to  be  determined  the 
next  cue  requests  that  a  frame  be  displayed  on  the  video 
disc  of  the  laser  output  with  the  diffraction  grating  in 
place.  When  this  has  been  done  continue  the  program  and  it 
will  plot  out  an  LSF  of  that  image  whereby  the  scale  of  the 
data  can  be  determined  in  raicroradians  per  point. 

The  next  cue  alerts  the  operator  to  ensure  that  a  frame 
of  background  video  is  displayed  from  the  video  disc  for  use 
in  the  difference  process. 

The  program  will  then  ask  for  the  number  of  total  frames 
that  the  user  wishes  to  average  for  the  measured  target  spot 
LSF.  The  value  should  be  in  the  range  from  1  to  600.  The 
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latter  value  being  the  greatest  number  of  frames  the  video 
disc  can  hold  at  one  time.  The  program  will  then  proceed  to 
record  and  process  images  from  the  video  disc  and  DS-30* 

The  program  will  provide  a  cue  when  it  wants  the  next  image 
to  be  sequenced  on  the  video  disc. 

The  program  will  alert  the  operator  when  it  is  about  to 
produce  a  plot.  If  the  HP-9862  plotter  is  not  ready  simply 
stop  the  program  until  a  piece  of  plotting  paper  is  in 
place . 

The  program  will  request  the  system  data  after  the 
measured  image  spot  data  has  been  plotted.  The  data  format 
is  not  critical  but  should  be  in  the  correct  units  as 
follows : 

Wavelength 

Diameter  of  objective  lens 

Obscur at ion /objective 

Scale  of  data 

Range  to  target 

cn2 

Extinction  coefficient 

If  the  source  pattern  is  to  be  computed  the  program 
requests  peak  amplitude  and  standard  deviation  of  the 
intensity  distribution. 
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The  program  will  then  continue  computations  without 
further  inputs  until  it  has  run  through  with  the  initial 
parameters.  The  I/O  parameters  will  determine  at  which 
steps  output  is  plotted.  After  one  run  through,  the  program 

p 

will  request  a  new  Cn  ,  range,  and  coefficient.  If  the 
input  range  is  positive  the  program  will  again  run  through 
from  the  atmospheric  calculations  with  the  new  range.  If 
the  range  is  negative,  the  branch  is  to  the  beginning  of  the 
prediction  phase  where  all  new  system  inputs  are  requested. 
If  the  range  input  is  zero,  the  program  ends. 
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ANNEX  B 


PROGRAM  VERIFICATION  AND  FLOWCHARTS 
I.  MEASUREMENT  PHASE 

Section  I  (steps  0  -  50)  consisted  of  program  setup 
information  and  the  DS-30  function  command  programming. 

This  section  was  easily  verified  by  observing  that  all 
initializing  functions  were  accomplished  and  the  response 
of  the  DS-30  could  be  visually  noted. 

Section  II  (steps  51  -  86)  consisted  of  transferring 
the  contents  of  the  four  quarters  of  the  DS-30' s  memory  to 
the  HP-9825  and  producing  an  LSF  of  each  image.  For 
testing,  an  image  consisting  of  all  black  (8  bit  binary 
ones)  or  all  white  (8  bit  binary  zeroes)  was  placed  in  each 
of  the  four  quarters  of  memory,  each  quadrant  was  then  read 
out  to  the  HP-9825.  The  values  received  by  the  HP-9825 
corresponded  to  the  contents  of  the  DS-30  memory  quadrant 
and  did  show  proper  changes  when  the  subject  quadrant 
memory  was  inverted  from  black  to  white.  The  LSF  was 
produced  by  averaging  the  pixels  in  each  horizontal  line  of 
video,  each  line  producing  a  point  value  as  the  image  was 
scanned  vertically  from  top  to  bottom.  Examples  of  LSF ' s 
produced  for  a  "black”  and  a  "white"  image  are  in  Figures  2 
and  3* 
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FIGURE  2.  "BLACK"  IMAGE  LSF 
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FIGURE  3.  "WHITE"  IMAGE  LSF 


Section  III  (steps  87  -  118)  subtracted  a  reference 
background  from  an  image  to  yield  a  difference  image,  then 
shifted  the  resulting  LSF  so  that  the  center  of  area  under 
each  image  was  centered  on  zero.  In  this  way  the  wander 
from  each  succeeding  image  could  be  determined,  removed, 
and  the  total  frames  used  could  then  be  averaged  to  show 
the  short  term  atmospheric  action  on  the  laser  spot  image. 
Proper  operation  of  this  section  can  be  seen  in  Figure  4 
where  two  images  of  slightly  different  position  and 
intensity  were  averaged  to  yield  an  LSF  halfway  between  the 
two  image  functions  as  expected.  Figure  5  shows  a  differ¬ 
enced  LSF  where  a  small  laser  spot  was  added  to  a  reference 
background.  The  result  is  almost  entirely  the  laser  spot 
with  a  small  amount  of  "video  noise"  caused  by  the  camera 
and  associated  electronics. 
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FIGURE  4.  AVERAGING  EXAMPLE 


FIGURE  5.  LASER  SPOT  LSF 


FIGURE  6A.  MEASUREMENT  PHASE  FL0V7CHART 
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II.  PREDICTION  PHASE 

Section  IV  (steps  119  -  139)  enters  the  parameters  of 
the  atmosphere  and  optics  needed  to  model  and  predict  the 
laser  spot  pattern  on  the  target.  This  section  also  yields 
the  calculated  Gaussian  laser  intensity  output  if  that  data 
is  not  measured  in  the  next  section.  The  output  of  this 
section  is  shown  in  Figure  7. 

Section  V  (steps  140  -  172)  consists  mainly  of 
branching,  data  reading,  and  subroutine  action.  These  sub¬ 
routines,  branching  and  data  input  will  be  verified  with 
proper  operation  of  the  program  as  a  whole. 

Section  VI  (steps  173  -  183)  calculates  the  diffraction 
limited  optics  effect  on  the  laser  output  and  produces  a 
radial  point  spread  function  of  the  intensity.  The  output 
for  a  point  source  should  be  an  airy  function  as  results 
show  in  Figure  8. 

Section  VII  (steps  184  -  204)  again  consists  mainly  of 
subprograms  and  also  calculates  the  product  of  the  calcu¬ 
lated  optics  and  laser  output  transforms.  The  output  is  in 
Figure  9. 

Section  VIII  (steps  205  -  220)  computes  the  MTF  of  the 
atmosphere  based  on  the  inputs  from  Section  IV.  The  output 
is  in  Figure  10. 

Section  IX  (steps  221  to  the  end  of  the  main  program) 
is  mainly  subprogram  action  in  computing  the  final  LSF 
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predicted  on  the  target  using  the  source,  optics  and 
atmospheric  MTF's.  It  also  provides  a  link  to  re-enter  the 
program  for  changed  parameters.  The  output  of  this  section 
is  listed  under  the  subprograms  used. 
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FIGURE  7.  EXAMPLE  OF  CALCULATED  LASER  OUTPUT 


FIGURE  9.  LASER  AND  OPTICS  MTF 
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FIGURE  113.  PREDICTION  PHASE  FLOWCEAR' 
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FIGURE  I1C. 


PREDICTION  PHASE  FLOWCHART 


III.  SUBPROGRAMS 

Section  X  (subprogram  ABEL)  produces  a  point  spread 
function  (PSF)  from  the  LSF  computed  for  the  target  laser 
spot.  The  input  and  resulting  output  is  shown  in  Figures 
12  and  13- 

Section  XI  (subprogram  LSF)  enables  the  program  to 
compute  a  line  spread  function  for  later  conversion  to  a 
MTF  when  the  input  is  calculated  by  Sections  IV  and  VI  as 
point  spread  functions  or  when  the  data  is  input  in 
Section  V  in  the  point  spread  function  format.  Sample 
input  and  output  is  in  Figures  14  and  15. 

Section  XII  (subprogram  B)  calculates  the  fraction  of 
the  power  inside  a  circle  of  a  chosen  radius  about  the 
center  of  the  laser  spot.  This  is  helpful  in  determining 
if  the  spot  is  intense  enough  on  the  target.  The  input  and 
output  are  Figures  16  and  17. 

Section  XIII  (subprogram  FXFORM)  is  a  discrete  Fourier 
transform  program  that  functions  both  as  direct  and  inverse 
depending  on  the  flag  set  in  the  main  programs.  Several 
Fourier  pairs  were  checked  to  ensure  correct  calculations, 
Figures  18  through  21. 

Section  XIV  (subprogram  SOURCE)  digitizes  an  im'-g"*  of 
the  laser  source  for  inclusion  in  the  prediction 
calculations.  It  is  a  storage  of  DS-30  video  from  the 
Vidicon  camera  input. 
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Section  XV  (subprogram  SCALE)  digitizes  an  image  with  a 
standard  diffraction  grating  for  calculating  the  scale  of 
the  data.  Figure  22  is  an  example  of  this  image. 

Section  XVI  (subprogram  OUTPUT)  yields  in  useable  form 
all  the  data  required  for  analysis.  This  subprogram  was 
used  for  all  the  output  in  the  plot  format.  Functions 
"MIN”  and  "MAX"  are  used  by  the  "OUTPUT"  subprogram. 

The  program  was  run  as  a  whole  numerous  times  with 
different  I/O  parameters  and  data  in  order  to  exercise  the 
branching  combinations  and  program  steps.  Table  1  lists 
the  I/O  parameters  available  and  their  uses.  This 
operation  was  simply  to  ensure  the  program  ran  as  was 
intended  . 
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FIGURE  13.  "ABEL"  OUTPUT  (PSF) 
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FIGURE  14.  "LSF"  INPUT  EXAMPLE 
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FIGURE  15.  "LSF"  OUTPUT 
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FIGURE  22.  SCALE  IMAGE  EXAMPLE 


TABLE  1 


INPUT/OUTPUT  PARAMETERS 


ARRAY  VARIABLE 

VALUES 

USE 

0(1) 

0 

Same  as  1  except  program  after  one 
time  through  ends,  not  allowingpfor 
computations  with  new  range,  C  , 
ext.  coeff.  n 

1 

Source  pattern  input  is  point 
spread  function. 

2 

Source  pattern  input  is  line  spread 
function . 

3 

Source  pattern  input  is  modulation 
transfer  function  of  laser  only. 

4 

Source  pattern  input  is  modulation 
transfer  function  of  source  and 
optics . 

5 

Source  pattern  is  to  be  calculated 
with  amplitude  and  sigma  of  inten¬ 
sity  distribution  as  factors. 

0(2) 

0 

Optics  diffraction  input  as  PSF. 

1 

Optics  diffraction  input  as  LSF . 

2 

Optics  diffraction  pattern  to  be 
calculated  as  PSF. 

0(3) 

0 

No  plot  of  source  line  spread 
function  is  required. 

>0 

Plots  source  line  spread  function 
(available  only  if  Q(1)=0,1  or  5). 

0(4) 

0 

No  plot  of  source  MTF  required. 

>0 

Plots  source  MTF  (available  only  if 
0(1)  does  not  equal  3  or  4). 

0(5) 

0 

No  plot  of  product  of  source  and 
optics  MTF  required. 

>0 

Plots  product  of  source  and  optics 
MTF ' s  (available  only  if  0(1)  does 
not  equal  4 ) . 
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TABLE  1  (CONT’D) 


ARRAY  VARIABLE 

VALUE 

USE 

0(6) 

0 

No  plot  of  optics  diffraction  PSF 
required . 

>0 

Plots  optics  diffraction  PSF 
(available  only  if  0(1)  does  not 
equal  4  and  0(2)=2) . 

0(7) 

0 

No  plot  of  optics  LSF  required. 

>0 

Plots  optics  LSF  (available  only  if 
0(1)  does  not  equal  4  and  0(2)  does 
not  equal  1). 

0(8) 

0 

No  plot  of  optics  MTF  required. 

>0 

Plots  optics  MTF  (available  only  if 
0(1)  does  not  equal  4). 

0(9) 

0 

No  plot  of  atmosphere  MTF  required. 

>0 

Plots  atmosphere  MTF. 

0(10) 

0 

No  plot  of  system  MTF. 

>0 

Plots  total  system  MTF . 

0(11) 

0 

No  plot  of  LSF  of  predicted  target 
spot . 

>0 

Plots  LSF  of  target  spot. 

0(12) 

0 

No  plot  of  target  PSF  required. 

>0 

Plots  target  PSF. 

0(13) 

0 

No  plot  of  the  fraction  of  power 
inside  circle  required. 

>0 

Plots  fraction  of  power  inside  a 
circle  of  radius  r. 

0(14) 

0 

any 

other 

value 

No  plot  of  calculated  source  PSF 
required . 

Plots  calculated  source  PSF. 
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TABLE  1  (CONT'D) 


ARRAY  VARIABLE 

VALUES 

USE 

0(15) 

0 

Plots  measured  LSF  of  target  spot. 

>0 

Used  as  indicators  of  output 
branching  only,  entry  value  not 
significant. 

0(16) 

1 

Inputs  source  with  diffraction 
grating  from  DS-30  for  scale  compu¬ 
tations  from  output.  No  output 
plotted. 

2 

Same  as  1  except  plots  the  above. 

any 

other  Does  not  do  above  input 
value 


ANNEX  C 


PROGRAM  LISTING 


0  :  dim  BS123]  ,C$[51  ,D$[71  fIS[272I  ,0$(5J,PSf5]  ,wR  (  .5  J ,  X$.  [  40  j  ,R( 
5121  ,1.(512] 

1:  du  0(16]  ,  AS,  ( 13]  ,  Y$.  (  4G]  ,25.(40] 

2:  files  0  rT  ,£  ,R  ,G  ,V  ,  W,  C,  I  ,r 

.3:  fmt  1,  zf  c 

4:  time  150 

5:  cfg  3;cfg  2;cfg  1 

6:  rread  7,l;sprt  7,R(*];rread  7,1 

7:  bur  "  DATA  ",  I  $ ,  3 

8:  rread  4,l;sprt  4,R[*] 

9:  beep 

10:  "INPUT  I/O  CONTROL  0(1]": 

11:  for  1*1  to  16;ent  0[I];next  I 
12:  if  0(1] *2;gso  "SOURCE" 

13:  clr  70  4 

14:  if  O [ lo 1 *1  or  0(16]*2;gsb  "SCALE" 

15:  beepjdsp  "ENSURE  VIDEODISC  ON  .THEN  PRESS  CONTINUE" 

16:  stp 

17:  beepyent  "  #  OF  FRAMES  TO  AVERAGE'*,  H 

18:  ent  "WINDOW  WIDTH  , FIRST  PIXEL" , r9 ," LAST  PIXEL", rlO 
19:  H+1*S 
20:  0-c3 

21:  for  0*1  to  H  by  4 
22:  for  1*2  to  9 

23:  char  (  47 )  schar  (5  5)  ichar.  (64)  icnac (80)  &  char.  (96)  *C$ 

24:  if  I*2;"93Q  U20"*D$ 

25:  if  1*3;"  020 

26:  if  1*4; "930  120 "-0$ 

27:  if  1*5  ;  "020  "♦DS 
28:  if  1*6;  "930  22G"*D$ 

29:  if  1*7;  "  020  *-*D$ 

30  :  if  1*8;  "930  320"*D$ 

31:  if  1*9; "020 "-0$ 

32:  if  len  (D$)  <4;  gto  38 

33:  char  (1)  &"2"&char (2)  &C$&char  (15)  &DS[5  ,"]*£$ 

34:  cli  7 
35:  wait  100 
36:  wrt  704,  B$ 

37:  cli  7 

38:  cnar  (46)  *C$.[1, 1] 

39:  cnar(l)&"2  "&cdat  (2)  &C$&cnar  ( 15)  S.DS  [  1 , 3 } -AS 
40:  cli  7 
41:  wait  100 
42:  wrt  704, A$ 

43:  cli  7 

44:  if  I#3  or  0>l;gto  47 
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45 

46 

« 

49 

50 

51 

52 

53 

54 

55 

56 

57 
53 

59 

60 
61 
62 

63 

64 
65: 
66 

67 

68 

69 

70 

71 

72 

73 

74 
7  5: 

76 

77 
7b 

79 

80 
31 
82 

83 

84 

85 

86 

87 

88 
89: 

90 

91 

92 

93 

94 


dsp  “BACKGROUND  RECORDED" 

wait  500;dsp  "CONTINUE  WHEN  NXT  FRAME  READY";  Stp 
if  D$# " 0 20 "  or  1  =  3  or  1*9  ana  0<5;gto  T9  * 

dsp  t,CONT.  rtHEN  NXT  F RAME  SELECTED"  ;  stp 
next  I 

if  ri-0+l<4;gto  53 
for  X=1  to  4 
gto  54 

for  X*1  to  H-Q+l 

"!  "&"G“&"A"&"P"&char.  <96)  +  P$ 

if  X#l;gto  57 

l*J;0*Y;gto  62 

if  Xff 2 ; gto  59 

sf  g  3 ;  2  5  6  ♦  J ;  0*  Y ;  g  to  62 

if  X#3;gto  61 

1«>J  ;  255-*Y  ;  gto  62 

256«-J;255*Y 

for  L*1  to  256 

dsp  L 

0*D; "  "♦0-$;  Y+1*Y 
Y*512+J-S 
for  1*2  to  6 
int  (S/16)  +T 
T-U 

S-T*16+I*16*f 

U-S 

0$&char  (.T)  -0$ 
next  I 

char  (1)  &"*"&char  (34)  &char.  (2)  sOS&char.  (3)  & P$&char.  (15)  *b 5 
fat  1,  z,  c 
cli  7 

wrt  704 . 1 , B$ 

buf  "DATA" 

tfr  704, "DATA", 256 

rds ( "DATA") »E; if  E»-l;jrap  0 

for  I*r9  to  rlO 

num(  I$.[I,I  ] )  *Z 

D+Z*D 

next  I 

D/(rlQ-r9) *R[L] 
next  L 

if  fig 3*0; gto  91 

"TAKE  DIFFERENCE  OF  IMAGE  AND  BACKGROUND": 
rread  3,l;sread  3,I{*J;for  1*1  to  256 
abs(R(I]-I(Ij) ♦Rfl] ; nex  t  I 
gto  93 

rread  3,l;sprt  3,R[*J 

gto  113 

0*G;0*F 

for  1*30  to  226 
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95:  R[  I]+G*G;next  I 

96:  for  1*30  to  226  ;R(I]  +  F*F 

97:  if  F>*G/2;gto  99 

98:  next  I 

99:  I* A 

100:  if  0*1  and  X*2;l*rl 

101:  rl-A*r2 

102:  for  J*1  to  256+r2 

103:  R[J-c2] *R[  J  J 

104:  if  J<256+r2; nex t  J 

105 :  for  K*25o+r2  to  256 

106:  0*R[KJ 

107:  next  K 

108:  sprt  7,r2 

109:  rread  4,l;sread  4,1 [ * 3 

110:  for  1*1  to  256 

111:  I  [I]+R(I]  *R(  I]  ;next  I 

112:  rread  4,1?  sprt  4,R[*] 

113:  next  X;next  0 

114:  "TAKE  AVERAGE  OF  ALL  FRAMES "  : 

115:  for  1*1  to  256; R[ I ]/ (H-l) *R( I ]  ;R[  I  513  -I  ]  ;  nex  t  I 
116:  rread  4,1; sprt  4,R[*] 

117:  if  0 [15] <1; gsb  "OUTPUT" 

118:  beep 

119:  ent  "RANGE", R,"SQCN",Q, "EXT  CCEF"  ,E,  "laAVLNGTfi"  ,  h,  "OBJ  LE 
NS"  ,0 

120:  ent  "OBS/OBJ",  B, "SCALE" , rll ; Sprt  6  ,R,  Q,  E  ,  W  ,0  ,B,  rll 
121:  "SYSTEM  DATA-W-WAVLNTH  IN  METERS  0*DIA  OF  OBJ  LENS  IN  ME 
TERS" : 

122:  "B*DIA  OF  OBSUR/DIA  OF  OBJ  LENS  rll*SCALE  OF  DATA  IN  MIC 
RORADS/PT"  : 

123:  "R-DISTANCE  TO  TGT  IN  METERS": 

124:  "0-SQUARE  OF  ATMOS  TURB  INDEX  IN  METERS~-2/3" : 

125:  " E* EXTINCTION  COEFF  (1/METERS)": 

126:  if  O [1 ] #5 ;gto  140 
127:  beep 

128:  ent  "AMPLITUDE", A, "SIGMA" ,C 

129:  for  1*1  to  256 

130:  rll*(I-l)*F 

131:  F‘>2/(2*C“2) -*3 

132:  if  G>13 ; gto  134 

133:  A*exp  (-G)  *I.[I] ;  I*K;  gto  135 

134:  0*1(1] 

135:  next  I 

136:  rread  5,l;sprt  5,I(*];rread  5,l;sread  5,R(*j 
137:  1*0(15]  ;if  0.[14]#G?gsb  "OUTPUT" 

138:  gto  149 

139:  "SOURCE  PATTERN  DATA  INPUT": 

140:  rread  2,l;sread  2,R[*] 

141:  for  1*1  to  256 
142:  I *K 
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143:  if  R[ I] <-500; gto  146 

144:  next  I 

145:  gto  149 

146:  for  J=I  to  256 

147:  0*R [JJ 

14  8:  next  J 

149:  0[1)-t 

150:  "BRANCH  ACCORDING  TO  TYPE  OF  SOURCE  PATTERN  DATA": 

151:  gto  152;if  I>l;gto  155;i£  I>2;gto  161; if  I>3;gto  206;if 
I >4; gto  152 

152:  "CONCERTS  SOURCE  POINT  SPREAD  FUNCTION  TO  LINE  SPREAD  FC 
N"  : 

153:  gsb  "L3F" 

154:  2*0 [15]; if  0[3]>0;gsD  "OUTPUT" 

155:  "CALCULATION  OF  FOURIER  .TRANSFORM  OF  LINE  SPREAD  FCN  OF 
SOURCE"  : 

156:  for  1*1  to  256 
157:  R[IJ*R[513-I] ;nex t  I 
153:  gsb  "FXFORM" 

159:  rread  5,l;sprt  5,R(*J 

160:  3*0  [  15  ]  ;  if  0[4]>0;gsb  "OUTPUT" 

161:  if  0 [  2]  *2;  gto  173 

162:  "OPTICS  DIFFRACTION  INPUT  IF  NOT  CALCULATED" : 

163:  rread  8,l;sread  8,R[*] 

164:  for  1*1  to  256 
165:  I*K 

166:  if  R ( I ] <-500 ; gto  169 

167:  next  I 

168:  gto  182 

169:  for  J*I  to  256 

170:  G*R[J];next  J 

171:  gto  182 

172:  "CALCULATION  OF  DIFFRACTION  LIMIT  POINT  SPREAD  FCN": 

173:  B~2*D; 1-D*H; 3.14159e-6* rll*0/W*Z 

174:  for  1*1  to  256 

175:  Z*  (1-1) *Y jif  Y >30 ;gto  178 

176:  (  'AIRY  '(  Y)-D* 'AIRY' <  Y*B)  )  *2/d “2*R[  I] 

177:  I*K ; gto  179 
178:  0*R[ I ] 

179:  next  I 

130:  rread  8,l;sprt  8,R[*J 

161:  4*0  [15]  ; if  0.[6]>0;gsb  "OUTPUT" 

182:  if  O[2]*l;gto  138 

183:  "CONVERTS  OPTICS  POINT  SPREAD  FCN  TO  LINE  SPREAD  FCN": 
184:  rread  8,l;sread  8,R[*] 

135:  gsb  "L3F" 

136:  rread  6,l;sprt  3,R[*] 

167:  5*0 [ 15] ; if  0[7]>0;gso  "OUTPUT" 

188:  "TAKES  FOURIER  TRANSFORM  TO  GET  OPTICS  MODULATION  TRANS F 
£R  FCN": 


58 


± 


lo9 

190 

191 

192 
193: 

194 

195 

196 

197 

198 

199 
200; 
201: 
202 
203: 
20  4 

205 

206 

207 

208 
209: 

210 

211 

212 

213 

214 

215 
21o 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 
227 
223 
229 


rread  8,l;sread  3,R[*] 
for  1*1  to  256 
R[I J*R[513-IJ ;next  I 
gsb  "FXFGRM" 
for  1*1  to  256 
if  R[l]=0;gto  196 
R [ I J/R [ 1 ] X  J ;next  I 
rread  8,l;sprt  o,R[*] 

6*0[15];if  0  [  3  j  >0; gsb  "OUTPUT" 

"CALCULATION  OP  THE  PRODUCT  OF  TWO  FOURIER  TRANSFORMS" : 

rread  8,l;sread  8,R[*];rread  5,l;sread  5,L(*] 

for  1*1  to  256 

R[I]*I [I]*R[IJ  ; nex t  I 

rread  10,1; sprt  10, R[*] 

7*0 [15]; if  0(5] >0; gsb  "OUTPUT" 
exp( -E*R) *X 

"ATMOS  MODULATION  TRANSFER  FCN  CALC-SHORT  TERM"; 

0*R*21.  6*A;  W*  (-.  333333 3.) *B ;  C*V;  0*L;X*R[1] 
for  1*2  to  256 

gto  214;  if  L<*0;gto  215 ;  if  L*0;gto  20  9 

976. 5625 *(  1-1 ) /r 11*T ; A*T*1 . 666 667 *(B- (T/O)  ~  .3  33333  3) *0 

sfg  14 

if  G>13;1*L 

if  G-V < 0 ;-l *L 

G*V; X*exp(-G) *R[ I ] ; I*K; cf g  14;gto  218 
0*R[ I ] ; gto  218 

R [K ] *cos ( ( I-K) *  3 . 1415 9/K) *  2*R [ I ] 
gto  21 7; if  I-3*K/2<0;gto  213 
I*K ;1*L 
next  I 

rread  9,l;sprt  9,R[*] 

8*0  [15];  if  O  [  9]  >0;  gsb  "OUTPUT" 

"TRANSFER  FCN  OF  SOURCE*OPT ICS*ATMOSPHERE" : 
rread  10,l;sread  10, R [*]? rread  9,l;sread  9,1  [ * ] 
for  1*1  to  256 
I [ I ] *R[ I ] *1 [ I J 
I.  [  I  ]  *R  [  I  ] 

R[  I] *R[513-I] 
next  I 

9*0[  IS]  ;  if  O  [  10]  >0;  gsb  "OUTPUT" 

"INVERSE  FOURIER  TRANSFORM  GIVES  TARGET  LINE  SPREAD  FCN" 


230:  sfg  l;gsb  "FXFORM"; I*K;cf g  1 
231:  10  *0[  15]  ;  if  Q[ll]>0;gsb  "OUTPUT" 

232:  "CONVERTS  LINE  SPREAD  FCN  TO  POINT  SPREAD  FCN  BY  ABEL  TR 
AWSFORM" : 

233:  gsb  "ABEL" 

234:  11*0[ 15] ; if  Q[12]>G;gsb  "OUTPUT" 

235:  "CALCULATES  THE  FRACTION  OF  THE  POWER  INSIDE  A  CIRCLE  Of 
RADIUS  R" : 
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236:  gsb  "B" 

237:  12*0(15]; if  Q(13]>0;gsb  "OUTPUT- 

238:  "READ  IN  NEW  RANGE  AND  Q  IF  RANGE  NEG.  READ  IN  ALL  NEW  D 
AT  A"  : 

239: "“STARTING  WITH  I/O  CONTROL  CARD, IF  RANGE=C  PROGRAM  STOPS 

«  . 


240:  "IF  RANGE  IS  POS. 
£  NEW"  : 


241 


"RANGE  AND  Q" : 


THE  PROGRAM  CARRIES  OUT  CALC'S  WITH  .TH 


242 

243 
24  4 
245 
24  6 

247 

248 

249 

250 

251 

252 

253 

254 


if  O[l]*0;gto  245 

ent  "RANGE", R, "SgCN" ,U, "EXT  COEF",E 
gto  206; if  R<=0;gto  119; if  R=Q ;gto  245 
end 

"ABEL" :R(l]*N;1.4R[l]-1.8R[2]+.4R[3]*R[lj;£or  1=2  to  K 
R  ( I J  ♦M 

. 4N+.2M-.6RI 1+1] *R[I] ;M*N;next  I 
for  1=1  to  K 

R(  I]  /  (2*V (  (I  +  .  1)  *2-11)  )  *R(  I] 
for  J*I+1  to  K 

R(  I]  +R[  J]  /V  (  (J  +  .  1)  *  2-1 1 )  *R[  I] 
next  J;R(  I]/ir*R[I  ]  ;dsp  I;next  I 
ret 


255 

256 

257 

258 

259 

260 
261 
262 

263 

264 


"LSF"  : for  1*1  to  K 
l*J;dsp  T;R(I]*Z 
/(I-I+JM)  *Y 

2* ( (1-f rc(Y) }*R{ int (Y)  ]  +  frc ( Y)  *R[int  ( Y)  +1] )  +z*z 

J+l*J;if  Y<K ; jrep  -2 

Z*R(lj;next  I; ret 

"B"  :.25irR(l]*R(l] 

for  1=2  to  K 

2irIR[I]+R[I-l]*R(I] 

next  I 


265 

266 

267 

268 

269 

270 

271 

272 
27  3 
27  4 

275 

276 

277 

278 

279 

280 

281 

282 

283 


if  R[K]*0;ret 

for  1*1  to  K;R[ I.]/R[ K]  *R[ I ] ; next  I;ret 
"AIRY" :if  pl<0; beep;asp  "ERROR  pl<0";stp 
if  pl*0  ;  l*r4;  ret  r4 
0*r5;if  pl>15;jmp  2 
20+lQpl-pl~(2/3)*r6;jmp  2 
90+pl/2*r6 

if  pl<5 ; 6+pl*rl2 ; jrap  2 

1. 4*pl+6Q/pl*rl2 

raax(  int(  rl2 ) ,  int  { 3-t-pl/4 )  )  *rl2 

for  M*rl2  to  r6  by  3 ; le-28*r8 ; Q*r 13*rl4 

sfg  10;if  M/2*int(M/2)  ;cfg  10 

for  J*1  to  M-2;  2*  (M-J)  *r8/pl-rl3*r!5  ;r  8*rl3 

rl5*r8;if  M-J-2*0 ; rl5*r 4 

cmf  10  ;  rl4+2*r8  *f  Igl0*rl4  ;  next  J 

2*r8/pl-rl3  *rl5 

rl4+rl5*rl4 ; r4/rl4*r4 

if  abs(  r4-r 5) -abs  (r 4*le-6)  < *0  ; 2*r4/pl*r 4 ;  re t  r4 
r4*r5 ;nex t  M 
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284:  beep;dsp  "Accuracy  not  obtained" ;wait  1000;ret  r4 

285:  "FXFORM" :rad;9-N;ina  I 

286:  h/2~(N-1)+T 

287:  for  M=1  to  N ; 2* (N-M) ♦rl6 

288  :  for  J*0  to  2*(M-l)-l;cll  'BI '( J,P,N-1) 

289  :  cos  (PT)  *C;sin(PT)  (l-2flgl)*P 
290:  for  I=2rloJ  +  l  to  2rl6J+rl6 
291:  R[I]*rO;R[I+rl6]»rl 

292:  I[I]*r2;I(I+rl6]*r3 

293:  rO+rl*C+r3*P*R[  I]  ;  r2+r3 *C-rl*P*I  [I] 

294:  rO-rl*C-r3*P»R[I+rl6] 

295:  r2-r3  *C+rl*P*  1 1 1+rl6] 

296:  next  I;next  J;next  M 

297:  for  1=0  to  2  N-l;cll  BI  *  ( I ,  J,  N)  ;  if  I-J>0;jmp  8 

298:  if  I*J ; jmp  5 

299:  R[I+1]/V(2-N)-P 

300:  I.[I  +  l]/y  (2“N)*Z 

301:  R[J+1]-R[I+1]  ;I(J+1]-L[I  +  1] 

302:  P-R[  J+l]  ;Z*I  [J+l] 

303:  R[I+1]/V(2~N) *R[I+1] 

304:  I[I+1]//(2‘N)*I[I+1] 

305:  next  I;deg;ret 
306:  "BI" :Q-*p2; pl»p4 
307:  for  Z*1  to  p3 
308:  p4/2-*p4;  2p2*p2 
309:  if  frc(p4) #0; p2+l*p2 
310:  int(p4)  ♦p4;next  Z;ret 
311:  "SOURCE" :sfg  2 

312:  beep;dsp  "DISPLAY  SOURCE  IMAGE  ON  VID  DISC"; wait  3000 
313:  Cll  'SCALE  ' 

314:  r read  2,l;sprt  2,R[*] 

315:  cfg  2; ret 

316:  "SCALE" :if  flg2;gto  318 

317:  beep;dsp  "DISPLAY  SCALE  IMAGE  ON  VID  DISC"; wait  3000 
318:  dsp  "THEN  PRESS  CONTINUE" ; stp 

319:  ent  "WINDOW  WIDTH,  FIRST  PI XEL " , r 9 , "LAST  PIXEL",rl0 
320 :  for  1=2  to  3 

321:  char(47)&char  (55)  &char  (64)  schar  (8(l)&char  (96)*C$ 

322:  if  1=2; "930  Q2Q"*D$ 

323:  if  1*3;  "020  "*D$ 

324:  if  len(D$)  <4  ;  jmp  5 

325:  char  (1)&"  2"&char  (2)  &C$&char  (15)  &D$[  5 , 7]*B$ 

326:  cli  7 
327 :  wait  100 
323:  wrt  704, B$ 

329:  char.  ( 46 ) *C $  [1, 1J 

330:  char  (L)&"  2  "&char.(2)  &C$&char.  <15)  s.D8{l,.3)-AS 
331:  cli  7 
332:  wait  100 
333:  wrt  704, A$ 
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334:  next  I 

335:  "1  "&"0"&"A"&"P"&char.(96)*P$ 

336:  1*J ; 0** 

337:  for  L=1  to  256 
338:  asp  L 

339:  0*D;  ""*0?;*+l** 

340:  **512+J*S 
341:  for  1*2  to  6 

342  :  int(S/16) *T;T*U;S-T*16+I*lo*T 
343:  U*S;0?&cnar  (.T)  *0? 

344:  next  I 

345:  cftar  (1) &"*"&char  (34)  &char  (2)  &0?&char  (  3)  &P?&char  (15)*B? 

346:  cli  7 

347:  wrt  704.1, B$ 

348:  buf  "DATA" 

349:  tfr  704 , "DATA" , 256 

350  :  rds  ("DATA") *E;if  E*-l;jmp  0 

351:  for  I*r9  to  rlO 

352:  num(  I  $  [I  ,1  ]  )  ♦Z ;  D+Z*D 

.353:  next  I 

354:  D/ (rl0-r9) *R(L] 

355:  next  L 

356:  if  f lg2; gto  358 

357:  if  0(16]  *2;  gsb  "OUTPUT" 

358:  ret 
359:  "OUTPUT": 

360:  beep 

361:  dsp  "STOP  IF  PLOTTER  NOT  READY  THEN  CONTI NUE " ;wai t  5000 

362:  sfg  14 'MIN '  (  256)  /  'MAX '  (  256)  *rl7  ;cf  g  14 

363:  if  rl7<-.l;scl  - . 75 , 1. 05 ,-30 , 264 ;plt  1,0,1; jmp  2 

364:  sc 1  -.15,1.05,-30,264; pit  1,0,1 

365:  for  1*10  to  0  by  -1 

366:  pit  1/10,0, 2;plt  1/10,256/150,2 

367:  pit  1/10  ,0,2; nex  t  I 

.368:  for  1*0  to  16 

369:  pit  0 , 161 , 2;  pit  1/150,161,2 

370:  pit  0,161  ,2 

371:  next  I 

372:  pen;csiz  .9, 1,1. 4 

373:  fxd  0;for  1*16  to  0  by  -l;plt  -.1,161-1,1 
374:  lbl  16(16-1) ;next  I 

375:  fxd  l;for  1*0  to  10;plt  I/10-. 06 ,-5 , 1 
376:  lbl  1/10 ; nex  t  I 

377:  if  0(15] *10  or  O [ 15 ]*0; "TARGET  LINE  NUMBER"**? 

378:  if  O  (15]  *1  or  0(15]  *11;  "  POSITION  VALUE"**? 

379:  if  0(15J *10; "PREDICTED  NORMALIZED  LINE  VALUE" *X? 

380:  if  O [ 15 ]-0;" MEASURED  NORMALIZED  LINE  VALUE"*X? 

381:  pit  -.12,51, l;csiz  1. 2 ,1, 1. 4 , 90 ; lbl  *? 

382:  pit  0 , -13 , l;csiz  1. 2 , 1, 1. 4 ; lbl  X? 

363:  csiz  .5,1, .7,0 
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384:  'MAX'(25b) *tl 

385:  if  r  7*0 ;  l*r7 

386:  for  J*1  to  256 

387:  pit  R(  J]/r7,257-J;next  J 

388:  pen;ret 

389  :  "MIN"  :A(l]-p2 

390:  for  1*2  to  pi; if  R[I]>p2;jmp  2 

391:  R[I]*p2 

392:  next  I; ret  p2 

393:  "MAX":Rtl]*p2 

394:  for  1*2  to  pi; if  R[I]<p2;jmp  2 

395:  R [ I J ♦p2 

396:  next  I  ;ret  p2 

*21194 
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ABSTRACT 


This  appendix  presents  a  comparative  analysis  of  common  path  (laboratory 
environment)  measurements  of  Cn2  for  a  telescope/slit  scan  system  and  a 
video  system.  Included  is  a  computer  program  in  algorithmic  language  that 
uses  a  digital  storage  oscilloscope  to  analyze  the  video  data.  Both  HeNe  and 
GaAs  lasers  are  used  in  the  comparative  measurements. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  increasing  use  of  lasers  and  laser  technology  for 
military  applications  has  brought  about  a  need  for  analysis 
of  the  laser  beam  in  its  environment,  the  turbulent 
atmosphere.  A  project  at  Naval  Postgraduate  School  dealing 
with  this  subject  is  continuing  and  is  the  main  topic  of 
this  thesis. 

The  patterns  produced  by  lasers  on  targets  have  inherent 
problems  that  include  broadening,  beam  wander  and  intensity 
fluctuations  brought  about  by  turbulence  in  the  atmosphere. 
These  effects  of  atmospheric  turbulence  or.  laser  propagation 
have  been  well  determined  (Ref.  1].  In  terms  of  the  Fried 
model,  [Ref.  2]f  C  n2,  the  refractive  index  structure 
constant,  has  been  determine!  to  adequately  express  the 
Modulation  Transfer  Function  (MTF)  ,  cr  the  Mutual  Coherence 
Function  (MCF)  for  the  atmosphere. 

A  system  has  been  developed  that  provides  a  measurement 
of  C  „  2  for  atmospheric  turbulence  along  the  optical  path 
through  which  a  laser  is  propa;ated.  The  system  employs  a 


vidicon  and  telescops  as  the  dstsctor  and  a  distant  laser  as 
the  source.  This  system  duplicates  the  slit  scanning  system 
presently  in  use  at  Naval  Postgraduate  School  [Ref.  2]. 
Tests  using  the  vidicon  equipment  have  been  previously 
completed  with  values  for  Cn2  on  the  order  of  10~15  m-2/3 
being  obtained  [Ref.  3].  Measurements  have  been  made  using 
the  same  experimental  set-up,  except  that  the  Tektronix  468 
Digital  Storage  Oscilloscope  are  used  for  digitization 
instead  of  the  Quantex  DS-30  Digital  Video  Processor.  These 
measurements  have  demonstrated  values  of  Zn2  of  comparable 
ac curacy. 

B.  PROBLEM 

By  using  the  digital  storage  capability  of  the  Tektronix 
468,  the  previously  recorded  laser  signal  is  used  as  an 
input  and  is  evaluated  by  modifying  the  program  developed  by 


Crager  [  Ref  .  4  ]. 

k  brief  ovarvi 

ew 

of  the  procedure 

is 

described  below. 

Detailed  explanti 

ons 

of  the  experiment 

and 

computer  program 

are  contained  i 

n 

Chapters  III  and 

IV, 

respectively. 

The  approach 

taken  is  basical 

iy 

the  same  as  that 

of 

Crager,  but  because  of  the  unavailability  of  a  disk  RDM  it 
is  necessary  to  store  both  raw  and  processed  data  cn 
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magnetic  taps.  Its  data  are  leaded  and  recorded  as 
necessary  by  the  HP  93  25.  The  basic  assumption  permit  ting 
analysis  is  that  the  horizontal  TV  scan  line  through  the 
laser  spot  is  considered  to  accurately  mirror  a  point  spread 
function  of  the  iaage.  Pinal  analysis  using  the  463 
oscilloscope  has  shown  this  to  be  a  valid  hypothesis.  The 
digitized  data  from  the  463  oscilloscope  agree  with  the 
previously  measured  lata  which  use  many  pixels  of  television 
data  digitized  by  the  DS-30  to  express  the  point  spread 
function . 

The  sequence  of  analysis  is  that  a  video  tape  recording 
is  made  of  the  TV  image  of  the  propagated  laser  beam.  The 
output  of  the  video  recorder  is  sent  to  the  Tektronix  468, 
where  the  derived  TV  scan  line  is  digitized,  averaged,  and 
stored.  The  HP  9825  records  the  digitized  data  and  produces 
a  line  spread  function  (LSF)  ,  by  integrating  the  point 
spread  function, 

y  (max) 

L3  F(x|  =  J PSF  (r)  iy  (1-1) 

y  (min) 

where 


1 1 


r  =  ^ixz+Y2)  ' 

and  computes  the  Fourier  transform  of  the  LSF.  The 
diffraction  limited  Fourier  transform  of  the  optics  is  n oh 
computed  if  the  MIF  of  the  optics  has  not  been  previously 
measured  experimentally.  Next,  the  program  finds  the  STF  of 
the  atmosphere  by  dividing  the  Fourier  transform  of  the  LSF 
by  the  Fourier  transform  of  the  optics.  Finally,  by  curve 
fitting,  the  program  computes  a  single  value  for  Cn2. 

The  program  now  predicts  the  size  of  a  laser  spot  on  a 
target  using  the  calculated  value  of  Cn2.  This  value  of  cn2 
is  used  to  calculate  an  ilTF  of  the  atmosphere  which  is  then 
multiplied  by  the  Fourier  transform  of  a  source,  and  the 
Fourier  transform  of  the  optics.  The  program  then 
calculates  the  inverse  Fourier  transform  of  the  products  of 
the  above  and  uses  the  Abel  transform  to  give  the  angular 
point  spread  intensity  distrioution.  From  this  data,  the 
fraction  of  energy  as  a  function  of  the  total  energy  within 
a  given  radius  a  is  calculated. 
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II.  THEORETICAL  considerations 

The  theory  of  laser  beam  propagation  through  a  turbulent 
medium  has  been  explained  by  Crittenden,  and  others,  and  is 
re-emphasized  here  for  continuity  purposes  [Ref.  2].  Since 
the  effects  of  turbulence  on  laser  weapons  is  of  major 
concern,  measurement  and  prediction  capabilicy  for  these 
effects  or.  laser  beams  is  entirely  relevant. 

Due  to  the  existence  of  reciprocity,  these  ideas  apply 

to  either  laser  designators  or  imagers  [Ref.  5].  In  the 

study  of  the  theoretical  model  by  Fried,  the  effects  of 

atmospheric  turbulence  are  ii  vestigated  [Ref.  6].  This 
0 

model  uses  the  idea  of  a  long  term  optical  transfer  function 
(OTF)  when  considering  atmospheric  turbulence.  The  long 
term  OTF  results  from  simply  ta<ir.g  an  image  of  sufficiently 
long  term  which  sees  effectively  all  possible  turbulence 
configurations. 

When  observing  tae  effect  of  diffraction  as  the  result 
of  the  finite  aperture  of  a  point  source,  it  car.  be  seen 
that  the  resulting  image  is  not  uniform.  Considering  the 
point  spread  functioa  to  have  the  same  shape  in  the  image 
plane  regardless  of  its  position,  the  image  function  is  the 
result  of  the  convolution  of  the  source  function  and  the 
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optics  diffraction  function.  The  convolution  theorem,  as 
described  in  Fourier  transfers  theory,  yields 

i  (v  (x)  ,  v(y))  *H  (v(x)  ,  v  (y)  |  *o  (v(x)  ,  v(y) )  (v  (x)  ,  v  (y) )  (2.1) 

where 

i(v(x),v(y))  =  Fourier  transform  of  the  image  function 

H  (v  (x)  ,v  (y)  )=  Fourier  transform  of  the  optics  diffraction 
function 

o(v(x),v(y))  *  Fourier  transfers  of  the  object  function 

B  (v  (x)  ,  (v  (y ))  =  Modulation  Transfer  Function  of 

the  atmosphere 

v  =  spatial  frequency 

A  point  sourca  such  as  a  laser  can  be  analyzed  in  two 
dimensions  using  a  point  spreai  function.  Through  Fourier 
transform  theory,  the  image  point  spread  function  is 
transformed  into  a  two  dimensional  OTF  of  the  optical 
signal.  This  problem  may  be  simplified  by  scanning  the 
image  point  spread  function  using  a  vidicoa  or  slit-scanning 
system.  In  applying  the  convolution  theorem,  the  Fourier 
transform  of  the  image  point  spread  function  is  multiplied 
with  the  Fourier  transform  of  the  optical  system  resulting 
in  the  Fourier  transform  of  the  overall  system.  The  Abel 
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•transfer a  described  by  Grisa,  is  applied  to  this  result  to 
re-transf or  a  the  oi a -dimension al  image  LSF  into  a  two 
dimensional  image  PSF  [Ref.  7]. 

As  demonstrated  by  Crittenden,  and  others,  a  numerical 
value  for  CQ2  may  be  obtained  by  curve  fitting  using  the 
following: 

M  =  exp(-21.49*Cn2*Z*f5/3*\~1/3)  (2.2) 

where 

M  =  MT?  of  the  atmosphere 
2  =  Range  in  meters 

f  s  f*v  =  Angular  spatial  frequency  in  cycles/radian 
F  =  focal  length  of  the  optical  system 
v  =  linear  spatial  frequency  in  cycles/meter 
A  »  wavelength  in  aecers 

Cn2  is  obtained  by  a  linear  regression  of 

In  (t!)  =  -21.49*Cn2*Z*fs/3*;ri /3  (2.3) 

where  Cn*  is  the  only  parameter. 
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•  EXPERIMENTAL  PROCEDURE 

The  experiment  is  performs!  using  two  different  laser 
sources,  helium-neon  (He- Ne)  and  gallium-arsenide  (3a-As). 
Figure  3.1  shows  a  block  diagram  of  the  experimental  set-up 
used  for  both  lasers. 

The  measurement  of  Cn2  along  the  optical  path  is  made  by 
using  a  viaicon  and  telescope  at  the  far  end  cf  the  corridor 
in  the  basement  of  Spanagal  Hall.  Atmospheric  turbulence  is 
produced  by  nine  overhead  hot  air  ventilators.  The  optical 
equipment  similar  to  that  described  in  reference  3,  consists 
of  a  6  inch  diameter  Cassegrain  telescope  with  a  90  inch 
focal  length  [  Ref .  3]. 

Using  neutral  density  absorption  filters  to  attenuate 
the  intensity  of  the  laser  beam,  the  telescope  is 
illuminated  at  its  input  aperture.  The  light  image  from  the 
telescope  is  split  by  a  beam  splitter  with  one  beam  sent  to 
the  vidicon  and  the  other  to  a  slit  scanner  for  comparison 
of  the  two  systems.  The  light  is  then  transformed  to  an 
analog  signal  and  recorded  on  a  Panasonic  (NV-9240)  video 
tape  recorder  for  processing.  The  data  from  the  slit 
scanning  system  are  recorded  on  a  precision  instrument  tape 
recorder  for  analysis.  After  recording  the  video  signal. 
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Fig.  3.  1  Laser  Experimental  Set-Up 


analysis  is  accomplished  using  a  Tektronix  468  Digital 
Storage  Oscilloscope  and  a  Hewlett  Packard  9825  Computing 
System.  The  linearity  of  the  Panasonic  tape  recorder  was 
demonstrated  in  Cragec's  thesis.  This  is  also  confirmed  by 
observing  the  images  as  real  time  and  recorded 
displays  [Bef.  3], 

The  468  oscilloscope  is  set-up  in  the  following  manner. 
An  input  signal  to  channel  A  is  located  on  the  scope  using 
the  non-storage  mode.  The  signal  is  then  displayed  using 
the  A  INTENSITY  switch  of  the  horizontal  display.  This 
intensified  zone  is  used  to  position  the  3  sweep  (delayed} 
to  the  desired  location  within  the  A  sweep  interval  to 
obtain  an  expanded  view  of  a  waveform  for  examination.  Once 
the  waveform  is  centered  on  the  scope,  the  horizontal 
display  is  switched  to  S  DLY’D.  This  is  done  to  facilitate 
the  digital  storage  circuitry  time  base  by  using  the  setting 
of  the  3  TI  ME/D  17  switch. 

The  waveform  is  expanded  la  time  by  decreasing  the  A 
TINE/DIV  switch  setting  and  moving  the  waveform  back  to  the 
middle  of  the  scope  using  the  delay  time  position  control 
dial.  When  a  representative  waveform  is  obtained,  the  B 
TIME/DI7  switch  is  then  used  in  conjunction  with  the  delay 
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time  position  control  dial  to  expand  the  tin-  scolo-  until 
the  single  central  maximum  of  the  waveform  is  centered  and 
one  horizontal  sweep  is  displayed. 

Once  this  waveform  is  satisfactorily  obtained,  the  AtfG 
storage  mode  is  selected.  The  468  will  average  the  input 
signal  for  a  selected  number  of  sweeps  and  display  the 
accumulated  wavefora.  All  data  for  this  thesis  use  32 
sweeps  for  each  average.  The  468  is  now  ready  tc  transfer 
data  when  interrogated  by  the  8?  9825. 

The  computer  is  the  controller  for  all  interfacing 
operations  with  only  a  aiaimal  amount  of  operator 
interaction.  The  operator  interface  is  mainly  to  ensure 
that  the  equipment  is  properly  set  up  and  to  select  if  data 
are  to  be  plotted.  Digitization  of  the  analog  signal  is 
accomplished  by  the  Tektronix  45^.  Processed  waveform  data 
are  transferred  fron  the  microprocessor  memory  to  the 
Storage  Display  RAM  [Ref.  8]. 

The  controlling  program  of  the  HP  9825  interrogates  the 
463  via  the  IEEE  488  interface  bus.  The  468  receives  the 
data  request  from  the  controller  and  sends  the  waveform 
messaqe,  both  preamble  and  data.  The  waveform  message  is 
stored  in  the  calculator  memory  for  further  processing. 
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When  the  message  is  completed,  the  468  concludes  wita  an  end 
of  instruction  terminator  and  the  controller  takes  control 
of  the  bus  again. 

Data  processing  of  the  digitized  signal  begins  with  the 
data  being  stored  on  magnetic  tape  for  further  use  by  the 
computer.  The  data  are  stored,  processed  and  plotted  by  the 
main  program.  Subroutines  are  called  as  necessary  for  their 
specific  uses.  The  two  signal  waveforms  are  recorded  and 
processed.  First,  the  signal  from  a  laser  beam  incident  at 
the  aperture  of  the  telescope  is  recorded.  Figures  3.2  and 
3.3  show  these  recorded  data  wiich  are  referred  to  as  the 
source  data.  Second,  for  calibration,  the  signal  from  a 
laser  beam  with  a  diffraction  grating  in  place  at  the 
aperture  of  the  telescope  is  recorded.  The  entire  system  is 
calibrated  by  using  a  grating  that  consists  of  closely 
spaced  vertical  bars  in  front  of  the  telescope.  This 
produces  a  diffraction  pattern  in  the  image  plane  and  allows 
for  the  calculation  of  the  scale  factor.  These  recorded 
date  are  referred  to  as  the  scale  data  and  are  plotted  in 
Figures  3.4  and  3.5. 

During  the  recording  process,  the  signal  from  the  laser 
beam  is  recorded  on  video  tape  while  the  signal  is  observed 
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Fig.  3.2  Plot  of  Source  Doto  for  He-Ne  Laser 
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Fig.  3.3  Plot  of  Source  Data  for  Ga-As  Laser 


q  Data  for  He~Ne  Laser 
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Fig-  3.5  Plot  of  Scale  Data  for  Ga-As  Layer 


o?.  “he  45  8  cscilicsc  o  pe.  This  procedars  is  important  ir. 

da t a  t  aK  u.n  g  since  xiss  r  axrgnment  end  texescope  z  oc  u  s  t.  n  g  arc 

often  very  sensitive  to  minor  movements  of  either.  a pon 
completion  of  recoriing  the  data  on  video  rape,  the 
experiment  is  concerned  with  the  compilation  and  analysis  of 
data.  The  first  item  to  be  determined  is  the  scale  factor. 
Once  the  scale  factor  has  bean  determined,  it  will  remain 
constant  throughout  the  calculations,  unless  the  telescope 
focal  length  is  changed. 

The  scale  factor  is  calculated  ir.  the  following  manner. 
The  number  of  points  between  the  peaks  in  the  plot  of  the 

scale  data  (diffraction  grating  in  place  in  front  of  the 

telescope  aperture)  is  measured.  When  the  distance  between 
peaks,  the  spacing  between  the  bars  in  the  grating,  and  the 
wavelength  cf  the  laser  are  known,  the  scale  factor  can  be 
calculated  from  the  relation 

sin  (theta)  =  A  /d  (3.1) 

sin  (theta)  ~  theta  (3.2) 

(small  angle  approximation) 

sf  =  (theta)Zndp  (3.3) 

sf  =  A / (d*nd?)  ( 3. 4) 
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where 


sf  =  scale  factor  in  radians  per  point 
A  =  wavelength  in  meters 

d  =  spacing  between  lines  on  grating  in  meters 

ndp  =  number  of  points  between  central  maximum 
and  first  order  diffraction  peak 

This  standardizes  the  data  for  the  abscissa  of  the  plots  ir. 
radians  par  point. 

The  program  now  takes  the  transferred  data  and  computes 
a  point  spread  function  as  shown  in  Figures  3.5  and  3.7. 
Figures  3.3  and  3.9  represent  the  point  spread  function 
after  integration  using  equation  (1.1)  to  obtain  a  line 
spread  function.  Next,  the  Fourier  transform  of  the  line 
spread  function  is  calculated.  These  curves  are  plotted  in 
Figures  3.10  and  3.11. 

The  diffraction  of  the  optics  is  computed  and  plotted  as 
shown  in  Figures  3.12  and  3.13.  Following  the  same  method 
as  before,  the  line  spread  function  of  the  optics  function 
is  calculated  and  pLotted  in  Figures  3.14  and  3.15.  The 
program  now  takes  the  Fourier  transform  of  the  optics 
function  LSF.  These  data  are  plotted  in  Figures  3. 16  and 
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Fig.  3.  7  Plot  of  Point  Spread  Function  for  Ga-As  Laser 


3.17.  The  ?3u:isU  transforms  of  the  system  and  conics  are 
divided  point  by  point  to  yield  the  HTF  of  the  atmosphere, 
as  shown  in  Figures  3.18  and  3.19.  The  HTF  of  the 
atmosphere  will  be  used  to  calculate  the  long-term  value  of 

Cn2- 

Cn2  is  calculated  by  linear  regression  of  in  (.1TF)  versus 
f5/3.  The  slope  of  the  regression  is  proportional  to  Cn2. 
Two  simultaneous  equations  are  solved  using  Cramer’s  rule 
and  the  matrix  ROM  of  the  HP  9825.  The  equations  used  are 

a*X  +  b*I  =  Y  (3.5) 

a*X*  ♦  b*X  =  Y*X  (3.6) 

where 

a  =  slope  of  the  curve 
b  =  intercept 

X  =  summation  of  (I*sf )  **(5/3) 

X**2  =  summation  of  (I  *sf)  **  (1 3/3) 

I  =  point  number  (total  number  =  256) 

T  =  summation  of  natural  logarithm  of  MTF 
Y*X  *  summation  of  products  of  two  values 
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C n2  is  obtained  from  the  above  information  using  equation 
(2.3).  This  will  yield 

Cn2  =  a/(-21.  49*Z*f**  (5/3)  *  **(-1/3))  (3.7) 

The  program  now  goes  to  the  prediction  phase  after 
calculating  C  „  2  for  the  atmosphere.  If  a  Saussian 
distribution  for  an  input  source  and  a  value  of  the  standard 
deviation  are  assumed  or  known,  the  resulting  source 
function  can  be  calculated.  A  plot  of  the  computed  source 
is  shown  in  Figure  3.20.  It  can  be  seen  from  this  plot  that 
the  half- width  at  half-maximum  is  approximately  u 
microradians.  As  before,  the  line  spread  function  of  the 
computed  source  is  calculated  and  plotted.  This  is  shown  in 
Figure  3.2  1.  The  Fourier  transform  of  this  data  is 
calculated  and  plotted  in  Figure  3.22. 

The  program  now  multiplies  the  Fourier  transform  of  the 
computed  source  with  the  transform  of  the  system  (including 
the  atmosphere)  and  plots  the  result  in  Figures  3.23  and 
3.24.  Plots  of  the  inverse  Fourier  transform  are  shown  in 
Figures  3.25  and  3.25.  Next,  the  Abel  transform  is  computed 
and  plotted  in  Figures  3.27  and  3.23.  The  Abel  transform, 
as  described  in  Chapter  II,  transforms  a  one-dimensional 
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line  opr sa d  function  into  a  two-dimensional  poire  spread 
function.  Finally,  the  fraction  of  power  inside  a  circle  of 
radius  R  is  calculated  and  plotted  as  shown  in  Figures  3.29 
and  3.30.  This  is  the  fraction  of  power  that  one  would 
expect  to  be  incident  on  a  target  using  the  measured  value 
of  atmospheric  turbulence  as  an  input  parameter.  A  detailed 
analysis  of  the  computer  program  is  the  subject  of  the  next 
chapter. 
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Fig.  3.8  Plot  of  Line  Spread  Function  for  Hu-Ne  Laser 
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Plot  of  Fourier  Transform  of  LSF  for  Ga-As  Laser 
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Fig.  3.  14  Plot  of  LSF  of  Optics  Function  for  He-Ne  Lasar 
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Fourier  Transform  of  LSF  of  Optica  Function  for  He-Ne  Loser 
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Fig.  3.  17  Plot  of  Fourier  Transform  of  LSF  of  Optics  Function  for  Ga-As  Laser 


Fig.  3.20  Plot  of  Computed  Source 
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Fig.  3.21  Plot  of  Line  Spread  Function  of  Computed  Source 
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Fig.  3.  22  Plot  of  Fourier  Transform  of  LSf  of  Computed  Source 
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Fig.  3. 25  Plot  of  Inverse  Fourier  Transform  for  He-Ne  Laser 
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3.28  Plot  of  Inverse  Fourier  Transform  for  Go- As  Loser 


\bel  Transform  for  He-Nu  Laser 
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Fig.  3.28  Plot  of  Abel  Transform  for  Ga~As  Laser 
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:ig.  3.29  Plot  of  Power  Fraction  Inside  Circle  of  Radius  R  for  He-Ne  Laser 


IV 


OMFUrSH  PROGRAM  13  ALGORITHMIC  LANGUAGE 


The  program  presented  here  is  written  in  the  algorithmic 
language  as  outlined  by  Graham  [Ref.  9],  It  is  intended 
that  the  language  in  use  hers  be  macaine  independent. 
Comments  appear  in  parentheses. 


Algorithm_Laser 
(Input  which  laser  is  used.) 

Input  Laser_Wav eiength 
ESTRING  <==  Typejjf  _Laser 

(Input  if  the  particular  laser  source  data  are  recorded. 
The  following  lines  determine  if  the  laser  source  data  are 
already  recorded.  If  the  source  data  are  act  recorded  -.her. 
the  program  goes  to  the  transfer  subroutine  to  allow 
interfacing  with  the  468  to  effect  a  data  transfer.  The 
parameters  passed  to  the  Plot  subroutine  have  the  following 
significance:  the  first  parameter  determines  the  horizontal 

index,  for  example  512  points  or  the  number  of  microradi  ar.s. 
The  second  parameter  tells  the  computer  which  file  the  data 
are  recorded  on  so  tie  data  can  be  loaded  into  memory.  The 
final  parameter  determines  the  labeling  of  the  horizontal 


axis  either  in  data  points,  aicroradiars,  cr  li.-.^s  per 
microradian .) 

If  Laser_Source_no t_recoried  then 
Call  Transfer  (Source_Fiie) 

Call  Plot  (512  ,Source_Fila,  0) 

End  If 

(Input  if  the  particular  laser  scale  data  are  recorded. 
Scale  data  are  recorded  with  the  diffraction  grating  in 
place.  If  the  scale  data  need  to  be  recorded,  the 
subroutine  Transfer  is  called  and  the  scale  data  are  stored 
on  a  separate  file.) 

If  Lase r_Scale_aot ^recorded  then 
Call  Transfer  (Scale_File) 

Call  Plot  (512  ,Scale_fiLe,3) 

End  If 

Input  Scale_?actor 

(Calculation  of  the  scale  factor  is  detailed  in  Chapter  III. 
Comparison  of  the  time  base  of  the  recorded  data  is 
necessary  to  ensure  that  the  data  recorded  have  the  same 
time  scale  on  +he  oscilloscope.  If  the  time  scales  are 
different  the  program  stops  and  displays  in  error  message. 
A  file  is  loaded,  the  waveform  preamble  is  searched  for  the 
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"XONITS"  or  "XINCR1'  o  sell  loscope  settings  and  the 
particular  data  are  stored.  This  is  done  for  each  of  one  3 
recorded  files:  source,  scale,  and  collimated  data. 
Collimated  data  refer  to  the  same  laser  recorded  under 
controlled  conditions.  The  laser  beam  is  projected  in  a 
collimator  and  recorded  on  tape.  Th:  s  data  could  be  used  to 
remove  any  non-linear L ties  of  the  video  recorder,  electronic 
equipment,  or  receiving  optics.  Data  elements  are  compared 
to  ensure  consistency  of  inputs.) 

Input  Laser_source 

3<==  po  s(Laser_Sou  rce,  'R:  •) 

HI  <»=  val  ( »  B+2‘  ,3+4' ) 

3<==  po  s  (Laser.  So  a rce,  ' IT: ) 

AS  THIN  G  <==  Laser  _  Source  ('  3  +  3  *  ,  •  3+4  ’) 

Input  (Laser_Sca  lei 

B< -=  po s(Laser_Sca le, '  E  :  •  > 

R2<==  v  al ( ' B+2* ,3+4') 

B<==  po s(Laser_Sca le, ' IT: ') 

BSTR IN G <==  Laser_3 cale ( »8+3' , ' 3+4  •) 

Input (Colli mat ad_3  ource) 

B<==  po s (Collimate d_Source,'R:  ') 

R 3<==  val('  B+2'  ,'  B+4) 

B<==  pos(Colliaate  d_So  urce, '  IT:  ') 
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CSTRING<=  =  Ccliiaated_ Source ('B+3', »B+4') 

If  R 1  =  R  2=R  3  than 

If  ASTRINS  =  BSTRING  =  ESTRING  then 
Input  RDATA(File_tfumber) 

Bn d  If 

Else 

Output  'ERROR  IN  TIME  3 GALE' 

Stop 
End  If 

A<=-B<  =  =C<==0 

(Array  is  initialized  to  0) 

Do  for  I  =  1  to  51 2 

YD  ATA  (I)  <  =  =  R  DATA  (I) 

End  Do 

(This  section  takes  the  data  from  the  waveform  message  as 
transferred  from  the  468  and  removes  any  DC  background  from 
ar.  image  as  well  as  subtracting  the  minimum  value  in  order 
to  "zero"  the  data.  This  makes  y=0  as  the  starting  point 
for  both  calculation  and  plotting  purposes.  Data  up  to  51 
and  greater  than  463  are  summed  and  stored  in  B. 

This  represents  the  "wings"  of  the  curve.  This 
process  removes  the  DC  background  from  an  image.) 

Do  for  I  *  1  to  51  2 
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k<~  A  +  YDATA(I) 


If  I  <  51  than 

B  <==  YDATA(I)  ♦  B 
End  If 

If  I  >  463  then 

B<==  YDATA(I)  ♦  B 
End  If 
End  Do 

(Output  for  debugging) 

Output  A, 3 

(Average  of  the  sum  of  data  values  in  the  "wings".) 

D<==  B/100 

(This  is  the  total  sum  of  all  data  less  -he  "wings"  value. 
A<==  A -  51 2  *  D 

(All  data  are  summed  and  stored  in  C.  Each  data 
is  compared  with  the  previous  total  divided  by  2. 

This  will  find  the  index  value,  I,  for  which 
the  maximum  value  exists.) 

Do  for  I  =  1  to  51 2 

YD  ATA  (I)  <  =  =  r  DATA  (I)  -D 
C<  —  YD  AT  A  (II  *  Z 
If  C  <  A/2  then 
J<==  I 
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End  If 
End  Do 

(Output  for  debugging.) 

Output  C, J 

(The  curve  is  shifted  so  that  the  maximum  value  starts 
at  the  origin  of  the  horizontal  axis.  This  is  called 
the  point  spread  function  and  is  recorded  as  RDATA. ) 
Do  for  K  =  1  to  256 
IDATA  (K)  <=  = 

YDAT  A  (  ( J-K ♦  1)  mod51  2  +  1 )  +  YDATA  ( (K* J- 1)  mod5  124-1) 

End  Do 

Do  for  I  =  257  to  512 

♦ 

RDATA  (I)  <-=  3 
End  Do 

Do  for  I  =  1  to  256 

RD  ATA  (I)  <==  I  DATA  (I) 

End  Do 

(Determines  if  tabular  or  plotted  data  are  desired.) 

If  PSF_OOTPUT_waat  ed  then 
Call  Table  (RDATA) 

End  If 

If  PSF_ PL3T_wanted  than 

Ca  11  Plot  ( 32 ,  RDATA,  1) 
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End  If 

(Calculates  the  line  spread  function  from  the  point 
spread  function.) 

Call  LSF 

(Determines  if  tabular  or  plotted  data  are  desired.) 
If  LSF_OtJTPUT_want  ed  then 
Call  Table  (RDATA) 

Else  If  LS?_PLOr_a anted  then 
Ca  11  Plot  ( 32 ,  RDATA,  1) 

End  If 

(Calculation  of  Fourier  Transform  of  LSF.) 

Input  S  DAT  A 

(The  maximum  value  of  the  array  T D A T A  is  at  1=1.  Thi 
is  put  back  into  array  RDATA  and  folded  over  sc  that 
IDATA  (257)  to  IDATA(512)  is  the  mirror  image  of 
IDATA  (1)  to  IDATA(255  ).  Then  this  array  is  recorded 
as  RDATA  and  the  Fourier  transform  calculated.) 

Do  for  I  =  257  to  512 
RDATA  (I)  <  =  =  3 
End  Do 

Do  for  I  =  1  to  25  6 

RDATA  (514-1)  <==  RDATA  ( I ) 
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End  Do 


3D  AT  A  {  2  57)  <  ==  3041.4  (256) 

Call  FX  FORM 

(Determines  if  tabular  or  plotted  data  are  desired.) 

If  FXFORM_OF_LSF_OOTP3T_wanted  then 
Call  Table  (RD  AT  A) 

Else  If  FXFORM_OF_LSF_?LOr_tfanted  then 
Call  Plot  (32,  RDATA,2) 

End  If 

(Dia  of  Cbccur/Dia  of  Obj  Lens  in  meters, RATIO*. 064/.  164) 
(This  is  the  measurement  of  the  ratio  of 
diameters  of  obscuration  to  telescope.) 

( R 1  -  Scale  of  data  in  microradians=1.60) 

(Diameter  of  telescope  Obj  Lens  in  meters,  OBJECT* 

0.164  me  tars) 

If  Lase  r_S  o urce  =  He-N  e  then 
Wa velength<==  6.328*e-7 

Elsa 

Wa velength<==  9.35*e-7 
End  If 

(This  begins  the  calculation  of  the  Fourier  transform 
of  the  optics  using  the  Airy  function.  All  of  the 
data  regarding  the  optics  are  recorded  on  a  separate  file.) 
(These  next  four  statements  tic  the  constants  in  the 
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argument  of  the  optics  function.) 


3 1  <==  Scale_Factoc 
D<==  SATI0**2 
H<==  1- D 

Z<==  PI*Scale_Factor*Dbject't  1  *e-6/Wavelength 
Do  for  I  =  1  to  256 
Y  <  =  Z*(I-1| 

If  Y  >  30  than 

RDATA  (I)  <==  0 

Else 

RDATA  (I)  <==  (AI3Y{Y)  -D*AIRY  (Y*R.\TIO)  )  **2/4  **2 
End  If 
K  <  =  I 
End  Do 

(Record  RDATA  on  separate  fiLe) 

(Determines  if  a  plot  of  the  optics  function  is  desired.) 
If  OPTICS_FONCTIDJl_PLDT  wanted  then 
Call  Plot  (32,  RDATA,  1) 

End  If 

(Calculates  the  line  spread  function  of  the  optics 
function. ) 

Call  LSF 

(Determines  if  a  plot  of  the  L5F  of  the  optics  function  is 
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desired.) 


If  PLOT_OF_LSF_OPTICS  wanted  then 
Call  Plot  (32,  RDAI  A,  1) 

End  If 

Do  for  I  =  257  to  512 
RD  AT  A  (I)  <=  =  D 
End  Do 

(Dana  are  ro  —  ced  over  a <5  —  an  ana  t.i a  Fcurrer  transform 
calculated .) 

Do  for  I  =  1  to  256 

RDATA  (514-1)  <  ==  RDATA  (I) 

End  Do 

RDATA  (257)  <  =  =  RDATA  (256) 

(Record  RDATA) 

(Fourier  transfora  of  the  L3F  of  Optics  is  calculated.) 
Call  FXFORH 

(Determines  if  plot  of  Fourier  transform  of  Optics 
function  is  desired.) 

If  PLOT_OF_FXFOR!! _  OF_D  PTICS  wanted  then 
Ca  11  Plot  (255  ,  RDAT  A,  2) 

End  If 

(Calculation  of  tae  quotient  of  two  Fourier  transforms.) 
(This  section  takes  the  Fourier  transform  of  the  system 
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and  divides  it,  point  by  point,  by  the  Fourier  tr 


r  a  r  s  r  ora 


of  the  optics.  This  yields  tie  modulation  transf? 
function  (HTF)  of  the  a  tmospi  ere. ) 

Input  RDAIA_OF_FXFORM_OF_S3aRCE 
Do  for  I  =  257  to  512 


RD  ATA  (I)  <  =  =  D 


End  Do 


Do  for  I  =  1  to  256 


RD  ATA  (5  14-1)  £  ==  R  DATA  ( I ) 


End  Do 

Input  R  DAI A_OF_FX F  0RM_0 F_DPI ICS 
Do  for  I  =  1  to  256 

R  D  AT  A  ( I )  <  =  =  IDATA  (I) /RDATA  (I) 

End  Do 

(Determines  if  a  plot  of  the  dlF  is  desired.) 

If  PL3TJ3F_HTF  jDP_ATX3S_vanted  then 
call  Plot  (32,  RDATA  ,  1) 

End  If 

(Calculation  of  Cn2  by  curve  fitting. 

This  procedure  is  described  in  Chapter  III.) 

(All  "R"  variables  used  in  this  calculation  are 
initialized  to  0.) 

R 1 7<  3=  3  18<  =  =  R19<  =  =  R20<==R21<  =  =  R22<==a23<=  =  E2U<==R25<  ==Q 


(Analysis  of  data  shows  than  computer  tima  is  wasted 
beyond  96  points.) 

Do  for  I  =  1  to  95 

(R18  -  when  summed  is  tha  (f(it)l°/3.  This  is  the 
total  of  tha  angular  spatial  frequencies  squared.) 

R  1  9<==  ((I-1|  *Scale_?aotor)  **  (10/3) 

R  1  7<  =  =  R17  ¥  R18 

(R19  -  when  summed  is  the  (r(ii)s/3.  This  is 
the  total  of  the  angular  spatial  frequencies.) 

R  1  9<= =  ( (1-1 1  *Scale_Paotor)  **  (5/3) 

R 2  0<==  Hi 9  *  R20 

(R21  -  is  the  total  number  of  points.) 

R2  1<«  I 

Do  While  RDAT  A  (I(  >  3 

(R22  -  when  summed  is  tha  total  of  the  product  of  the 
natural  logarithm  of  each  point  with  R19.) 

R 22<==  1  n  (RD  ATA  (It  J  *R  19 
R23<  =  =  R22  ¥  R23 

(R24  -  when  summed  is  tha  totaL  of  the  natural  logarithm 
of  each  point.) 

R24<  =  =  l  n  (RD  ATA  (I)  ) 

R25<3=  R  24  ♦  R25 
End  Do 
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d  Do 

(Output  matrix  values  for  deougging. ) 

Output  R17,  R20,  R21  ,  R23,  R25 
(Assign  value  to  matrix. 

This  section  takes  the  above  calculated  values  and 
sets  up  a  matrix.  The  matrix  aquation  solved  is 

B  =  A-i*C.) 

(Create  3  matrices  A  ( i ,  ^ )  ,3(2,1)  ,  C  (.1,1/.) 

A  (1,1)  <=  R20 
A  ( 1 ,  2)  <  —  R  2 1 
A  (2,1)  <—  R  17 
A  (2,  2)  <—  R  20 
IUVERT_  HAT  RIX_A 
C  (1,  1)  <=  R  25 
C  (2,  1)  <  — R23 
B<  ==  A  *  C 

(Output  matrix  values  for  debugging.) 

Output  B  (1  ,  1)  ,B  (1,  2) 

(Calculation  of  Cn*  using  equation  (2.  3) . ) 

R22<==  8(1,  1)  / ( -2 1 .49* Range* R 2  0*  (Wavelength**  (-.  333  3  3)  )  ) 
Output  »CNSQ=',R22 

(The  program  now  starts  the  prediction  phase.  It 
calculates  a  source  assumed  ;  o  have  a  Saussian 
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distr ibu  tion  by  means  of  the  computed  source  pattern 
below.  Using  similar  Fourier  transform  techniques 
the  program  uses  the  calculated  value  of  Cn2  and  predict 
the  power  incident  on  the  target.) 

(Computed  Source  Pa tt er n : A  =  A3 *e xp  (-x**2/2* Sigma** 2) 

(This  is  an  "arbitrary"  Gaussian  source  pattern 
with  a  standard  deviation  for  Sigma  set  equal  tc  2.) 

AO  <==  A  _Z  Z  2  Q 

C<==  S igma*Scaie^P actor 

Do  for  I  «  1  to  25  6 

F <  =  (1-1)  *5cale_Factor 
3<  =  F  **2/  (2  *  C**2 ) 

Do  While  3  >  13 

2DATA  (I)  <==  0 
End  Do 

2D  AT  A  (I)  <==A3  *exp  (-G) 

End  Do 

(Record  RDATA.  This  records  the  computed  source.) 
(Determines  if  plot  of  computed  source  is  desired.) 

If  P LOT _0F_ COMP UT E  D_S3  (JRCE_w  an  ted  then 
Call  Plot(32,  RDATA, 0) 

End  If 

(Calculates  the  line  spread  function  of  the  computed 
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scurca. ) 


Cali  LSF 

(Determines  if  a  plat  of  the  LSF  of  the  computed  scurcs 
is  desired.) 

If  PLD  T_OF_LSF_OF_COMP  _SO(JRC  E  =  1  then 
Ca  11  Plot  (32,  RDAT  A  ,  0) 

End  If 

(Calcuxa  cron  of  rha  Fourrsr  oransrcrw  --  LoF  or  ccutpu  re 
source. ) 

Do  for  I  =  257  to  512 
RDATA  (!)<  =  =  3 
Er.  d  Do 

Do  for  I  =  1  to  256 

RDATA  (514-1)  <==  R  DATA  (I) 

End  Do 

RDATA(2  57<==  RDATA  (256) 

Call  F  X  FORM 

(Determines  if  a  plot  of  the  Fourier  transform  of  the 
LSF  of  the  computed  source  is  desired.) 

If  PL0T_0F_FXF3RJ!_0F_C0MP0rS D_SO(7RCS_wan ted  then 
Call  Plot  (255  ,  RDATA,  2) 

End  If 

(Calculation  of  the  product  of  two  Fourier  transforms.) 
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the  Fourier  transform 
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with  the  dTF  cf  tie  atmosphere.) 
Do  for  I  s  1  to  256 
R  D  AT  A  ( I )  <-  = 


ED  ATA  (I)  _COMP  DTED_SD[JRDE*RDATA  (I)  _FX  FK_OF_MTF 
2n  d  Do 

(Product  of  FXFQRhS  cf  Sour ce* A tmospher e*0 ?t ics . ) 


RDATA (± ) <== 


HD  AT  A ( 1 ) _ PRODUCT _D F_2_ F X FU S *  H D AT A  (1) _FXFM_OF_OPTICS 
Do  for  I  =  2  to  2  56 

RD  ATA  (I)  <=  = 

HD  ATA  (I)  _?RDDUCT_0F_2_FXF;1S*RDATA  (I)  _fxf:i_o?_d?tics 
?.  D  AT  A  (51U-H  <  ==  HD  AT  A  ( II 
End  Do 

RDATA  (257)  <==  RDAT  A  (256) 

(The  above  data  predicts  what  the  Fourier  transform 
of  the  entire  system  is,  using  the  calcuiiced 
value  of  C n2  for  the  atmospheric  turbulence.) 

(Determines  if  a  plot  of  the  result  of  Fourier  transform 
is  desired.) 

If  PLOT  OF  FXFORM  PRODUCTS  wanted  then 


Da  11  Plot  (32,  RDATA, 2) 


(Inverse  Fourier  transform  gives  target  LSF. ) 

Call  IN  VERSE_FXF0  3  M 

(Determines  if  a  plot  of  the  inverse  Fourier  transform 
is  desired.) 

If  PLOT_OF_INVERSE_FXFROM_wanted  then 
Call  Plot  (32,  RD  AT  A ,  0) 

End  If 

(Converts  target's  one  dimensional  IFF  to  a 
two  dimensional  PSF  by  Abel  transform.) 

Call  ABEL 

(Determines  if  a  plot  of  the  ADel  transform  is  desired.) 

If  PLOT_OF_ABELJCFORM_wanted  then 
Ca  11  Plot  (32,  RDAIA,0) 

End  If 

(Calculates  the  fraction  of  power  inside  circle  of  radius 
R.  This  predicts  the  fraction  of  power  that  will  be 
incident  on  target.  ) 

RD  AT  A  ( 1  )  <=  =  0.25*PI*RDAIA(1| 

Do  for  1*2  to  256 

RD  AT  A  (I)<  =  =  2  *PI*  RD  ATA  ( I)  t-RDATA  (1-1) 

End  Do 

(Determines  if  a  plot  of  the  fraction  of  power  is  desired.) 
If  PLOT  OF  POWER  wanted  then 
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Call  Plot  (32,  RDAr  A ,  0) 

End  If 
Ena  Laser 

(The  following  subroutines  are  used  in  the  Algorithm 
Laser .) 

Subroutine  JTr  a  nsf  ?r  (File_N umber) 

(The  transfer  subroutine  gets  caw  data,  both  Dreamble 
and  data  from  the  waveform  massage  sent  by  the 
468.  It  processes  the  message  by  finding  the  minimum 
value  of  the  array  then  subtracting  this  value  from 
each  element  in  the  array.  This  "zeros"  the  array.) 

Output  'Ensure  Egaipnent  set  up  properly' 

Output  'Continue  when  ready' 

S<==  (-1) 

Do  While  E  =  -1 

("DATA"  is  the  I/O  buffer  where  the  data  from  the  468  are 
sent.  The  status  of  the  buffer  is  read  while  data  are  being 
transferred.  When  the  transfer  is  complete,  the  interface 
is  cleared.) 

Data<==  Trans  ferr ed_Dat  a 
If  Data_Tr aas fer_Complete  then 
E<==  0 
End  If 
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2nd  Do 

Clear, Inter  face 
N<==  1 

{This  is  where  header  format  stops  and  data  begin.) 

A<==  po  s(Data, '  S' ) 

Do  for  I  =  A+2  to  687  by  16 
Do  for  J  =  1  to  16 

If  J«-I>5  88  t  hen 
RJ<==  0 
N<  =  =  N+1 
End  If 

BJ<»®  nu  m  (Data  (1  + J)  ) 

(The  numerical  values  of  each  element  are  stored  in  YDATA. 
YDATA  (NT)  <==  HJ 
N  <  ==  SM 
E  nd  Do 

End  Do 

E<==  minimum  (TDArA) 

Do  for  I  =  1  to  5 1  2 

YDATA  (I)  <==  YDATA  (I)  -E 
RDATA  (I)  <==  Y  DAT  A  (I) 

End  Do 

Store_RDAIA  (Fil e_S  umber) 
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Return 


End  Transfer 

(This  subroutine  calculates  the  Line  Spread  Function 
for  the  previously  recorded  point  spread  function 
by  using  equation  (1.1). | 

Subrout ine_LSF 

(Depending  on  which  file  number  has  been  passed  down 
from  the  calling  subroutine,  either  scale  or  source 
data  are  loaded.) 

If  File_Number  =  Scale_Data  then 
Input  Scale_Data 
RD  ATA<==SCALE_DAr  A 

Else 

Input  RDAIA 
End  If 

Do  for  I  =  1  to  51 2 
IDATA  (I)  <  =  =  3 
End  Do 

Do  for  I  *  1  to  256 

IDATA(I)<  =  =  3  DATA  (I) 

End  Do 

(Plots  have  shown  that  computer  time  is  lost  and  no 
valuable  information  is  gained  beyond  about  24  points.) 
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Do  for  I  -  1  to  2« 


J<=  1 
Output  I 
Q<  =  IDATA  (I) 

Do  While  8  <  2a 

R<==  S2RT  (I*I  +  J*J) 

Q<== 

Q  +  2*  ( (1  -  fraotion  (3)  )  *i  (ir.t  (?.)  )  ♦ fraction  (?)  *1  (int  (?)  ♦  1|  | 
J<==  J*1 
End  Do 

ID  ATA  (I)  <  =  =  2 
End  Do 

Do  for  I  »  1  to  21 

RD  ATA  (I)  <  =  =  I  DATA  (I) 

End  Do 

If  File_Nuaiber  =  Scale_Data  then 
Store^Sxral  e_0  at  a 

Else 

Store_RDATA 
End  If 
Return 
End  LSP 
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(This  subroutine  calculates  tha  Fourier  transform 
of  the  given  data  using  the  Cooley-Tukey  Algorithm. 
If  the  inverse  statement  is  true,  then  the  inverse 
Fourier  transform  is  calculated.) 

Subrout ine_FXFOEM 
Set_Had ian_Mode 

(2**9=512  which  is  the  number  of  points.) 

N<==  9 

Do  for  I  =  1  to  51 2 
IDATA  (I)  <==  3 
2nd  Do 

If  File_Number  =  Scale^Data  then 
Input  Scale_Data 

Else 

Input  RDAIA 
End  If 

T<==  PI/2**  (N-1 ) 

Do  for  J  =  3  to  2**  (M-1)  - 1 
(31  is  a  bit  inversion  subroutine.) 

Call  3I(J,P,N-1) 

C<=  cos  (P*D 

If  ISVERSE_FXFORM  =  True  then 
Fl?7<==  1 
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?lg  7<==  0 


End  If 

P<==  sia  (P*D  *(1-2*Flg7) 

Do  for  I  =  2*R0*J*1  to  2*R0* J+RQ 
31  <  -=  3  DATA  (I ) 

R2<==  RDATA(H-RO) 

33  <.  ==  IDATA  (I) 

R'4<==  IDATA  (I  +  RO) 

ROATA  (I) <  =  =  3 1+R2*C*R4*P 
IDATA  (I)  <--  33*R4*C-R2*P 
HDATA  (I  +  B0»  <■*  31-R2*C-R4*P 
IDATA  (1*30)  <  =  =  R3-34*C*32*P 
End  Do 
End  Do 
Output  M 

End  Do 

(This  section  is  for  re-ordering  the  block.) 

Do  for  I  «  0  to  2 *  *N-1 
Call  BI(I,J,N) 


Do  While  I-J  <  or  =  3 
Do  While  I  *  J 

?<==  RDATA  (I*  1)  /SQRT  (2**N) 
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Z<=  =  IDATA  (If  1)  /SQRT  (2**N) 


HDATA  (1  +  1)  <==  RDATA  (Jf  1 1 
IDATA  (Jfl)  <  =  =  IDATA(Jfl) 

RDATA  (Jf  1)  <  =  =  P 
IDATA  (Jfl)  <  =  =  Z 
End  Do 

RDATA  (I»-  1)  <==  RDATA  (1+1)  /SQRT  (2**N) 
IDATA  (If  1)  <=  =  IDATA  (1+1) /SQRT  (2**N) 
End  Do 

If  Pile_Nuaber  =  Scaie_Da-a  then 
Store_So ale_Data  (File_Number) 

E  lse 

S  tors_RD  AT A ( File^M  a  m  her) 

End  If 

End  Do 

Set_Degrees_Mode 
Return 
Ena  FXFORS 

(This  subroutine  fakes  the  binary  number  ?1  conn 
P3  bits  and  inverts  it  end  for  end,  e.g.  310111 
111010.) 

Subroutine  BI(P1,P2,?3) 

?2<3=  0 


aining 

become 
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?U<  =  =  P  1 


Do  for  Z  =  1  t:  ?3 
?  4  <==  P4/2 
?  2  <==  2*P2 

If  fraction(?4)  If  0  than 
?  2<==  ?2  1 
End  If 

P 4  <==  int  (P'4 1 
End  Do 
Return 
End  31 

(This  function  calculates  Airy (x)  =2*J  1  (xi  /z  where  Ji  (j| 
is  the  Bessel  function  of  order  one.) 

Function  AIRY(P1i 
If  PI  <0  then 

Output  •  ERROR -AR3  OtlENI  LESS  THAN  D  ' 

Step 
End  If 

If  PI  =  0  then 
R  4  <«  1 

Return  R4 
End  If 
R5  <==  0 
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If  PI  >  15  then 


R6<==  90  *  PI/2 
R  1  2<  =  =  1.4*?1  ♦  60/PI 
End  If 

If  PI  <  5  then 

R6  <==  20  ♦  1 3  *P 1- P 1 **2/3 
R  1  2<  =  =  6  *  PI 

Else 

R6<==  20  <■  1  0  *P  1- P 1  **2/3 
R  1  2<-~  1 . 4 *P  1  +  50/PI 
End  If 

R1  2<~-  maximum  ( lit  (P12)  ,  int  ( 3  +  P 1/4) ) 
Do  for  :i  =  R 12  to  F.6  ay  3 

r  a  <==  i*9~23 

R 1 3<==R14<==  0 
If  M/2  =  int  (M/2)  then 
Flg10<==  0 

Else 

FlglO  <  =  =  1 
End  If 

Do  for  J  -  1  to  1-2 

R  1 5<==  2  *  ( M- J )  *R8/P1-R13 
R  1  3<  =  =  R  8 
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R  3  <==  a 1 5 


If  M-J-2  =  0  then 
R4<  ==  R  1  5 
End  If 

If  PlglO  =  3  than 
Fig  1 0<=  =  1 

Else 

Fig 10<=  =  0 
Er.d  If 

R14<==  a  14+2*R8*FLg10 
End  Do 

R1  5<*»  2*R 8/P  1-R1  3 
R14<  =  =  R14+R15 
R4  <==  R4/R1  4 

If  (abs  (R4-35) -abs(R'4*1  *e-6)  )  <  or  = 
R4<=  =  2*  R4/P 1 
End  If 
Return  R4 
R  5  <==  R4 
End  Do 

Output  'ACCORACr  (JOT  OBTAINED* 

Return  R4 
End  AIRY 


0  then 
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(The  subroutine  Abel  takes  a  one-dimensional 
line  spread  function  and  calculates  a  two 
dimensional  point  spread  function.) 

Subroutine  ABEL  (File  _  Number) 

N<==  RDATA  (1) 

RDATA  (  1)<==  1.4*RD  ATA(  1)  -1.3*RDATA  (2)  ♦  .  4+EDATA  (3) 
Do  for  I  =  2  to  5* 

M<—  RDATA  (I) 

RDATA  (I)  <  =  =  .  4*N+.  2*M-.  5*  RDATA  (1*1) 

N<=  M 
End  DO 

Do  for  1*1  to  54 

RDATA  (I)  <  =  =  RDATA  (I)  /(2  *SQRT  (  (I  +  .  1)  *  *2-1*1)  ) 
Do  for  J  =  1  to  54 
RDATA  (I)  <=  = 

RDATA  (I)  +RDATA  (J)  /SQRT  (  (J*.  1)  **2-1*1) 

End  Do 

R  D  ATA  (I)  RDATA  (I)  /PI 

Output  I 
End  Do 
Return 
End  ABEL 
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(This  subroutine  prints  out  data  in  tabulae  fora.) 
Subroutine_TABLE(File_Nuaber| 

Input  R DAI  A  (File_S umber) 

Do  for  I  *  1  to  32 

Do  for  J  =  1  to  15 
Output  RDATA  (  16  (1-1)  «-Ji 
Er.d  Do 

Output  RDATA  (16*1) 

Return 
End  TABLE 

(This  subroutine  plots  the  desired  data.) 

(Each  time  a  plot  is  called,  tils  subroutine 
plots  the  particular  graph  and  labels  it. 

The  horizontal  axis  is  labeled  according  to  the  first 
and  last  parameters  passed  by  the  calling  subroutine.) 
S  ubroutine_?LOT  (Pi  ,  Fi 1 e_H  umber,  P  3 ) 

Input  RDAT  A  (File_!1  umber) 

A<==  0 
B<==  ? 1 

C<  ==  minimum  (RDATA  ) 

D<==  ma  ximum  (RDAr A ) 

XM IN <==  -.  1  *  (B-A) 
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XM Aa<=  =  B+. 05* (S-A) 

YMIN<=  =  C- .  1*  (D-C) 

YM  AX<=  *  D+.05*(D-Z) 

E<==  B 
F<==  10 

(Osing  the  value  of  PI,  the  horizontal  increment  fc 
plotting  and  labeling  is  determined.) 

If  PI  =  512  then 
3<  =  64 

Else  If  PI  =  256  then 
3<  =  32 

Else  If  PI  =  64  taen 
3  <  —  8 

Else 

3  <  —  4 
End  If 
Plot  3,C,  1 

(This  lines  and  places  tic  macks  on  the  horizontal  axis.) 
Do  for  I  «  E  to  0  by  -G 

Plot  I*  (B-A)  /  E+  A,  C ,  2 
Plot  I*  (B-A)  /E+A,  C*  ( D  —  3 )  /  150 , 2 
Plot  I  *  (B-  A)  /  E*  A, C  ,  2 
End  Do 
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(This  lir.es  and  places  tic  marts  cn  the  vertical  a:cis.) 

Do  for  I  =  0  to  F 

Plot  A, I* (D-C) /F+C, 2 
Plot  A+  (B-A)/150,t*(D-C)  /F+1,2 
Plot  A,I*(D-C)  /FH  ,2 
End  Do 

Charact  er_Size  1-2,1, .7 ,0 
(This  numerically  libels  the  vercical  axis.) 

Do  for  I  =  F  to  0  by  -1 

Plot  A-.075*(B-A)  , I* (D» I) /F  +  C  ,  1 
Label  I/F 
End  Do 

(This  sets  up  a  character  string  for  labeling  the  vertical 
axis. ) 

If  P3  =  0  then 

ASTRING<==  'DATA  POINTS" 


L<=  1 

Else  If  P3  -  2  then 

ASTRING<«  *  L  INES/MITRD  RADI AN • 

L<  *»  1  /  (2*L| 

End  If 

(This  numerically  libels  the  horizontal  axis.) 
Do  for  I  =  0  to  E  by  3 
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Plot  A  +  (I/E-.  025)  *  (B-A>  ,  C-  .025*  (D-S)  f  1 


Label  I*L 
End  Do 

If  P3  #  0  then 

L<=  Scale.Factoc 
End  If 

(This  labels  the  horizontal  axis.) 
Plot  .4*(B-A)+A,.35*(D-C) +C,  1 
Label  A  SIR  IMG 

Plot  07*  (B-A)  +  A,  .3*(D-C)  ♦2,  1 
Charact  er.Size  =  1.2,1,  .7,90 
ASTSISG<*«  * NORMALIZED  INTENSITY* 
Label  ASTHING 

Charact er_Size  =  .5, 1,1. 5,0 
(This  plots  data.) 

I<==  0 

Plot  I ,  naximum  (RDA  TA  (I)  )  ,1 
Do  for  I  *  1  to  PI 

Plot  1“  1 , EDA T  A  (I) 

End  Do 

If  P3  *  0  then 

Input  PLOT  .LABEL 
ISTRINGC**  PLOT  LABEL 
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End  If 

Plot  .  o*B,.9*D,1 
Label  •  PLOT  OF' 

Plot  .  6*B,.  87*0,1 
Label  ISTRING 
Plot  .  6  *B, . 84*0, 1 
Label  ESTRING,  •  Laser' 
Return 
Er.d  PLOT 
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V.  SUMMARY 


The  work  reports!  is  this  thesis  supports  the  model 
predicted  and  measured  in  Crittenden,  and  others,  for  the 
long  exposure  case  [Ref.  2 ]«  The  vidicon,  in  replacing  the 
mechanical  slit  scanning  system,  shows  no  degrading  of  the 
signal  data.  It  also  demonstrates  that  a  good  approximation 
for  the  point  spread  function  may  be  made  by  recording  a 
single  TV  line  through  a  laser  spot.  This  line  is  then  used 
to  calculate  the  one-dimensional  line  spread  function.  The 
linearity  of  the  video  tape  recorder  is  seen  in  the  results 
of  the  KTF  and  CQ2.  The  463  oscilloscope  is  the  workhorse 

9 

for  the  entire  systea.  It  affectively  displays,  stores, 
digitizes,  and  transfers  data  on  a  real  time  basis. 

The  overall  systea  does  not  have  the  capability  of  the 
system  described  by  Crager.  However,  the  relative 
simplicity  cf  the  structured  programming  technique  coupled 
with  a  lower  equipment  cost  demonstrates  that  comparable 
measurements  and  data  evaluation  using  this  equipment  can  be 
made.  Further  investigation  should  include  use  of  the  Data 
Precision  6000  digital  waver  tm  analyzer  thereby  allowing 
data  to  b-e  sampled  and  analyzed  on  a  near  real  time 
basis. 
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ANIO  A  -  PROGRAM  LISTING 


0: 

1: 

2: 

3: 

4: 

5: 

6  : 
7: 

8  : 
S : 
10; 
11 
12 
13; 
14: 
15 
16; 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 
.50 
.51 
.32 

33 

34 

35 
.56 

37 

38 


"HE-NE: SOURCE  ON  9  , SCALE  ON  7,&  COLLIMATED  ON  5": 
"GA-AS: SOURCE  ON  8  , SCALE  ON  6  .&COLLIMATED  ON  4": 
ent  "SELECT  CODE  FOR  PRINTER" , A ;dev  "print"  ,A 
£mt  l,z,c;tmt  2,£8  .1  fz;rmt  3,t8.1 
aim  15111024]  ,R[512]  ,*[512]  ,CS[2G]  ,D$[20]  ,  ES  f  2Q  ] 
aim  I  [512]  ,AS[32]  ,B$.(2G]  ,A[2,2]  ,fi[2,l]  ,C[2,1] 
bur  "DATA"  ,  1$  ,.5;ina  *  ,1  ,R;0*R*W*V*u  ;  l*L;Dur  "I",  IS, 
beep; asp  "ENTER  LASER  USED";wait  1500 
ent  "HE-NE=1,GA-AS=C"  ,W 
it  a=1;  "HE-NE "«-E$;  9*0;  jrrp  2 
"GA-AS "*£$; 8*0 

beep;asp  ES&  "  SOURCE  DATA  RECORDED?" ; wai t  1500 
ent  "l=YES,Q=NO"  ,V 

it  V=0  ana  W=l;  stg  0  ; ell  *  'TRANSFER  '(C)  ;  Jitp  2 
i£  V.=G  ana  W=G;s£g  0  ;cll  ''TRANSFER  '  (Q) 
beep;ent  "PLOT  OF  SOURCE? , 1* YES , 0“NC"  ,R 
It  R=l;cll  "'PLOT ' (512,0,0) 

beep; asp  E$&  "  SCALE  DATA  RECORDED? "; wait  15C0 
oeep;ent  "1=*ES , Q*NO" , U 

it  U=0  ana  W*l;s£g  0;cll  "TRANSFER  '  (7*0)  ;jn:p  2 
it  U=0  ana  W=0; stg  0;cll  "TRANSFER '( 6*0) 
beep ;asp  "WHAT  IS  SCALE  FACTCF?" ; wait  luOO 
ent  "SCALE  FACTOR*?  ;G  GETS  PLOT"  ,S 
it  S=G  ana  U#0; y-2*0;cll  "PLOT"  (512 ,0, C) 
it  S=G  ana  0=0 ;  ell  "PLOT"  ( 5 12  ,0 , 0 ) 
it  S=0 ;  3mc  .-3 
S*L 

it  w=l;  9*T*0;j3mp  2 
:  8*.T*0 
lat  T , I  $ ,R( *] 
pos ( IS , “R: " ) *B 
val ( I  $ (B+2 ,  B+4] )*rl 
pos  (I$,"  IT  :  ". )  *B 
I  S  (B+.3  ,B+4  ]  *AS 
lot  T-2 , 1  $ ,R [ *  J 
pcs(  I$,"R:  ".)  *B 
val(IS[B+2,B+4])*r2 
pos  ( I S  ,  "  LT :  ".)  *B 
IS(B+3,B+4]*BS 
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39:  laf  T-2,I$,R[*] 

40:  pos ( 1$ , "R: " } ♦B 
41:  val  (1$  f  S+2  ,B+4]  )  -*-r3 
42:  pos  ( 1$ ,  "IT : ".)  ♦B 
43:  I5(B+3,B+4]  -C$ 

44:  if  r2  =  rl  and  A$.=B$  ana  r2  =  r3  ana  BS=C$;jmp  2 

45:  beep;asp  "ERROR  IN  TINE  SCALE" ; stp 

46:  "E" :ldf  Q, 1$ ,R [* ] ; ina  Y;ara  R+Y ;0+A+B+C 

47:  for  1=1  tc  512 

48:  A+Y  [  I  ]  *A  ;  if  I <51;B+Y{ I] *B 

49:  if  I>463;B+Y('l]-*B 

50:  next  I;prt  A  ,B 

51:  B/100-D; A-512*D-A 

52:  tor  1=1  tc  512; Y[I J I J ; C+Y ( I j < 

53:  it  C<A/2;I»J 
54:  next  I;prt  C,J 
55:  for  K=1  tc  256 

5t>:  Y [  ( J.-K+L)  ir.ca512+l]  +  Y(  (K+J-l)  ir.cc512+l]  -I[K] 

57:  next  K;ina  R;ara  I*R 

58:  "IMAGE  POINT  SPREAD  FCN"*I$;rcr  10, I$,R[*j 
59:  beep;ent  "PSF  CU.TPUT  =  1  AND/OR  CCNT"  ,Z 
60:  beep; ent  "PSF  PLGT=1  AND/OR  CONT" ,Y 
61:  if  Z=l;cll  "TABLE"  (10) 

62:  if  Y*l;cli  "PLOT"  ( 22 , 1C  ,1 ) 

63:  ir.in(I  [*]  )*r22;max(I  [*].)  -r23 
64:  prt  r22,c23 

65:  oeep;ent  “LSF  GUTPUT  =  1  AND/OR  CCNT”  ,Z 

66:  beep; ent  "LSF  PLOT=l  AND/CR  CONT”  ,Y 

67:  ell  "LSF';  "IMAGE  LINE  SPREAD  FCN"-I$;rcf  lU,I$,R[*] 

68:  it  Z=l;cll  "TABLE '(10) 

69:  if  Y»l;cll  "PLOT" ( 32 , 10 ,1) 

70:  "CALCULATION  OF  FXFORM  OF  LSF":laf  10 , 1$ ,R( *] ; ina  I 
71:  ara  R-*I 

72:  I ( 1 ] *R ( 1] ; for  1=2  tc  256 
73:  I  [I]  *R[I.)*R(514-I]  ;next  I 
74:  R(256] *R[257] ;rct  1Q,I$,R{*1 

75:  cll  "FXFORM";" FXFORM  OF  IMAGE  LSF" -I $ ; rc f  10,I$,R[*j 
76:  beep ;ent  "FXFORM  OF  LSF  OUTPUT  =  l  AND/CR  CONT"  ,Z 
77:  beep;ent  "FXFORM  OF  LSF  PLOT  =  l  AND/OR  CCNT"  ,Y 
78:  if  Z=l; cll  "TABLE" (10) 

79:  if  Y»l;cll  "PLOT"  (32 ,10,2) 

80:  "B=D I A  OF  CBSCUR/DIA  OF  OBJ  LENS  IN  METERS": 

81:  "B=0 . 064!/.  164 "  : 

82:  " r 1=SCALE  OF  DATA  IN  MICRORADIANS , rl»l. 60" : 

83:  "0=DIA  OF  OBJ  LENS  IN  METERS ;0= 0 . 1 64M" : 

84:  "W=WAVELENGTH  IN  METERS": 
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65:  "CALC  CF  DIFFRACTION  LIMIT  POINT  SPREAD  F CN" : 
fab:  ina  R;it  w=l; b . 328e-7*N ; jnp  2 
87:  9.05e-7-*w 

86:  .  064/.  164-*B;  .  164-*0  ;L-*rl 

89:  B  “2-*D  ;  1-D+H ;  n*rl*0*le- 6/W+Z 

90:  for  L=1  tc  256 

91:  Z*  (I-l)-*Y;it  Y>30;gtc  +3 

92:  ('AIRY'(Y)-D*'AIRY'(Y*B)  )  ~2/H  “2-*R  [  I ] 

93:  I-*K;gto  +2 
94:  0 -R [ I ] 

95:  txa  .5;asp  R(I];next  I;"OPTICS  FUNCTION "♦I  $ 

96:  ref  11 , 1  $ , R { * ] 

97:  beep; ent  "OPTICS  FCN  PLOT=l  AND/OR  CONT"  ,Y 
98:  it  Y=l;cll  ' 'PLOT '( 32 ,11,1 ) 

99:  sig  5;cll  ’ 'LSF  ';  "LSF  OF  OPTICS  FCN " -I  $ ;  ret  11,I$,R[*J 
100:  beep ;  ent  "PLOT  OF  LSF  OF  OPTICS  FCN  =  1  AND/CR  CCNT",Y 
101:  ir  Y  =  l;cll  "PLOT '<32,11, 1) 

102:  sfg  5;  ina  I;ara  R-*I ;  I.[  1]  ♦R[  1] 

103:  for  1=2  to  256 ; I  ( I J  ♦R  ( I J  -*R [  514-1  ]  ;  next  I 
104:  R[25b]  -R[257]  ;rcr  jl,I$,R[*] 

105:  ell  'FXFCRM';  "FXFCRM  CF  LSF  OF  OPTICS"  ♦II? ;  ter  11,I$,R(*] 
106:  beep  ;ent  "PLOT  OF  FXFM  OF  GPTICS  =  1  AND/OR  CONT",Y 
107:  if  Y=l;cll  * 'PLOT  '  ( 256 , 11, 2 ) 

108:  "CALCULATION  CF  QUOTIENT  OF  TNO  FOURIER  TRANSFORMS " : 

109:  lot  10 ,15, Rt*] ; ina  I;ara  R-I 

110:  iaf  ll,I$,Rl*j 

111:  for  *1=1  tc  256 

112:  I[I]/R[L]-R[I]  ; nex t  I 

113:  \MTF  CF  SYSTEM"-*!?;  rcr  10,15  ,*[*] 

114:  beep  ;ent  "PLOT  OF  RTF  OF  SYS=1  AND/OR  CONT",Y 
115:  it  Y=l;cll  ' 'PLOT ' ( 32 , 10 , i) 

116 :  Q-*rl7-*rl6-*rl9-*r20-*r21-*r2  2-*r2  3-*r2  4-*r2  5 

117:  tor  1=1  tc  96 

118:  ( ( I - 1) *L ) “ (10/3) -rl8 

119:  rl8+rl7-*rl7 

120:  ( (I-1)*L) “ (5/3) -rl9 

121:  rl  9+r20-*r  20 

122:  I*r2 1 

122:  if  R ( L] <=0; g to  "CC" 

124:  (In  (R(I]  )  -*r24)  *rl9  +  r  23-*r23 

125:  r24+r25-*r25 

126:  "CC" : next  I 

127:  prt  rl7,r2G,r21,r23,r2  5 

128:  r  2 0 -* A ( 1 , 1  ]  ;r  21-A [  1 .2] ;  rl 7-A [ 2  ,1]  ;  r2 0-A (  2, 2  ] 

129:  r25-*C[l,l]  ;  r23*C.[  2,1] 

130:  inv  A-*A;mat  A*C-*B;flt  5;aprt  B 
131:  "CALCULATION  OF  CNSQ"  : 

132:  B  [1,1]/  (-21.4  9*145*r2  0*h"  (.-.  3  3  33  3)  )  *r2  2 
133:  prt  "CNSQ*"  ,r22;  stp 


91 


124:  cfq  ;laf  2 

125:  "TRANSFER"  :  beep;  asp  "ENSURE  EgUIP.  SET  UR  PROPERLY " 

126:  wait  1500 

127:  asp  "PRESS  CONTINUE  WHEN  READY" ; Stp 
128:  asp  "wait ing “ ; wai t  1500 
13S :  but  "DATA" ; 0*2 
140:  tfr  703, "DATA" 

141:  asp  "transferring" 

142:  rds(  "DATA"  )  +E;  if  E=-l;  Dir.p  0 
143:  clr  703 

144:  asp  "setting  output" 

145:  pos  (IS,"%")*A;1-N 

146:  tor  I=A+2  tc  687  by  16 

147:  for  J  =  1  tc  16;i£  J+I>688; C+rU ; jmp  3 

148:  num <I$[I+J] ) -rJ 

149:  rJ-*Y(N] 

15o:  N  +  1-*N 
151:  next  J 
152:  next  I 

153:  irin  (  Y  ( * ]  )  +E;ror  1  =  1  tc  512; Y [ I ] -E*Y [ I J ; next  I 

154:  ina  R;a ra  Y-*R;  rcr  pi, 15  ,R(*]  ;  ret 

155:  "LSF"  :  if  rlg5;laf  11  ,I§  ,R[  *]  ;  jcp  2 

156:  laf  10,I5,R[*J 

157:  ini.  I;ars  R-*I;tar  L=1  tc  24 

158:  1*0 ; asp  I;LlIJ+C 

159:  v ( I* 1+0*0 ) *R 

16  0:  2*  ( (1-rrc  (R) )  *1  [  int(  F)  1  +trc  <R)  *  L[  ir>c  (R)  +1] )  +u*w 
161:  0  +  1*J ;  it  R<24;  jir.p  -2 
162:  <0*L  [  I  ] 

16  3:  next  I ;  ina  R;ara  I*R; if  rig 5;  ref  11 , 1$  , R (*.]  ;cry  ;ret  pi 
164:  ret  pi 

165:  "FXFCRM" : raa;9*N; ina  I ;  if  rlg5;la£  11 , 1 5  ,  R  [  *  ]  ;  jir.p  2 
166:  laf  10 , 1  $ , R [ * ] 

167:  ti/2“  (14 —  1 )  *T 

168:  for  M»1  tc  N y  2“  (N.-M)  *rQ 

169:  for  0=0  tc  2“  (M-l).-l;  ell  '  '31  '  (0  ,P  ,N-1 ) 

170:  cos  (P*T) *C ; sin(P*T) * ( l-2*t lg 7) *P 
171:  tor  I=2*r0*0+1  tc  2*rG*J+rG 
172:  R(I]*rlyR[I+r0]*r2 
173:  I[I]*r3;I(I+rQ]*r4 

174:  rl+c 2*C+r4*P*R( I] ; r  2+r  4*C -r 2*  P*I [ I ] 

175 :  rl- r 2*C-r4*P*R( 1+  rC] ; r3-r  4*C  +  r2*  P*  I  [  I  +  r0  ] 

176:  next  I; next  J;asp  M;next  to 
177:  for  1=0  to  2“N-l;cll  'fll  '( I  ,J  ,N  ) 

176:  ir  I-J>G;gto  "BE" 

179:  if  1= J ;g to  "INC" 

180:  RlI+lJ/V(2“N)-P;I[I  +  l]/y  (2~N)*Z 
181:  R[0+1]-RII+1];I[0+1]*I[I+U 
182:  P*R[0+1J ; Z*I (0+1) 
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153 

184 

185 

186 
187 
IS  8 
lb  9 

190 

191 

192 
193; 

194 

195 

196 

197 

198 
195 
2  CO 
2ul: 
202 

203 

204 

205 

206 
207 
2  06 
209 
21G 
211 
212 

213 

214 

215 

216 

217 

218 
215 
220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 


"INC"  :R  [  I  +1  ]  /  y  (2~N)*R(I+1]  ;I  f  1  +  1]  /y  <2~M*I  I  H-l  I 
"B£":next  I 

it  rlg5;rcr  11, 18  , R  (  *]  ;aeg  ;  ret 
a  eg ; ret 

"BI "  :  0  ♦p2 ;  p  l-*-p4 
tor  Z=1  tc  p3 
p4/2-*p4  ;2*p2-*-p2 
it  r  rc  (  p4)  *0;  p2+l-»p2 
int  ( p  4 }  -*p4 
nex t  Z ; ret 

"AIRY"  :  if  pl<G ; beep ;asp  "errcr-argunientcu"  ;stp 

it  pl=0;  l-*r4;  ret  r4 

0-*-r5;if  pl>15;3mp  2 

2G+lG*pL-pl~2/3*r6  ;  jir.p  =2 

80+pl/2+r6 

ir  cl<  5;  6+pl*rl2;  jirp  2 

1 . 4  *pl+6  0/pl*r 12 

ir.ax(int  irll)  ,.int  (  3+pl/ 4 } )  +rl2 

tor  N=r  12  tc  r6  by  3  ;le-2o  -»r  8 ;  C*rl3*r  14 

src  lb;  ir  R/  2=int (i-./  2)  ;crg  lo 

ter  J=1  tc  y-2 ; 2* (K-d ) *r 8/pl- rl3*rl5 ; rb-r 13 

rl5*r8;it  M-J-2=G; r 15+r4 

ernf  10;  ri.4+2*  r8*flglC-*r  14;  nex  t  j 

2*rb/pl-r  13*rl5 

rl4+rl5-rl4;  r4/rl4*r4 

it  ats  (r  4-r5).-aos  (r  4*le-6)  <»0;  2*r4/pl*r4;  r  et  r4 
r4-*r5;next  to 

Deep; dsp  "ACCURACY  NOT  GSTAI RED " ; wt l t  150G;ret  r 4 

'"IABLE  "  :  lat  pl,I$fR[*J 

tor  1=1  tc  32;rcr  J=1  tc  15 

wr  t  "pr int. 2"  ,R [ 16 ( I -1) +0 ] ; nex  t  J 

wrt  "print.  3"  ,R[  lbl]  ;  next  I;  ret  C*£ 

"PLOT":  lat  p2,I$,R(*] 

0  -*-A ;  pi  +6 ; m in  (  R  (  *  ]  )  *C ;  nta  x  (R  [  *  ]  )  *C 

SCI  A-.l(B-A)  ,B+i .  C  5  t£-A }  ,C-.1(D-C)  fO  -h  .  G  5  (D-C) 

B*E; 1G*F 
it  pl=512;64»G 
if  pl=256;32-*G 
it  pl=64;8*G 
it  pl=32;4-G 
pit  B,C,1 

tor  I=£  tc  G  by  -G 

pit  I  (B-A)/E+A  ,C  ,2 

pit  I  (B-AJ/E+A  ,C+ (Dl-C)/150,2 

pit  I (B-A) /E+A ,C,2 

next  I 

tor  1=0  tc  F 

pit  A , I ( D-C) /F+C , 2 

pit  A+(B-A) /150, I (D-C)/F+C ,2 
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2  32:  pit  A, I (C-C) /F+C, 2 
233:  next  I;pen 
224:  c £12  1.2,1,.. 7 
235:  fxo  1 

236:  for  I=F  tc  0  ty  -1 

227:  pit  A-.075(B-A) , I (D-C)/F+C ,1 

238:  lbl  I/F 

225:  next  I 

240:  fxc  0 

241:  ir  p3=0; "DATA  P0INTS"-A$ ; 1+L 

242:  if  p3=2;  "LINES/MICRORADIAN"-»A$;  1/(2*L)  -L 

243:  for  1=0  tc  E  by  G 

244:  pit  A+(I/E-.025i)  (B-A)  ,C-.025(D-C)  ,1 
245:  lbl  I*L 

246:  next  I;  if  p3#G  ana  p3*2;  "MCRCRADIANS"*AS 
247:  if  p3=2;S-L 

246:  pit  .4(5-A)+A,.-.G£(D-C)-«:;lbl  A$ 

245:  pit  07  (B-A)  +A  ,  ..3  (D-C )  +C;  csiz  1.2,1, .7,50 

250:  "NORMALIZED  III  TENSITY" -A  $;  It  1  A$ 

251:  csiz  . 5 , 1 , 1 . 5 , C ; O+I 

252:  pit  I ,max ( R[ * ] ) , 1 

253:  for  1=1  tc  pl;plt  1-1 ,R [ I ] ; nex t  I 

254:  csiz  1.2,1,.. 7,0 

255:  if  p2#0;gtc  +2 

256:  beep  ;ent  "PLOT  LABEL?"  ,1$ 

257:  pit  . 65 , . SC , 1; lbl  "PLOT  OF" 

258:  pit  .65, .8/5,1; lbl  1$ 

255:  pit  .  65  , .  34C  ,  1;  itl  ES,"  LASEF" 

26C:  pen;ctg  ;ret  0  +  Y 
*26625 


0:  “COMPUTED  SOURCE  PATTERN"  : 

1:  ina  R;  10C0-»-A;i2*L-*C  ;  tr  k  l;faf  10 

2:  tor  1=1  tc  256 

3:  L* ( L—  3 ) ♦  F 

4:  F“2/(2*C*2)*G 

5:  if  G>13;jmp  2 

6:  A*exp(.-G)  ♦R  [  I]  ;I*K;  jmp  2 

7:  0  -*R  ( I  ] 

8:  next  I 

S :  "COMPUTED  SOURCE" *15 
10:  ref  10 , 1  $  ,  R  [  *] 

11:  beep ;ent  "PLOT  OF  COMPUTED  SOURCE  =  l  AND/CR  CCNT"  ,Y 
12:  ir  Y  =  l; ell  '  'PLOT  '  { 32 , 1 0 ,  Q) 

12:  "CONVERTS  PSf  TO  LSF": 

14:  cll  'LSF' 

15:  "LSF  OF  COMPUTED  SOURCE" *1$ 

16:  ref  10 , 1  $  ,  R  [  *  ] 

17:  beep;er.t  "PLOT  Or  LSF  OF  COMP  SOURCE  =  l  AND/OR  CCNT"./ 
18:  if  Y  =  l;cll  "PLOT' (32,10,0) 

19:  "CALCULATION  Of  FXFORM  OF  LSF  OF  COMPUTED  SOURCE": 

20:  ina  I;ara  R-*I 

21:  I  ( 1  ] *R ( 1J 

22:  ter  1=2  tc  256 

23:  L[I]-R[IJ-R(514-Ij  ,-next  I 

24:  R[25b]-R(257] ;rcr  10, !$.*(*] 

25:  cil  *  'FXFORM' 

26:  "FXFORM  OF  LSF  OF  CALC  SOURCE" ♦IS; rcr  iG,I$,R[*] 

27:  beep  ;ent  "PLOT  OF  FXFM  OF  CALC  SRCE  =  1  AND/CR  CCNT",Y 
28:  if  Y  =  l;cll  PLOT '( 256 ,10 , 2) 

29:  "CALCULATION  OF  PRODUCT  OF  2  FXFOKMiS": 

30:  laf  10 , 1  $ , R [ *)  ;  ina  I;ara  R+I;trx  0 

31:  lar  10,I$,R[*J 

22:  for  1=1  to  256 

33:  R  ( I  ]  *L[I  ]-R(  i]  ;next  I 

34:  ref  10,I$,R[*] 

35:  "TRANSFER  FCN  OF  SOURCE*ATMCSPHERE*OP.TICS"  : 

36:  laf  11 , 1  $ , R ( * J ; ina  I;ara  R*I 
37:  I ( 1 ] *R [ 1 ] *R [ 1 ] 

38:  tor  1=2  tc  256 

29:  R[I] *  I [ I ] *R 1 1 ] *R [514-1] 

40:  next  I 

41:  R(256] -R ( 257 ] 

42:  "FXFM  OF  SGURCE*ATMOS*OPTICS " ~I $ 

43:  ref  10,I$,R{*] 

44:  beep ;ent  "PLOT  OF  FXFM  PRCDUCTS  =  1  AND/OR  CONT"  ,Y 
45:  if  Y=l;cll * 'PLOT ' ( 32 , 10 , 2) 
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46:  "INVERSE  FXFCRM  GIVES  TARGE 1  LSF " : 

47:  srg  7;cil  '  'FXFCRM '  ;cfg  7;"  INVERSE  FXFCRM"* 1$ 

48:  ref  10,I$,R(*] 

49:  ceep  ;enc  "PLOT  OF  IKV  FXFM=1  AND/ OR  CCNT"  ,Y 
50:  if  Y=l?cll  "PLOT  '(32,10,0) 

51:  "CONVERTS  LSF  TO  PSF  BY  ABEL  TRANSFORM": 

52:  cil  'ABEL "ABEL  XFORM  PSF"*I$ 

53:  ret  10,I$,Rl*] ; sfg  0 

54:  beep;ent  "PLOT  OF  ABEL  XF0RM=1  AND/CR  CCNT",Y 
55:  if  Y=1 ;cll  * 'PLOT' (  32 , 10  ,  C )  ;stp 

56:  "CALCULATES  FRACTION  OF  POWER  INSIDE  CIRCLE  OF  RADIUS  R" 
57:  .  25*  n*R [ 1 ] *R [ 1] 

58:  for  1=2  tc  256 

59:  2ir*I*R[I]+R(L-l]*R[I] 

60:  next  I 

61:  "POWER  FRCTN  INSIDE  CIRCLE" -I  $ 
o2:  ref  10,I5,R[*] 

63:  oeep; ent  "PLOT  OF  POWER=I  AND/CR  CCNT" ,Y 
64:  if  Y»l;cil  "PLOT ' (32 ,10,0) 

65:  stp 

66:  "LSF":ir  flg£;luf  11  ,1  $  ,R  [  *  j  ;»3ir:p  2 
67:  Ibf  10,1$, Rl*] 

68:  ina  I;ara  R*I;tor  1=1  tc  24 
6  9:  l*J;cisp  I ;  I  ( I  ]  *0 
70:  y(I*I+J*d)*R 

71 :  2*  (  (1-r  rc  (R) )  *1  [  int  (  R)  ]  +r  rc  (R)  *1.  [  ir.t  tR)  +1]  ) 

72:  J+  1*J  ;  if  R<24;jfr,p  -2 
73:  will] 

74:  next  I;  ina  R;ara  I*R;  if  rlg5,*rcr  1 1 , 1  $  ,  R  [  *  ]  ;  c  rg  ;ret  pi 
75:  ret  pi 

76:  "FXFCRM" : rac ; 9*N ; ina  I ;  ir  flg5;lar  11 , 1  $  ,R  ( *  j  joir.p  2 

77:  lot  lu  ,  I  $  ,  R  [  *  J 

78:  7r/2~(  N-l)  *T 

79:  ter  M=1  to  N; 2~ (N-M) *rG 

30:  ter  0  =  0  tc  2"  (R-l).-l ;  ell  ‘  'BI  '  ( J  ,P  ,N-1) 

81:  ccs  ( P*T) *C ;si n ( P*T)  *(l-2*rlg7)-P 
82:  fer  I=2*r0*d+1  tc  2*rO*J+rO 
83:  R[IJ*rl;R(I+rQ] *r  2 
84:  I [ IJ *r3; L(I+rC] *r4 

85:  rl+r2*C+r4*P*R(I] ; r 3+ r 4*C- r2* P*I (I) 

66 :  r 1- r2*C-r4*P*R [ I+r0] ; r3-r4*C+r 2*P*I { I+r0 ) 

87:  next  I;  next  J  ;asp  M;  next  M 

88:  for  1=0  tc  2“N-l;cll  'Bl'(I,J,N) 

89:  if  l-0>0;gto  "BE" 

90:  If  L=J;gto  "INC" 

91:  R  ( I  +  1]  /y  ( 2  ~N )  *P;I(I+l]/y  (  2~N)  *Z 
92:  RfJ+1]  *R(  1  +  1]  ;I  [J+l]  *I[I+1] 

93:  P*R ( J+l J ; 2*1 (J+lj 

94:  "INC" :R ( 1+1] /V ( 2~N ) *R [ I +1 ] ; I ( 1  +  1] /y ( 2*N) *L ( 1  +  1] 
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95:  "EB" : next  I 

96:  ceg;ret 

97:  "BI"  :G*p2;pl*p4 

98:  for  Z=1  to  p3 

59:  p4/2-*p4;  2*p2+p2 

100:  if  f rc ( p4 ) #0 ; p2+l^p2 

101:  int(p4)+p4 

102:  next  Z;ret 

103:  "ABEL "  :  ina  R;lcf  10  ,1  $  ,R [ *]  ;  R [  1]  -N 
104:  1.4*R[l]-1.8*R[2]-h.4*R(3]-R[l] 

105:  for  1=2  tc  64 
10b:  R  ( I ] *M 

107:  .4*N+.2*M-.6*R[I+1] *R[I] ;M-N;next  I 

108:  for  1=1  tc  64 

109:  Rfl  ]/  (2*y  (  ( I-h.l)  *2-1*1.)  )  -R  [  I  ] 

110:  for  J=I+1  to  64 

111:  R[I]+RlJ]/Y((JHi.l)  "2-1*  L)*R[  I] 

112:  next  J 

J13 :  if  R [ L] < . 0 1 ;0 ♦R( I ] 

114:  R  (I  ]  /  u*R  ( IJ  ;asp  I;next  I 

115:  tor  1=65  tc  512; 0»R (I ] ; next  I; ret 

116:  "PLOT" :lat  p2,I$,R[*] 

117:  0-A;pl-E;min(K(*l ) -C ; max (R [ * ] ) *C 

116:  scl  A-.l(B-A)  ,B-h.G5  <B-A)  ,C-.1(D-C)  ,D+.0S(D-C) 

119:  B-£;10-F;if  f IgO ; pr t  A ,B,C,C 

120:  it  pl  =  512;  6  4-*G 

121:  if  pl=256;  32-*G 

122:  if  pl=6  4;  b-*G 

123:  if  pl=32;4-*G 

124:  pit  B,C,1 

125:  for  L=E  tc  0  by  -G 

126:  pit  I (B-A) /E+A  ,C  ,2 

127:  pit  I  (E-A)  /E+A  ,C+ (D.-C )  /  1 50 , 2 

123:  pit  I(3-A)/E+A,C,2 

129:  next  I 

130:  for  1=0  tc  F 

131:  pit  A, I (D-C)/F+C,2 

132:  pit  A+(B-A)/150,I (D-C)/  '■’+€, 2 

133:  pit  A, I(D-C)/F+C,2 

134:  next  I; pen 

135:  csiz  1 . 2, 1, .7 

136:  txa  1 

137:  for  I=F  to  0  by  -1 

138:  pit  A-.G75(n-A) ,1 (D-CJ/F+C , 1 

139:  lbl  I/F 

140:  next  I 

141:  fxa  0 

142:  if  p3=2;  "LINES/MICRCRAD IAN " *A$ ; 1/ ( 2*L, -L 
143:  for  1=0  tc  E  by  G 


97 


144 

pic 

A-t-  ( 1/ E- .  G25)  (8 

-A) 

,C- .025{D 

-C) 

14  5 

lbl 

I*L 

l4o 

next 

I; if  p3#2; "MICRORADIANS " -A$ 

147 

ir  p 

3*2;  S*L 

148 

pit 

.  4  ( B-A)  +A„-.G5 

(D- 

C)  +C ; 

lbl 

A$ 

149 

pit 

-.07 (B-A) +A , -3 

(D- 

C)+C; 

csiz 

1. 

150 

“NORMALIZED  INTENSITY 

“  *A$; 

ibl 

A$ 

151 

csiz 

. 5,1,1.5,0;0* 

I 

152 

pit 

I ,  ma  x  ( R  { *  J  )  ,1 

153 

for 

1=1  to  pi; pit 

1-1 

,R(I1 

;  nex 

t  I 

154 

csiz 

1.2,1,..  7,0 

155 

pit 

. 6£ , . 9D , 1 ; lol 

"PLOT  OF 

M 

156 

pit 

.68 , . 87D ,1; lbl 

15 

157 

pit 

.68, .840, l;lbl 

E$ 

,"  LASER" 

156 

"RET 

":pen;ctg  ;ret 

0- 

Y 

*  175  24 


,.7,  SC 
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